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In writing the following pages I have had two objects in 
view : first, to present the fundamental principles and pro- 
cesses of the Method of Least Squares in so plain a manner, 
and to illustrate their application by such simple and practical 
examples, as to render it accessible to Civil Engineers who 
have not had the benefit of extended mathematical training; 
and secondly, to give an elementary exposition of the theory 
which would be adapted to the needs of a large and con- 
stantly increasing class of students. 

Hence the work falls into two parts, the first practical 
and the second theoretical, but each illustrating and sup- 
plementing the other. The numbering of the articles renders 
reference from one to the other easy ; and the more thorough 
acquaintance the engineer makes with the second part the 
better will he adjust his observations, while it is only after 
much exercise with practical problems that the student can 
become thoroughly grounded in the theory. 

Should the book, then, be taken up by students un- 
familiar with the subject, let me suggest to them, that even 
if their aim be only to acquire a knowledge of its theory, 
the shortest and best way to do it is to become first familiar 
with the practical applications of Fart I.; this attained, the 
rest follows naturally and easily. 
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IV PREFACE. 

As I have not written fot mathematical experts, they 
will doubtless find considerable in the book at which to 
grumble- The idea of mean error does not appear. The 
term *' equations of condition" has been, in accordance with 
the sensible German practice, divided into "observation 
equations** and "conditional equations** {Beobachtungsglei- 
ckungen and Bedingungsgleichungen), and each is used in its 
proper place. Gauss' development of the law of probability 
of error has been followed as the best adapted to an element- 
ary presentation, and if this be objected to as defective, 
I claim at least the credit of knowing and of pointing out 
(Art, iyQ) just what and where those defects are. 

In preparing these pages I have consulted and freely 
used all the works upon the subject within my reach. 
The list of Literature and the historical notice at the end 
of the book will be of interest and value to all. 

If this little elementary work should meet with a 
favourable reception from the scientific public, it may be 
followed by another containing extended applications of the 
method to higher geodetic surveying, and to numerous other 
]>robIenis arising in physical science, which have here been 
necessarily left unnoticed. 

MANSFIELD MERRIMAN. 

-T> BciENTiFio School, 
KiiM' Haven, Conn., U.S.A., Feb. 6, 1876. 
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PART I. 

THE ADJUSTMENT AND COMPARISON OF ENGINEERING 
OBSERVATIONS, 

CHAPTER I. 

Il^RODTJCTION. 

Art. 1. When a quantity is observed with a view to 
determining its magnitude a number is obtained as the result 
of the operation. This number expresses how many units 
and parts of a unit are, according to the measurement, con- 
tained in the observed quantity, and is hence a measure of 
its magnitude. The word observation will be used in this 
book to express such numerical measures as well as the 
operation by which they are deduced. 

Every engineer is cognizant of the fact, that, when several 
observations or sets of observations are made to determine 
the magnitude of a quantity (for example, the length of a 
line) the results do not agree. Since the quantity can have 
only one value all of these results cannot be correct, and 
each one of them can be regarded only as a.n approximation 
to the truth. The absolutely true value of the quantity in 
question, we can never obtain or at least be never sure that 
we have obtained, and instead of it we must accept and use 
a value, derived from the combination of our observations, 
which may perhaps not exactly agree with any one of them, 
but which however shall be the most probable value (Art. 9). 

Part I. of this work will be devoted to the presentation 
and illustration of the methods now in use alnong scientific 
men for the adjustment and comparison of observations, and 
the determination of the most probable values of observed 

f s^ jjl^ Digitized by VjC ^ 



2 INTBODUCTION. 

quantities, with sufficient mathematical reasoning to ex- 
emplify the main principles of these methods in a clear light 
to the general reader. Part II. will contain a full develop- 
ment and discussion of the theories upon which they are 
founded. The first is designed more particularly for the use 
of those "who [have little time to devote to the niceties of 
mathematical theory, but who are desirous of learning the 
fundamental principles and practical applications of the 
science; the second is intended for the use of students 
taking up the subject from a theoretical point of view. Each 
part, however, will be to some extent dependent upon the 
other, and hence for convenience of reference the correspond- 
ing articles and formulae will be marked by corresponding 
numbers. 

Errors of Observations. 

2. CoTtstant Errors are those which always under the 
same circumstances have the same value, and whidh there- 
fore strictly speaking are not enrors, but the results of law. 
As such we may mention : theoretical erroi*s like the effects 
of refraction in increasing the size of a vertical angle, c^ 
the effects of temperature upon the length of rods used in 
measuring a base line,, which effects when their causes are 
understood may be computed beforehand, and are hence no 
longer to be classed as errors ; instrumental like those arising 
from an incorrect graduation of the limb of a theodolite, or 
an imperfect adjustment of the line of collimation in a tele- 
scope, which may be also removed by calculation or by a 
proper mode of using the instrument; and personal errors 
which arise in very delicate observations and are due to 
the habits of the observer, who may for instance in reading 
a vernier always give the number of records too large or too 
small by a constant quantity, and which may be corrected by 
the application of a "personal equation-" 

Such errors being capable of elimination need no furliier 
consideration in the discussion of our subject 

8. Mistakes are a class of errors committed by inex- 

oerienced and occasionally even by the most skilled observers, 

from mental confusion. As such are; mistakes in 
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measuring atl angle arising from sighting at the wrong 
signal ; mistakes in reading an angle by noting 54^ instead 
of 46^, etc. This class of errors sometimes admits of cor- 
rection by comparison with other sets of observations : it will 
receive no-further notice in this book. 

4. Accidental Errors are those that still remain after all 
constant errors and mistakes have been carefully investigated 
and eliminated from the numerical results* Such, for ex- 
ample, are the errors in levelling arising from sudden ex- 
pansions and contractions of the instrument or from the 
effects of the wind, or those often observed in sighting acro^ 
a river arising from the anomalous and changing refraction 
of the atmosphere. More than all, however, such errors arise 
from the imperfections of our touch and sight, which render 
it impossible for us to handle our instruments delicately, or 
estimate accurately small divisions of their graduation. These 
are the errors which appear in our numerical results, however 
carefully the measurements be made, and which form the 
subject of the following pages. Although at first sight it 
would seem that such iiTegular errors could not come within 
the province of mathematical investigation, it will be seen 
in the sequel that they are governed by a wonderful and very 
precise law, viz. the law of probability (Art. 11). 

Principles of Probability, 

6. We must therefore, as preliminary to our subject, 
state and exemplify the mathematical definition of probability 
and of the words "most probable" which we have used in 
Artl. 

If a coin be tossed up into the air, it is said in common 
language that the chances are equal that it will turn up 
hesid or tail, or that the occurrence of head or tail is equally 
probable. So if a die marked in the usual manner be 
thrown, it is said that the odds are one to five in favour of 
throwing the ace, or five to one against throwing the ace. 
This is the statement in popular language of the idea of 
probability; its mathematical expression is but slightly 
different. In throwing the coin we recognise that there are 
two possible cases, either head or tail may turn up, and one 
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4 INTKODUCTION. 

is as likely to occur as the other ; and hence we express the 

pi'obability of throwing a head by the fraction ^ , and the 

probability of throwing a tail also by ^ . So with the die 

there are six equally probable cases, one of which may be 
the ace, and hence the probability of throwing the ace in 

1 5 

one trial is ^ , and the probability of not throwing it is 77 . 

The probability of the occurrence of an event is, then, always 
a fraction whose denominator denotes the whole number of 
possible ca^es (each supposed equally probable) and whose 
numerator denotes the number of cases favourable to its occur- 
rence. 

Thus if there be in a bag 30 red and 20 white balls, and 

a ball is drawn at random, it will be either red or white ; 

both events are equally possible but not eqvxilly probable; 

.30 
for the probability of drawing a red ball is -^ and that of 

20 3 2 

drawing a white ^ , or ^ and - respectively. 

The mathematical expression of the probability of the 
occurrence of an event is, then, a numerical measure of our 
degree of confidence that it will occur. As a fraction may 
have any value from to 1, so our confidence may range 

from impossibility to certainty ; a small fraction like tttoT 

1023 
denotes a very small probability, and a large one like ^7^57 

denotes a very large probability or almost a certainty. 

6. If there be 18 red balls in a bag, and one is drawn 

18 
but, the probability that it is red is to = 1, that is, it is 

certain that such a one will be. drawn. Therefore : Unity is 
the mathematical symbol for certainty. Hence if the proba- 
bility of the happening of an event be known, the probability 
that it will not happen is unity minus the first probability. 
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•I 

Thus tlie probability of throwing an ace in one trial is ^ , 

6 

15 

the probability of not throwing it is 1 — 7? = t; , as shown 

o D 

otherwise in the preceding article. 

7. If there be in a bag 20 red, 16 white, and 14 black 
balls, and one is to be drawn out, the probability that it will 

be red is ^jr , that it will be white ^, , and that it will be 

14 
black -^. If however we ask the probability of drawing 

either a red or white ball, we have 36 favorable cases out of 

the 50 total cases, and the answer is ^7; or ^^^ + 777: * Hence, 

50 50 50 

if an event may happen in different independent ways, the 
probability of its happening is the sum of the separate proba- 
bilities. Thus in tossing a coin, the probability of throwing 

either head or tail is - -f - = 1, that is, one is certain to be 

thrown. 

8. Let there be two bags, one of which contains 7 black 

and 9 white balls, and the other 4 black and 11 white balls. 

The probability of drawing a black ball from the first bag is 

7 . 4 

:r^ , that of drawing one from the second ^z-^ . What now is 
16 ° lo 

the probability if I draw from both bags at the same time 

that both balls drawn will be black ? Since each ball in the 

first bag may form a pair with each one in the second, there 

are 16 x 15 possible ways of drawing two balls; and since 

each of the 7 black balls may be associated with each of the 

4 black balls to form a pair, there are 7x4 cases favorable 

to drawing two black balls. The required probability is 

7x4 . . . 74 

hence ^tfi — tp and since this is equal to =-^ x ^r-;:; we have the 
16x15 ^ 16 15 

principle that the probability of the happening of several in- 

dependent everds is equal to the product of their respective 

probabilities. 
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6 INTRODUCTION. 

Thus if three dice be thrown, the probability that all will 

be aces is 7; x ^ >< ^ = stt^ > * small fraction. 
000 zlo 

9. Suppose two coins to be thrown up at the same time, 
they may both turn up heads or both tails, or one may be a 
head and the other a tail. We wish to determine the re- 
spective probabilities. Let us call the two coins A and B, 
then the cases which ^lay happen are 

^1. A head and B head, 

2. A head and B tail, 

3. A tail and B head, 

4. A tail and B tail j 

and each of these cases being equally likely to occur, has for 

its probability t- Hence the probability that both will be 

heads is -. , that one will be head the other tail 7 + 7 = 0, 
4 4 4 2 

and that both will be tails 7. The sum of the probabilities 

7 + 5 + 7 is unity, as should be the case (Art. 6), since one of 
the events is certain to occur. 

In like manner if there be ten coins thrown up at the 
same time there may occur the following groups, having the 
respective probabilities as annexed. (For an easy method of 
computing these probabilities for 10 or any number of coins, 
see Part II., Art. 9.) 

All 10 coins may be heads, ....^» probability tk^ 
9 may be heads and 1 tail, .........^.. 



1024 

1024 
120 



10:^4* 
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6 may be heads and 4 tails, probabiKtj 
5 S 


210 

1024' 

252 


4 




■ 1024' 
210 


3 , 


...7 


' 1024 ' 
120 


2 


..8 , 

..9 


' 1024 ' 
45 


1 

All may be tail*.. 


'1024' 

10 
■ 1024 ' 

1 

> 



The sum of all these probabilities is unity, since one of 
the groups is certain to occur. In speaking of this case in 
popular language, it is said that it is very improbable that 
all the coins will be heads, in mathematical language we 

say that the probability is small, ^^-jK^y so both common 

and mathematical reasoning recognise that the group of 
5 heads and 5 tails is the most probable. 

In the theory of probabilities, therefore, we say, an event 
is the most probable when the a priori {or theoretical) proba- 
bility of its happening is greatest among all the probabilities of 
all the possible events, that is, when that probability/ is a max- 
imum. Thus if two persons of exactly equal skill should 
play 51 games of cards^ it is possible but not probable that 
one of them will win all the games ; the most probable case ia 
as every one will recognise, that one will win 25 and the 
other 26 games. 

We advise the reader, before proceeding further, to apply 
the above principles of probability to the solution of the 
following problems. • 

Problems, 1. A lotteiy has 250 tickets, 25 of which are 
prizes, and the remainder blanks. What is the probability 
of not drawing a prize ? 
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8 INTRODXTCnON'. 

2. A bag contains 100 balls, 50 are red, 12 are white, 
and the rest are black. What is the probability of drawing 
either a white or a black ball in one trial? 

3. The probability of throwing an ace with a single die in 
two trials is ^. What is the probability of not throwing it? 

4. The probability that A can solve a certain problem 

.2 . . 1 

is K, the probability that B can solve it is 7. What is the 

3 

probability that it will be solved if both try ? Ans. - . 

TO 

5. A bag contains 3 red, 4 white, and 5 black balls. 
Required the probability of drawing 2 red balls, the ball first 

drawn not being replaced before the second trial. Ans, ^ . 

6. Six coins are tossed up. What is the probability that 

15 
4 of them will be heads and 2 tails ? Ans. 7-: . 

64 

7* A coin is tossed Up four times in succession. What 
is the probability that all will be heads ? 

The ProbaJbility Curve. 

10. We may now proceed to consider the law of the 
probability of errors of observation. 

If a person accustomed to the use of the rifle, shoot 500 
times at a target, all the balls will not hit the central bulFs- 
eye, and some perhaps will not even hit the target. The 
deviation of each bullet from the centre of the target is an 
err(yry and furthermore an accidental error (Art. 4) produced 
by changes in the wind, imperfections in the aim of the 
marksman, etc. ; for all constant errors, such as the effect of 
gravitation, are assumed to be eliminated, in the sighting of 
the rifle. An examination of the bullet marks on the target 
shows us, however, that these errors are arranged around the 
central point in a very regular and symmetrical manner, 
''^'irst we observe that small errors (that is, deviations from 
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the centre) are more frequent than large ones ; secondly that 
they are arranged symmetrically around the centre so that at 
equal distances, above, below, and in every direction from that 
centre, the number of marks in a square inch is the same ; 
and thirdly we recognise the fact that very large errors, for 
instance, a deviation from the centre of half a mile, do not 
occur. Further, we observe that the greater the skill of the 
marksman, the nearer are the marks to the centre point. 

Again, suppose an engineer to measure an angle a hundred 
times, each time with equal care. Tl^e readings will disagree, 
the dififerences between the true valmB of the angle and each 
of his results will he an error ; and we recpgnise that these 
errors, like those of the marksman, will be subject to three 
laws ; viz. 1st, small errors will be niore frequent than large 
ones ; 2nd, errors of excess and deftciency (that is, results 
greater and less than the true value) will be equally numerous, 
and 3rd, large errors like those of 2" do not occur. Further, 
the more skilled the engineer in the measurement of angles, 
the nearer will be his readings to the true value of the angle, 
and the smaller will be his errors* 

These three axioms or fund^amental laws form the founda- 
tion of the science of the adjustment of observations. 

11. In any set of carefully made observations, then, the 
probability (Art. 5) of ^ small error is greater than that of a 
large one, the probability of an error in excess of the true 
value is the same as that of an equal error in deficiency, and 
the probability of a very large error is ziero. Thus the pro- 
bability of an error is Si function of the error, so that if x repre- 
sent any error, and y its probability, we may express the 
relation between these quantities by the equation 

(1) y=fi^), 

which is read, y equals a function of pc, that is, y is dependent 
upon X for its value. 

If we take then y as an ordinate, and x as an abscissa, we 
may regard this as the equation of a curve, which must be of 
a form so as to agree with the three laws adopted above; 
viz., its maximum ordinate OA must correspond to the error 0; 
it must be symmetrical with respect to the axis of Y, since 
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equal positiye and negative errors are equally probable, as x 
increases numerically the value of y must rapidly decrease. 




and when x becomes very large y miMt be zero. Fig. 1 repre- 
sents such a curve, OP and OM being errors, and PB and 
MG their corresponding probabilitiea. Further, since dif- 
ferent observations have different degrees of accuracy, each 
set of observations will have a distinct curve of its own, 
whose particular dimensions will depend upon the precision 
of the measurements, but l^hose general form must be that 
QfFifrl, 

In Part IT. we Lave given GaVSS* method of determining 
the equation of this curve frotn the principles of the theory 
of probabilities (Arte. 5, 9), the fundamental axioms of Art. 10, 
and the rule of the arithmetical mean or average (Art 23). 
Its equation is 

or as it is usually writtcfn 

(12) y = e^*V 

in which e denotes the base of the Napierian system of log- 
arithms or 2*71828, and c and h are constants which depend 
upon the accuracy of the observations, and whidi determine 
the dimensions of the curve for any particular case. Regard- 
ing for the present both c and h as unity, and attributing 
values to x, we may by logarithms find the cprresponding 
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values of y. The following table gives a few of these values, 
from which the reader may construct a curve to scale and 
observe its correspondence in form with that given in Fig. 1. 





y^r'' 


1 




OR 


» 


X 


y 


^ 


l-OOOQ 


*l-8 


0-0392 


*0-2 


0-9608 


*2-0 


0-0183 


^0-4 
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This is called the Probability Curve^ and its equation is 
the mathematical expression of the law of the probability of 
error, upon which the whole science of the adjustment and 
comparison of observations^ as here developed, is founded. 

12, This curve has many curious properties, only one of 
which we can notice here. 

If we measure an angle 500 timea there will probably 
be 250 readings greater than the true value and 250 less. 
Subtracting each reading from the true value we have the 
errors, 250 being positive and 250 negative. Now if we lay off 
each of these errors as an abscissa x, and erect its correspond^ 
ing ordinate y, we shall have a curve of the form of Fig. 1,. 
Now the property of the curve to which we wish to call 
attention is this : If we represent the total area between the 
curve and the axis of X by 500, the area OABP on the 
right of the axis Y denotes the number of positive errors 
less than OPy the area PBCM denotes the number greater 
than OP and less than Oif, and the area CMX denotes the 
number greater than OM: likewise the area OABP on the 
left of the axis Y denotes the number of negative errors 
less than OP, and so on ; also, if OP^OP, the area PBABP 
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exhibits the number of errors less numerically than OP, In • 
Fig. 1 we have drawn the ordinates so that the areas OABP, 
PBGM, and MGX represent 125, 45, and 80 errors respec- 
tively. 



Probability qj Errors, (22) F'= -^ J e"*"^ d (hx). 
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If, then, we express the total area of the curve by unity, 
the areas PBABP, MGAGM^ etc. corresponding to the sue- 
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cessive values OP, OM, etc. will be fractions, proportional to 
the number of errors less than those values of a?. The preced- 
ing table, which may be computed by the methods of Part II., 
is taken from the Berliner Jahrbuch for 1834, and gives the 
areas on both sides of the axis y corresponding to successive 
numerical values of hx, x being the error, and h the measure 
of precision, referred to in Art. 11. 

13. To use this table it is necessary to know the value 
of the constant h. Methods will be hereafter given (Arts. 
24, 38) by which its value may be found for any given ob- 
servations. Granting for the present that it may be deter- 
mined, the following example will show the use of the table, 
and exemplify the accordance of theory and practice. 

Bessel in his Fundamenta Astronomiae discusses 470 ob- 
servations made by Bradley upon the right ascensions of the 
stars Sirius and Altair, and determines the measure of pre- 
cision to be A = ^y7:^^^. Now let it be required to find the 

number of errors less than 0"*2, 0"'4, 0"'6; etc. For the 
number less than 0"'2 we must take x = 0"*2, and 

and for the number less than 0"'4, 0"*6, etc., hx will be the 
successive multiples of 0*361 6 by 2, 3, etc. We find then 
from the table, 

for X ^ 0"-2 with hx = 0-362 the area P' == 0-39102, 

... a; = 0:4 Aa? = 0-723 P' = 069372, 

...^ = 0-6 ^ = 1085 P' = 0-87511, 

...aj = 0-8 hx=:l'U7 P' = 0-95926, 

.. x^VO ^^7=1-808 P' = 0-98946, 

. ...x = co hx = oo P' = 100000. 

Now from the property of the curve above noticed tbese 
values of P' are proportional to the number of errors less 
than the corresponding value of x. Hence multiplying each 
of these numbers by 470, the total number of errors, and 
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comparing them with the actual number of errors as given 
by Bessel, we have 

theoretical no. of errors less than 0""2 is 184, actual no. was 1S% 

0-4..* 326,.. »...i....318^ 

0-6. ..412, , ..405, 

0-8. ..451, 445, 

...r0..*465, 462, 

over 1*0 ... 5, 8. 

The agreement between theory and experience though 
very close is not exact, partly, perhaps, because what we have 
called the actual number of errors were computed not from 
the true values of the right ascensions of the stars but from 
their most probable values as deduced from the 470 obser- 
vations. 

The number of errors between any two given limits is 
of course found by simple subtraction; thus, in the above 
example there are 142 theoretical errors between 0"'2 and 
0"-4, 86 between 0''-4 and 0"-6, 39 between 0"-6 and 0"-8, and 
14 between 0"-8 and r-0. 

Problems, 1. A marksman shoots 500 times at a target. 
If his skill is such that when the errors x are measured in 
feet h is unity, what are the number of bullet-marks between 
two circles described from the centre whose radii are 1 and 
2 feet ? Am, 76. 

2. An angle is measured 100 times, and the value of 

the measure of precision found to be A = t7/. How many 

errors are greater than 2" and less than 4"? How many 
greater than 8" ? 

This subject is further considered in Art. 55. 

The Method of Least Squares. 

14. Suppose now a number of observations be made to 
determine the values of certain quantities, for example, the 
lengths of lines. However carefully the measurements be 
made there will be disagreement in the results, and hence 
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we can never be sure that any adjustment that we may 
make will give us the absolutely true values of the lines. 
The most we can do is to determine approximate tesults, 
which shall be the most prdbahU values (Art. 9), and more- 
over be rendered the most probable value by the existence 
erf the observatioDfi themselves. 

The preceding principles of the probability of error atfbrd 
us a general rule for the determination of the most probable 
values of observed quantities. To deduce it in its simplest 
form, let us suppose that the observations are equally good, 
that is, that they are made with exactly the same precii»on ; 
then the constants c and k (Art. 11) which measure that 
precision will be the same for all observations. Let us 
designate by s, t, w, etc. the quantities whose values are to 
be derived from the measurements, and by M^y M^, J/,, etc. 
the results of those measurements, which are made either 
directly (Art. 18) upon s, t, it, etc., or indirectly (Art. 19) upon 
other quantities related to them. If our measurements were 
perfect, M^y M^, etc would be absolutely true values, and 
there would be no discordance in the results found for s, t, ii, 
etc. ; but being imperfect the results do not agree, and hence 
errors (Arts. 4, 10) exist. Let those errors, which are the 
differences between the measurements Jf^, M^, M^, etc., and 
the corresponding trv^ values be denoted by x^, x^, cc^, etc. 

Now if a? be any error and y its corresponding probability, 
the law of the probability of errors (Art 11) gives us 
(12) y = ce-A"-^. 
Hence for the errors actually committed we hav6 
Prob. of the error w^^y^^ te^^^^\ 

a*, =yj = ce"'^*^>*, 

— x^^y^^ce-^^'^, etc. 

Now by Art. 8 the probability of committing all these 
errors is the product of these respective probabilities : hence 
it is 

yiViVz ^*c* = ce-*"^* X ce"*^«' X ce"^^^*^ x etc. 

Designating this product by P, the expression becomes 

(23) P^cxcxcK etc. x e-**(^"+^«'+^>*+***+etc.), 
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in which c, e and h are constants. Now what the true values 
of 8y t, u, etc. are, we cannot hope to find, and must hence be 
content with determining their most probable values ; and the 
most probable system of values is that which has the greatest 
probability (Art. 9). Each of the errors, a?j, a?^* ^8> ^^> ^s 
dependent upon (is a function of) the quantities s, t, u, and 
the most probable system of values for the latter corresponds 
to the most probable system of errors. 

The most probable system, then, is that for which P is a 
maximum (Art. 9), and P in the above expression will be a 
maximum when the exponent of 6 is a minimum^ that is, when 

(24) ' x^ + x^ + a*,* + x^ + etc. = a minimum, 

and as these terms are the squares of the differences or 
errors, we have the principle that, the most probable values of 
quantities, which are the object of measurement, are those 
which render the sum of the squares of the errors a minimum, 

15. To illustrate the application of this principle, let us 
suppose that the measurements ilfj, M^y J^, etc. are made 
directly upon the same quantity M whose true value is z. 
Then we have committed the errors 

(«-Jlf,),(«-if^,(^-iQ.etc. 

Now by the preceding principle the most probable value for 
z is that which renders 

(24) (^-JlfJ>+(;,-JlfJ«+(^-ifJVetc. 

a minimum. If n be the number of observations, the last 
term of this expression will be {z — MJ\ Applying the usual 
method for determining maxima and minima we differentiate 
the expression thus: 

2{z'-M;jdz + 2{z'-M;jdz'h2{z--M;jdz + ... + 2{z'^M,)dz; 

place it equal to zero and divide the equation by 2dz, givin<» 

{z^M,) + {z^M,) + (^- Jf^) + ... + (^-ifj = 0. 

Solving this equation and denoting the resulting value 
of -s^by z^y we have 

^ ^ ' * 8 n 

*o '■ > 

n 
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that is, the most probable value of z is the average of the n 
measurements (Art. 23)*. 

From the principle of Art. 14 arises the name ''Least 
Squares " to express the method now universally in usq for 
the adjustment of observations. 

The Comparison of Observations. 

16. In Art. 10 we observed that the greater the skill of 
the marksman the nearer are the bullets to the centre of the 
target, and also that the more skilled the engineer, the less 
will be the deviations of his readings from the true value of 
the angle. In Art. 11 our law of the probability of error, 

(12) y = ce-*'^, 

contained two constants, c and h, dependent upon the pre- 
cision of the observations. The accuracy of different sets of 
observations may then be compared by comparing their mea- 
sures of precision given by h, or by comparing other constants 
related to h. The one usually employed for this purpose is 
called the probable error r, which we define as the error which 
has such a valvs that the mirriber of errors greater thorn, it is 
the same as the number less than it Thus in the probability 
curve, Fig. 1, the total area XBABX denotes the total num- 
ber of errors in a set of observations (Art. 12), and if the area 
PBABP be one-half of that total area, one-half of the errors 
will be less than OPy and hence by our definition OP is the 
probable error. Now since r is of such a value that if r = OP, 
the area PBABA is \ when the total area is unity, we have 
onljr to find from the table (Art. 13) the value oihx^hr for 
which P'=0*5. By interpolation between the values Aa?=0'46 
and Aa;=0-48, we find that JP'^O-S when hoc = 0-4769, Hence 
we have 

(29) hr = 0-4769, or r = ^-—^ , 

♦ The aboTe is not strictly a demonstration of the law of the average or 
arithmetical mean, for the equation of the probability onrye (Art. 11) is 
deduced npon the assumption that in direct observations the average is the 
most probable true value ; and hence the above is reaUy reasoning in a 
circle. It serves, however, to illustrate to the reader that the principle of 
least squares is in agreement with the universally adopted method of taking 
the arithmetical mean of equally good observations. See Art. 66. 
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from which we see that the values of h land r are alwtiyft 
reciprocally proportionaL 

To render more definite otir conception of the measure ctf 
precision h and the probable error r, let us consider the case 
of two sets of observations made with diflferent degrees of 
accuracy. Let the measure of precision of the first be Ai, and 
of the second h^) then for the probability of errors in the 
first set we shall have a curve whose equation is 

(19) y«A,t7r-i^"W 

and for the second li Curve given by 

(19) y«A,i7r-i€-^n 

in which % is a very, very small constant, independent of h, 
and TT is the number S'1416. Now let us suppose that the 
second set is twice as accurate as the first, so that A^ » h and 
A, » 2A ; tiien the equations will be 

y«feiiV"ie~*^ and y = 2At7r*i€"**^, 

If for purposes of comparison we consider h and i as unity 
and compute the values of y corresponding to successive 

Fig. ^. 
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numerical values of a?, we may obtain two tables similar to 
that of Art. 11, from which the two corresponding cuiTes 
may be plotted. Fig. 2 shows these curves, XB^A^B^X being 
the one for the set of observations whose measure of precision 
is Aj or A, and XB^fi^X the one for the set whose measure 
of precision is A, or 2A. These curves show at a glance the 
relative probabilities of corresponding errors in the two sets; 
thus the probability of the error is twice as much in the 
second as in the first set, the probability of the error OP^ is 
nearly the same in each, while the probability of an error, 
twice as large as OP^ is much smaller in the second than in 
the first set. Now if we draw the lines P^B^t P^B^, so that 
the areas Pfi^Afi^P^ and P^^JS^^ are respectively one-half 
of the total areas of their corresponding curves, the line OP^ 
will be the probable error of an observation in the first set, 
and OP^ the probable error of one in the second set. If we 
represent these by the letters r^ and r^, we must in each case 
have the constant relation 

(28) Vi = 04769, V, = 0-4769; 

and since A, is twice A^, it follows that r^ must be one-half 
of Tj, as is represented in Fig. 2. 

Thus we see that the probable error is an error of such 
a magnitude that it is as likely as not that any assigned 
error will exceed it or fall short of it. Hence the probability 

that any error x taken at random will be less than r is ^ , 

and that it will be greater than r also ^. It is then an even 

wager that any error x will be greater or less than the pro- 
bable error r. 

Referring now to the case of the 600 bullet-marks on 
the target, let us imagine a circle described from the centre 
which shall include exactly 250 of those marks. The radius 
of this circle is the probable error of the marksman. If 
he is to aim and shoot once more, it will be an even chance 
that the bullet will strike within the circle, that is, that his 
error x will be less than his probable error r. If another 
rifleman, less skilful than the first, shoot also 500 bullets at 
the target^ and we draw another larger circle which includes 
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250 of them, the radius of this circle will be his probable 
error. If the circle of the first marksman be 4 inches and 
that of the second 12, we recognise from^our constant rela- 
tion between h and r, that the precision of the first is three 
times that of the second ; and while we can afford to wager 
one to one that a shot of the first man will fall within 4 
inches of the centre, we can only for the second man afford 
to make an even wager that it will fall within 12 inches of 
the centre. 

Methods will be hereafter given (Art. 24) by which A, 
and consequently r, may be determined for any given ob- 
servations. For the present we give the following numerical 
example, illustrating more fully the practical use of the pro- 
bable error. 



17. Let an angle be meastli'ed ten 'times equally care- 
fully by a theodoUte reading to 10", and again be measured 
the same number of times with a transit reading only to 1'. 
Suppose the results to be the following: 



Observation. 


By Theodolite. 


By Transit. 




f •> ff 


' ia . Q f •,* i.. 


1. 


24 13 40 


24 13 


2. 


13 10 


14 


3. 


13 30 


15- 


4. 


13 40 


13 


5. 


14 


12 


6. 


13 20 


13 


7. 


. 13 30 


13 


8. 


13 40 


12 


9. 


13 60 


14 


10. 


13 40 


15 


Sums... 


242 16 


242 14 


Averages... 


240 13' 86" 


240 13' 24" 



It will at once be perceived that the average of the 
theodolite observations is the more accurate and reliable, 
H(m much more precise and reliable will be shown by a 
comparison of their probable errors. From thp formulae of 
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Art. 24, the probable error of the mean of the theodoKte 
observations is 3"*1, and of that of the transit 13"'8, so that 
the results may be written 

by the theodolite, 24? 13' 36" ± 3"-l ; 

by the transit, 24^ X3' 24" ± 13"-8. 

The meaning of these results is, that our* confidence in 
the work of the theodolite is such that we could make an 
even wager that the trm value of the angle is 

between 24" 13' 36" + 3"! and 24" 18' 36" - 3"-l, 

but that as far as the transit work is concerned we could 
only afford to bet one against one that it Is 

between 24" 13' 24" + 13"-8 and 24" 13' 24" - 13"-8. 

The range of probable error in thd^-^ first is only about 
one-fourth of that in the second, and hence we recognise 
that the theodolite average is four times as precise as that 
of the transit. Also, if these ten observations were to be 
again repeated, with equal carefulness, it will be an even 
wager that the new averages will dififer from the true value 
of the angle by the respective quantities 3"*1 and 13"'8. 

In like manner we may find that the probable' error of 
a single observation is by the theodolite 9"'7, and by the 
transit 43"*5, and these observations signify that before the 
observations were taken, errors greater and less than those 
(quantities were equally probable. That this is sustained- 
in practice the reader may assure himself by assuming the 
averages in the above example as the true values of the 
angle and then computing the errors, and he will find in 
the one case 4 errors greater and 6 less than 9"'7, and in the 
other 4 greater and 6 less than 43"*5 ; which, considering that 
this is not an actual example but one written down at 
random, is a remarkably close agreement. 

The probable error then furnishes us with the means of 
comparing observations. The smaller that quantity the 
better and more reliable are the measurements. Since the 
probable error r is inversely proportional to the measure of 
precision h (Art. 16), we see that if the probable error of 
one set of measurements is one-half .that of a second set, 
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its precision is twice as great; a glance at the curves of 
Fig. 2 shows this graphicafly to the eye. The probable error 
gives moreover an absolute measure of our confidence in the 
accuracy of the work. 

Kinds of Ohservations. 

18. Direct ohservations are those which are made di- 
rectly upon the quantity whose magnitude is to be deter- 
mined. Such are, measurements of a line by direct chaining, 
of an angle by direct reading with a transit, etc They occur 
in the daily practice of every engineer. 

19. Indirect ohservations are not made upon the quantity 
whose size is to be measured, but upon some other quantity 
or quantities related to it. Such are, measurements of a 
line through a triangulation by means of a base and observed 
angles, of an angle by regarding it as the sum or difference 
of other angles, the determination of the diflference of level 
of two points by readings upon graduated rods set up at 
different places, the determination of latitude by observing 
the altitude of stars, etc. In fact the majority of observa- 
tions in engineering and physical science generally belong to 
this class. - 

20. Conditioned observations may be either direct or 
indirect, but are subject to some rigorous requirement or 
condition. As such may be mentioned : the three measured 
angles iii a plane triangle must be so adjusted that their 
sum shall be exactly 180^ the sum of all the percentages in 
a chemical analysis must equal 100, the sum of the nortnings 
must equal the sum of the southings in any traverse which 
begins and ends at the same point, etc. 

21. Measurements which are subject to no such rigorous 
conditions are called independent, meaning thereby that the 
observed quantities have no mutual dependence, so that all 
systems of values are in thought equally possible, and a 
variation of the value of one quantity need not necessarily 
affect the values of the others. In conditioned observations, 
however, all systems of values are not in thought equally 
possible, but only those can be admitted which exactly satisfy 
the rigorous conditions. 
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As an illustration of these classes let us consider the 
angles A OB and BOG having their vertices at the same 
point (Fig. 3). If we set up a transit at the point 0, 
and measure the angle AOB or 50(7, each of these measure- 
ments is a direct observation. If however in order to deter- 




mine these angles we damp the instrument so that, its zero 
point has some arbitaxy direction OM, and then measure the 
jtngles MO Ay MOB and MOG^ the observations are indirectn 
Moreover, whether observed directly ox indirectlv, the values 
obtained for AOB and BOG are independent of each other, 
iPut if we measure the angles AOB, BOCsnidAOG by either 
of the above methods, these observations are conditioned, or 
subject to the rigorous requirement that when finally ad- 
justed AOB plus BOG must equal AOG, and no system of 
values can be adopted for these three angles which does not 
exactly satisfy this condition. 

These kinds of observations suggest the division of our 
subject into three pai*ts, and we s^l accordingly treat in 
Chapter II. of the adjustment and comparison of direct ob- 
servations upon a single quantity; in Chapter III. of in- 
dependent observations, either direct or indirect, upon several 
quantities ; and in Chapter IV, of conditioned observations. 
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CHAPTER 11. 

DIBECT OBSERVATIONS UPON A SINGLE QUANTITT. 

22. The principle of the arithmetical mean has been 
from a very remote antiquity employed to obtain the ad- 
justed value of an unknown quantity, which is the object of 
direct measurement. Its universal acceptance as an axiomatic 
rule for adjusting such observations, and the simplicity of the 
process itself, would seem to require that it should be placed 
at the foundation of any system for the combination of nu- 
merical measures. "We must however warn the reader that 
the principle of the average is not applicable to any observa- 
tions except those of equally good direct measurements upon 
a single quantity, and that it cannot be used for the combina- 
tion of either direct or indirect observations upon several 
related quantities. A single illustration will show the danger 
of its indiscriminate use, and satisfy the reader that care is 
required to limit it to its proper sphere. 

Let A and B be two points whose elevations above a 
given datum are to be determined, and suppose, in order 
that the observations may be equally good, that these three 
points are situated at the vertices of an equilateral triangle. 
Now let a level be set up between and A and measure- 
ments taken which show A to be 10 feet higher than 0, then 
let it be moved to a point between and B, and readings 
taken which show -B to be 15 feet higher than 0; and lastly 
let it be set between A and B, and readings taken showing B 
to be 4 feet above A. 

These observations, it will be at once seen, are discordant, 
and we ask, how shall they be adjusted? Many persons 
would, we are afraid, proceed to apply the law of the average 
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thus : for the elevation of A we have by one route 10 feet, 
and by the other 11 feet» and hence A must be taken as the 
average 10^ feet ; also by one route JB is 15, and by the other 
14 feet, and hence £ should be called 14^ feet. This however 
is entirely wrong, as a moment's consideration will show, for 
by this adjustmeut the pieo^ureipent bet^yeen and A has 
been corrected from 10 to lOJ^, that between and B has 
been corrected from 15 to 14^, while that between A and B 
remains 14^ — lOJ = 4, or has not been corrected at all. The 
principle of the average does not therefore apply to such 
observations, or in fact to any except equally good observations 
upon cb single quantity. The proper adjusted heights of A 
and JB, we may mention, are to be determined by the methods 
of the next chapter and are lOJ ai>.d 14| feet respectively, 
results which suppose the correction of each observation by 
an equal amount. 

The Arithmetical Mean. 

23. The average or arithmetical me;9n has always been 
accepted and used as the best rule for combining direct 
observations of equal precision upon one and the same quan- 
tity. The following reasoning may be given to justify its use. 

Suppose a certain. angle is to be measured, and let its 
true but unknown value be z. A single reading taken by a 
transit gives the value M^. This being the only measurement, 
its value M^ must be accepted as the wiost probable result, 
and used instead of the true value ^gr. Now let a second 
reading be taken with the same instrument and under the 
same circumstances, which gives tho value M^. Then since 
there is no reason for preferring one result to the other, we 
consider that errors in excess and deficiency from the true 
value are equally probable (Art 10), aijid hence combine the 
results so that the diflferences z — M^ ajid z — M^ are nume- 
rically equ^, and this gives 






and this value of z we must 3<;cept as the most probable and 
must use it as representing, as far as our observations go, the 
true value. 
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26 WRECT 03SEEVATI0KS UPON A SINGLE QUANTITY^ 

In like manner if we have n observati(His giving .the 
values Jfj, Jf,, JkTj .•. Jf^, errors greater and less than the 
true value are equally jMrobable, and in a large number of 
observations -will be equally numerous^ and the algebraic 
sum of all such errors should be zero, or 

The solution of this equation will give us, as far as our n 
measurements are concerned, the value of z, but those mea- 
surements being limited in number and discordant, we cannot 
regard it as absolutely true, but only as probably true, and 
we call it the most probabk value. Designating it by ^o we 
have from the above equation 



(29) ^^^ M^ + M,+M, + ... + M, 



Hence in any series of equally precise observations upon a 
single quantity, the. arithmetical mean is the most probable 
value, that is, the most accurate value deducible from those 
observations. See Art. 66- 

Problem. The bearing of a line is taken five times with 
a solar compass giving the values, 

N.12'K, N;7'il, :N.2'W., K12'W., N.5'E.; 

what is the adjusted bearing? v 



Probable Error of a Single Observati(m and of the 
Arithmetical Mean. 

24. Let one and. the same quantity, for example an 
angle, be measured n times with equal care ; let z^ be the 
average or arithmetical mean of the measurements 

ifj, if,, M^, etc. 

Let each measurement be subtracted from the average z^, and 
let the diflferences or residuals be called v^, v,, v,, etc., so that 

(31) v,^z,-M^,v^^z^-M^,v^^z,^M^,etc 
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Also let 2t;* denote the sum of the squares of these 
residuals^ or 

2ii» = v^« ^ t;^« + t;,» 4- ... + v/. 

Then, as proved in Part IL, the prohahle error of a single 
observation is 

(42) r = 0-6745y^, 

and the probable error of the arithmetical mean z^, is 

(47) r. = ^-=0-6745^_^. 

To find then the probable errors r and r^, we have to find 
the average, subtract from it the several measurements thus 
forming the residuals, then to square each of these and add 
the products, and then to use the above formulae. 

The reader should particularly observe the distinction 
between the residual v and the error x (Arts. 11, 14) ; the 
former is the diflference between the most probable result and 
a single observation; the latter is the diflference between the 
true value and an observation, and of course can never be 
exactly known since we can never be certain of having found 
the true value. Thus if z be the true value of an angle, and 
«g the average of the observations ifj, if,, M^, etc., the errors 
are (z — M^), (z — M^, etc., while the residuals are (Zf^ — MJ, 
(z^-M^, etc. 

• 25. The probable error gives us the means of comparing 
the precision of our observations, and of determining the 
degree of confidence which we can place in them. An illus- 
tration of its use has already been given in Art. 17. We 
append another example worked out in full. 

In the U, S. Coast Survey/ Iteport for 1854, the following 
24 measurements are given of the angle Quaker-Pocasset- 
Beacon, taken at the station Pocasset in Mass., each observa- 
tion being the result of an equal number of repetitions and of 
the same precision or weight. 

We first determine the average or the most probable 
value of the angle by adding the readings and dividing the 
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28 DIRECT OBSERVATIONS UPON A SINGLE QUANTITY. 

sum by 21. This is. 116' 43' 4y'-64. Then subtracting the 
first reading from this we have the residual t;^ = 5*19 which 



No. 


Reading. 


«. 


"»*• 




, ^ 






1. 


116 . 48 . 44-45 


5-19 


26-94 


2. 


,60-65 


-0-91 


•83. 


3. 


60-95 


-1-31 


1-72 


4. 


48-90 


. 0-74 


-65 


6. 


49-20 . 


0-44 


•19 


6. 


48-85 


0-7« 


-63 


7. 


47-40 


2-24 


5-02 


8. 


47-7^ ■ 


1-89 


8-57 


9. 


61*05 


-1-41 


2-00 


10. 


47-85 


1-79 


3-20 


11. 


SO -60 


-0-96 


•92 


12. 


48-45 


1-19 


1-42 


13. 


61-75 


-2-11 


4-45 


1*. ' 


4900 


0-64 


•41 


15. 


52-35 


-2^71 


7-34 


16. 


61-30 


-1-66 


2-75 


17. 


61-05 


-1-41 


2-00 


18. 


61-70 


-2K)6 


4-24 


19. 


49'-05 


0-69 


•35 


20. 


60-65 


-0-91 


•83 


21. 


49-26 


0-39 


•15 


22. 


46-75 


2-89 


8-35 


23. 


49-25 


0-39 


•15 


24. 


63-40 


-3-76 


14-14 


Meanzo 


= 116« . 43' . 49"-64 




Si;^=92 15 



we place in the polumn headed v: and in like manner taking 
the difiference between the mean apd each reading we fill out 
that column with the values v^, v,, Vj, etc Then by a table 
of squares we form the numbers »j" = 26 94, v^* = 0*83, etc., 
which we place in the column ^;^ Adding these we have the 
sum 2v'= 92*15. Then as w = 24we have, from the pre- 
ceding formuljB, the probable error of a single observation, or 



r= 0-674 
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BIRECT OBSERVATIONS TTPON A SINGLE QXTANTITY. 29 
and the probable error of the mean z^ is 

Hence the adjusted value of the angle may be written 

116'*3'49"-64!±0"-275. 

If theil these 24 observations were to be repeated under 
the same circumstances, it would be an even wager* that 
there would be 12 errors greater and 12 less that 1"'349, and 
also an even wager that the mean would diflfer from the true 
value of the angle by 0"'275. Our confidence then in the 
above mean 116*43 49"*64 is such that we regard 0"-275 
as the error to which it is liable, that is, it is an even wager 
that the mean is within 0"*275 of the true value, and of 
course also an even wager that it exceeds the true value by 
that amount. 

26. From the above values given for the probable errors 
we observe that 

If we denote the measures of precision corresponding to 
r and r^ by h and h^, we have (Art. 16) 

,„Q, 0-4769 , 0-4769 

(28) r = -^^andr, = -^; 

and inserting these in the expression above, we find 

(46) h,^h\/n, 

that is, the precision of the arithmetical mean increases as the 
square root of the number of observations. 

In order then to make the value of the angle in the 
above examples twice as precise, that is, make the probable 
error of the mean one-half as large, we must have four times 
as many observations, or 96. Let the reader test this rule 
by taking at random any six of those observations, and find 
the probable errors. The probable error of the mean will 
be approximately twice that given above, or 
2xO"-275 = 0"-55, 
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while the probable error of a single observatio& will remain 
nearly the same as before, or l"'So. 

Problems. !• A base line is measured five times with 
a steel tape reading to hundredths of a foot, and also five 
times with a chain reading to tenths oi^ a foot, with the 
following results: 

By the tape. By the chain. 

74117 feet 741-2 feet 

74109 ... 741-4 ... 
741-22 ,.. 7410 ... 
74112 •„ 741-3 ... 

74110 ,^ 7411 ... 
What are the averages And their probable errors ? 

Ana. By tape, 741 14 i OOIC. 
2. A line is measured five tjmes, ajid the probable error 
of the mean is 0-016 feet. How many additional measure- 
ments of the same precision are necessary in order that the 
probable error of the mean shall be only 0*004 feet? 

Weights of Observations. The General Mean. 

27. We have thus far considered our observations as 
equally precise, that is, as made with tlie same instrument^ 
and under exactly the same circumstances. We now come 
to the case of observations of unequal precision, or as 
we usually say of unequal vjeight. The sense in which we 
use this word we will try to make clear by a practical 
illustration. Suppose a line to be measured 20 times 
with the same chain. It) measurements giving the value 
9342 feet, 7 giving the value 934'0 feet, and 3 giving 
934-4. The adjustment of these results is effected by the 
above method of the average, by writing the first value 
10 times, the second 7 and the third 3, adding the numbers 
and dividing their sum by 20 ; this process is evidently the 
same as multiplying each value by the number of times it 
occurs, and dividing the siim of the three products by 20, or 
10 X 934-2 ^7x 934-0 -f 3 x 934-4 
20 

= 934o + ^^^^^:||-^^* = 9341C. 
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Now if the 10 measurements, instead of giving each time 
934-2 feet, had given 10 results whose mean was that 
number, 10 would be the weight of 934*2, or considering 
934*2 as a single observation, 10- is said to be its weight 
Also 7, 3, and 20 are the weights of 9340, 934*4, and 93416: 
the first being the equivalent of 7, the second of 3^ and the 
third of 20 observations, eadi equally good« 

By the term weights, then, we mean numbers related to 
the accuracy of the measurements, so that an observation of 
the weight 8 is to be regardied as equal to 8 observations of 
the weight 1 ; and observations havmg the same weight are 
to be considered as equally good. The average of n equally 
good observations, being the equivalent of those observations, 
has thug a weight of fi. 

. The combination of weighted observations is to be made 
by the process illustrated by the above numerical example ; 
thus if g., g^, g^ etc. be the weights of the measurements 
M^, M^y Jfj, etc., the adjusted value is 

or, the most probahle value of the measured quantity is found 
hy multiplying each observation by its weight, and dividing 
tike sum of the products by the sum of the weights. This value 
^is called the General mean to distinguish it from the arith- 
m^ical mean z^, which is only to be used when all the 
weights are equal. 

Weights should be carefully distinguished from measures 
of precision (Art. 16): the former are relative numbers, which 
are usually so taken as to be free from fractions; the latter 
are absolute quantities. The relation between them will be 
shown in the next article. 

JProblems. L An angle is measured 20 times with the 
same theodolite. The mean of 6 readings is 27^ 34' 32", the 
mean of 10 is 27* 34' 40", and of the other 4 is 27° 34' 48". 
What is the ^justed value ? Ans. 2T 34' 39"*2. 

2. The bearing of a certain line is taken as N, 89° 45' W., 
and as S. 89* 45' W. If the weight of the first observation 

jogle 
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is 13, and that of the second 2, what is the most probahle 
value of the bearing ? Ans. N. 89® 49' W, whose weight is 15, 

28. To establish the relation between weights and pro- 
bable errors, let us consider n observations upon the same 
quantity, giving the results Jfp J^, M^, etc. whose measures 
of precision are A-,, h^, A,, etc. If z is the true value of the 
quantity, the errors are (^9 — ifj —fl?^, (-2: — if J = ar^, etc. 
Then from Art. 11, 

Prob. of the error a?^ =yj =* d^ ^-^i^^^^ 

' - «^« = y2 = Ca« % 



«^t=ys=c.«"*^. 



etc., etc. 

From Art. 8 the probability P of committing all these errors 
is the product of the probabilities y^, y,, y^, etc., or 

(25) P = c,c,C3...e-^*^'^'+*«'^«'+*^*-'^«*^>; 

and the most probable value of z is that for which P is 
a maximum (Art. 9), and in order that P should be a 
•maximum, the quantity, 

(26) A,V + A,V+V^3« + etc., 

must be a minimum. Inserting in this the values of x^, 
iCj, a?g, etc. in terms of z, we have 

A," {z - M^y + A, (« - M^y + A3 (« - MJ" -f etc. = a minimum. 

Differentiating it with reference to z, dividing by 2dz and 
placing the derivative equal to zero, we have 

V(^-J/J + V(«-^s)+V(^-^3)+etc. = 0. 

Solving this equation, and denoting the resulting value of 
z by Z, we find the most probable result to be 

. . Vif, + Vitf, + h,\¥, + etc. 

The value of Z given by this expression must be the same 
as that given by the general mean in Art. 27. Comparing 
the two expi:essions, we see that 
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or, the weights of observations are proportional to the squares 
of their measures of precision. 

A direct application of these principles can be made to 
the veiy common case of diflferent sets of observations of 
unequal precision, arising from measurements by different 
instruments of one and the same quantity. Let the number 
of observations in the first series be n,, in the second n^, in 
the third w., etc.; let the mean as given by the first series 
be 2?,, by the second z^, by the third «,, etc.; and let the 
weights of these be ff^, ff^, ff^, etc.; and the corresponding 
measures of precision be h^/h^^\^ etc. Then from the above 
principle 

(49) 9^'9»- 9, ■■■■ K •- K ■ K- 

Inserting in these the values of A^, A,, etc., which are 
determined in Part II., we have (after striking out the 
common factor 2) 

in which St;" denotes the sum of the squares of the residuals 
in the first set, 2v"* in the second, etc. (Art. 24). We have 
then in such»a case to find the average of each set, next the 
weights of these averages by the proportion just given, and 
then the most probable value of the measured quantity by 
the general mean (Art. 27) of the several arithmetical 
means, that is, by the formula 

(50) 5- ^A+.y2^2 + .y8g8 + etc, 

91+92 + 9,-^^^' 

To illustrate this very common case we give the follow- 
ing example and problems. 

A certain line was measured by three different surveying 
parties, using three different chains ; the first party measured 
it 5 times with the results given in column I. below, the 
second 6 times with results as in column II., and the third 
4 times as shown in column III. What are the relative 
weights of the three means and the most probable length of 
the line ? 
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We first determine the three means z^^ z^, and z^: then 
subtracting each observation from its mean, find the residuals 



I. 


if. 


v'J. 


II. 


«•. 


«"«. 


III. 


tr"'. 


V"*, 


6110 


2 


4 


4980 


120 


14400 


6105 








.5090 


22 


484 


6100 








6100 


6 


25 


6140 


-28 


794 


6220 


-120 


14400 


6110 


-6 


25 


6100 


12 


144 


6160 


-60 


S600 


6106 








6120 


-8 


64 


6040 
6100 


60 



8600 









Zi = 6112 


2i;'*=U80 


2a =6100 


2t;"a= 36000 


Z8=6105 


Si;'"»=60 



which are placed in the columns v\ v", and v". From a table 
of squares we take the squares of these numbers and place 
them in the columns headed v\ v'", and v"^: and by their 
addition find the sums of those squares. For the first series 
of measurements Wj = 6, for the second n^ = 6, and for the 
third ^3 = 4. Designating the weights of the means by 
ffv ffiy ^^^ ffs* w® ^^v® fr^^ ^^® above proportion 



9i ' 9t "' 9^ •• lARn 



5 X 4 
1480 



6x5 4x3 
36000 • 50 ' 



or by reduction 

9x ' 92' 9z •• 600 • 37.: 10656; 

hence the mean z^ is equivalent to 600 observations of the 
weight 1, z^ to 37, and z^ to 10656 ; the second is hence the 
least, and the third the most reliable. The most probable 
value of the line then is 

„ 600 X 5112 + 37 X 5100 + 10656 x 5105 ^- ^. „^ 
^= 600 + 37 + 10656 =5105-36- 

Problems. 1. An angle is measured four times with a 
theodolite, six times with a transit, and five times with a 
sextant, giving the observations : 
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By the theodolite. 


By the transit. 


By the sextant. 


6" 17' 5" 


6" 17' 


6" 17' 20" 


6 17 10 


6 16 . 


6 17 


6 17 


6 15 


6 17 40 


6 17 5 


6 19 


6 16 50 


-,«v ^^'^ 


6 17 
6 18 


6 17 10 



What are the relative weights of the means, and the 
most probable value of the angle ? 

Ans. Z=6nr6"-36. 

2. Two sets of measurements are made upon the same 
angle by different observers with the same instrument. Each 
takes 5 observations, as follows : 



First Obserrer. 


Second Observer. 


47° 23' 40" 


47" 


' 23' 30" 


23 45 




23 40 


23 30 




23 50 


23 35 




24 


23 40 




23 20 



What are the relative weights of the observers, and the 
adjusted value of the angle ? 

Ans. g^ : g^ :: 100 : 13. 

29. Observations may also be combined when their 
probable errors are known. From Art. 28 we have the pro- 
portion 

(49) g, : g, : g, :: V - K • h\ 
and also from Art. 16 the relation 

(28) K:K-K---p--^^-^.-> 

'l '% 'i 

hence by comparison of the two proportions 

/XON 111 

(53) g, :g,:gs :: ri • -« • p; 

'i 'a 'a 

that is, the weights of observations are inversely proportional 
$0 the squares of their probable errors, HeUce having com- 

3-41 
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puted the probable errors of two arithmetical means, the 
relative weights can be immediately found, and the two 
means then combined by the general mean. 

Thus in Art. 17 we have a case where the same angle is 
measured by a theodolite and by a transit, giving the means 
and probable erro^rs, 

by the theodolite, 24M3' 36" ± 3"1; 

by the transit, 24' 13' 24" ± 13"-8: 

Designating the weights of these by g^ and g^, we have 
from the above principle, 

9i '9%'- gTp • 13^ •• ^^^ • ^ ^^^^^7' 
and the most probable value of the angle then is 

Z= 24« 13' + ^^^''^^125^''^^^ = ^*' ^^' ^^''**- 

Problems. 1. "What is the most probable value of the 
base line in Prob. 1 of Art. 26, as given by the two sets of 
observations? Ans. 741*146. 

2. Three sets of observations upon the same quantity 
give the following averages and probable errors, 

803-4 + 0-4; 803-2 ±03; 8031 + 0-08. 

What is the adjusted value of the quantity ? 

Probable Error of the General Mean, 

30. Having determined the adjusted value of obsei-va- 
tions of unequal weights by the preceding methods, we next 
inquire what is our degree of confidence in that result, or 
what is its probable error (Art. 16). 

Let n be the number of observations, or sets of measure- 
ments, j7i, g^, g^, etc. their relative weights; let -^ be the 
general mean as found by Art. 27, and O its weight. Also 
tet Vj, Vj, Vg, etc. be the residuals or differences between 
the general mean and each observation, and Xgv^ the sum^ 
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^9x^x -^ 9Jo^ -^ g^v^ -{- (dio,, that is, the sum of the quaptities 
formed by multiplying the square of each residual by its cor- 
responding weight. Then, as proved in Part II., the probable 
error of an observation whose weight is unity is 



(65) r = 0-6745 y^l^, 



Having found r, the probable error of the general mean, 
or of any observation, may be found by the principle of 
Art. 29. If (? be the weight of the general mean and R its 
probable error, g^ the weight of a given observation and r^ 
its probable error, we have 

(53) G : 5^x : 1 :: ;gi : ^ ' ^> 
from which we find 

(56) i = -^ and r. = -^; 

that is, the prohahle error of an observation whose weight is 
known, is equal to the probable error of an observation of the 
weight 1, divided by the square root of the given weight The 
weight G of the general mean is always equal to 

5^1 + 5^2+5^8 + etc. (Art. 27). 

If we have then several sets of measurements as in 
Art. 28, whose averages are z^, z^, z^, etc., we have to find 
their relative weights^,, g^, g^y etc., and their general mean 
Z. Subtracting each average from the general mean gives us 
the residuals v,, v,, etc.; multiplying the square of each by its 
weight and adding the products gives us the sum ^gi^. Then 
n being the number of sets, the above formulae famish us with 
the probable errors. Further, if we have found the probable 
errors r^, r^, etc. of the averages 2?^, «,, etc., the principle of 
Art. 29 gives us for the probable error of the general mean 



(54) 



-B = '-VS = *"«\/|=***-' 



which may in some cases be more convenient to use than the 
formulaB above. The following examples and problems will 
illustrate the application of the formulae. 
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1. Suppose that the observations in the example of 
Art. 25 are given as in the following column z^ the mean of 
the first five being 48"*81 with the weight 5, the mean of the 
next four 48"-76 with the weight 4, and so on. Then the 
operation for finding the probable error of the general mean 
is thus exhibited: 



g- 


0. 


«. 


««. 


flrt^. 


5 
4 
5 
3 
2 
5 


/ // 
116 . 43 . 48-81 
48-76 
49-53 
51-56 
50-38 
49-84 


0-83 

0-88 

0-11 

-1-92 

-0-74 

-0-20 


0-69 

0-77 . 

0-01 

3-69 

0-65 

0-04 


3-45 
3 08 
0-06 
11-07 
1-10 
0-20 


G=24 


^=1160 43' 49-64 


2iru«=18-95 



The general mean ^of course here agrees with the average 
found in Art. 24, and its weight O is the sum of the several 
weights or 24 the number of single observations. Subtracting 
each average firom the general mean we have the residuals in 
the column v, and from a table of squares we place their 
squares in the column v* ; multiplying each of these by its 
corresponding weight we have the quantities in the column 
g'^ whose sum is 18*95. Then, n being 6, we have from the 
formula 



^0-6745 



^/ 



18-95 



= l"-32. 



This is the probable error of an observation of the weight 
unity, and should hence agree with that of a single observa- 
tion 1"'35, as found in Art. 24. The discrepancy is due to 
the small value of n (see Prob. 2), The probable error of the 
general mean is then 
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which agrees suflSciently well for most practical comparisons 
with that found before. 

2. In Prob. 1 of Art. 28 the probable error of the theodo- 
lite mean is 1"*4, of the transit 23"*3, and of the sextant 5"*8, 
and their relative weights are nearly 288, 1 and 16. Hence 
the probable error of the general mean is by our second 
formula 

The first formula gives a less result, as w is only 3. Strict 
agreement in such results cannot be expected, since the 
theory upon which the formulae are deduced supposes n to be 
a large number, and in our examples illustrating their appli- 
cation we are obliged to choose cases involving but few 
observations. 

JProblems, 3. What is the general mean of the observa- 
tions on the base line in the question of Art. 26, and its 
probable error ? Ana. 74iri46 + 0012. 

4. Eight observations of a quantity give the results 769, 
768, 767, 766, 765, 764, 763 and 762, whose relative weights 
are 1, 2, 3, 4, 5, 6, 7 and 8. What is the probable error of 
the general mean, and the probable error of each obser- 
vation ? 

Becapitulation. 

31. We have now given and illustrated the methods for 
adjusting and comparing direct observations upon a single 
quantity. They fall under three heads : 

Ist. K all the observations are equally good or of equal 
weight, the average is the most probable result (Art. 23). 
The degree of confidence which we can. place in one of those 
observations, or in the average, is shown by their probable 
errors (Art. 24). 

2nd. If there are several sets of observations made under 
difierent circumstances or by diflferent instruments, the most 
probable result is given by the general mean (Art. 27), and 
in order to obtain that result the relative weights of the 
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several averages must be deduced by the proportion of 
Art. 28. The precision of the several means and of the final 
general mean is also shown by their probable errors (Art. 30). 

Srd. If we have single observations which are known to 
be of unequal precision, we have no means of finding their 
weights as in the preceding case, but must assign to them 
such weights as they seem to deserve in our judgment. 
Thus if an augle be measured once by a theodolite reading 
to 20", and once by a transit reading to l'> we recognise that 
the first is the more reliable, and should, in finding the 
mean, give it a weight of about 9 times as much as the second. 
By taking series of observations with two such instruments, 
their relative weights may be found and recorded for use in 
cases where single measurements arise. When such know- 
ledge does not exist, weights may be assigned, and then the 
general mean found as in Art. 27. The assigning of weights 
in such case is of course a matter requiring experience and 
judgment. 

The combination of direct observations upon different hnt 
related quantities is considered in the two following chapters. 
The most simple cases of such adjustments, such as finding 
the length of a base line, which has been measured in several 
portions, are referred to in Arts. 41 and 42, in connection with 
the discussion of their probable errors. 
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CHAPTER III. 

INDEPENDENT OBSERVATIONS UPON SEVERAL QUANTITIES. 

32. In Arts. 18 and 19 we divided observations into 
direct and indirect, the former being made upon the quanti- 
ties to be determined, and the latter upon other quantities 
related to them. The line of demarcation between the two 
is however not very distinct, nor is it necessary that we 
should be able to point it out clearly. For practical purposes 
independent observations need only to be distinguished as 
those upon one quantity, which have been already considered, 
and as those involving Toore than one quantity, which we are 
now to take up. The methods and formulae of the preceding 
chapter which treat of one quantity, although for convenience 
given separately, are indeea but particular cases of those now 
to be developed. 

Independent observations, whether direct or indirect, 
which ai*e made to determine the magnitudes of quantities, 
are generally represented by equations which we shall call 
observation equations. To illustrate how they arise let us 
consider the following practical case. Let Fig. 4. 

O (Fig. 4) represent a given bench-mark, 
and 8^ T, U three points whose elevations 
above are to be determined. Let five 
lines of levels be run between these points 
as indicated by the dotted lines of the 
figure, giving the following results : 

Observation 1. 8 above = 10 feet, 

% T fif= 7 ... 

3. T 0=18 ... 

4. T U^ 9 .., 

5. U'belowiSf= 2 ..• 
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It will be at once perceived that the measurements are 
discordant ; if we take observations 1 and 2 as correct, the 
height of S is 10 feet and T is 17 ; if 2 and 3 are correct, 8 
is 11 and Tis 18 feet, etc. : and in general it wiU be found 
impossible to find a system of values which will exactly satisfy 
all the observations. If we designate the elevations of the 
points S, T and Z7by the letters 5, t and w, the observations 
furnish the following equations : 

« = 10, 
<-5= 7, 
<=18, 
«-tt= 9, 
fi — M= 2, 
each one of which is an approximation to the truth, but all 
of which cannot be correct. The number of these equations 
is 5, the number of the unknown quantities is 3, and hence 
an exact solution cannot be made by algebraic processes. 
It being impossible then to find values of 5, f, and u which 
will satisfy all the equations, we must be content with deter- 
mining their moBt probable values (Art. 9). 

So, in general, independent observations upon several 
quantities give rise to independent observation equations 
greater in number than the unknown quantities to be deter- 
mined ; and our problem is^ to find, out of the many systems 
of values, all equally possible, which may be assigned to the 
unknown quantities, a system which is the most probable, and 
hence the best. 

Solution of Observation Equations. 

33. "We take up first the case of measurements of equal 
precision or of equal weight. Let IT, if^, Jlfg, etc. be the 
numerical results of the observations which are n in number, 
and which are made upon quantities related to the quantities 
S{ t, u, etc., w^hose values are to be found. Let the observa- 
tions give rise to the following observation equations : 
a^s + hf + c^u + etc. = if^, 
(66) a^s + bj, + Cjti + etc. = M^, 
a^s + Ij^ + CgM + etc. = i/g, 
6uC», ei/C., 
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in whicli a^, a,, b^, \, etc. are the known coefficients of the 
unknown quantities. The number of the measurements is n, 
and if the number of unknown quantities were also n the 
solution of the equations could at once be made. But as we 
have seen the number of the latter is usually less than n, and 
an exact solution is impossible. We must therefore be content 
with finding the most probable values of s, t, u, etc. 

Whatever system of values is chosen for 8, t, u, etc. it will 
not exactly satisfy each of the above equations, since the 
measurements M^, M^, M^, etc. are imperfect. If then we 
consider s, t, u, etc. as representing the true values of the 
quantities, the above equations may be written 

a^s + hf + cju, + etc. — ilf^ = x^, 
(67) a^8 + 6,< + c,w + etc. - M^^x^y 

etc., etc., 

that is to- say, they do not reduce exactly to zero for any values 
of the unknown quantities, but leave small difierences or 
errors x^, x^, x^, etc. Now by our general principle of Art. 14 
the most probable values of the unknown quantities are those 
which make the sum of the squares of the errors a minimum, 
and hence the most probable values of 8, <, w, etc. are those 
which make x^ + x^ + x^ + etc. the least possible. We pro- 
cede to develope a method for finding those quantities^ 

Let us first consider what is the most probable value of 
the unknown quantity 8. As we need only to regard 5, let 
us denote the terms in the above equations, independent of 8, 
by the letters N^,N^,If^, etc. Then they become 

a,8 + N^ = x^, 

a^8 + N^ = x^^ 

a^s-hN^^x^, 

etc. 

Squaring both terms of each of these equations and 
adding the results, we have 

(a,8 + N;)^ + {a^8 + N;)^+{a^8 + N;f + etc.^x^^ + x^*+x^'+eic. 

According to the principle above stated this quantity is 
to be made a minimum to give the most probable value of 8. 
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-Differentiating it with respect to 5, placing the first 3iflFeren- 
tial coefficient equal to zero, and dividing by 2, we have 

a^ (a,» -f JV^,) + a, (a,« + iV,) + ttg (a3« + JV;) + etc. = 0, 

and this is the equation which furnishes us with the most 
probable value of a. Hence we have the principle : To form 
the equation which gives the most probable value for one of the 
unknown quantities a^ s^we multiply each of the observation 
equations by the coefficient of s in that equation and add the 
results. In the same way, to find the equation for t we 
multiply each observation equation by the coeflScient of < in 
that equation and add the results. The equations thus 
formed are called normal equations; they will be the same in 
number as the number of the unknown quantities, and will 
be satisfied by only one system of values of the unknown 
quantities, which will therefore be the most probable system. 

34, For illustration, let us suppose that four measure- 
ments upon three unknown quantities have given the obser- 
vation equations 

5-^+2i*= 3 (1), 

3s + 2«-5tt= 5 (2), 

45 + « + 4w = 21 '. (3), 

-5 + 3^+3^ = 14 (4), 

from which it is required to find the best system of values 
of Sy t, and u. To form the normal equation for s, we must 
multiply each equation by the coefficient of s in that equa- 
tion and add the results; hence multiplying equation (1) by 1, 
equation (2) by 3, equation (3) by 4, and equation (4) by 
— 1, we have 

5-< + 2w= 3, 
95 + 6^-15^=15, 
165+4e + 16w = 84, 
«-3«-3w=-14; 
and adding these, we have the first normal equation 

275+6^ = 88 (5). 

We must now perform the same operation for t : multi- 
plying equation (1) by — 1, equation (2) by 2, equation (3) 
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by 1, and equation (4) by 3, and adding the results, we have 
the normal equation for t, viz. 

68 + 15t+u=-70 (6). 

In like manner, we multiply equation (1) by 2, equation 
(2) by — 5, equation (3) by 4, equation (4) by 3, and add the 
results to form the normal equation for u, viz. 

< + 64u = 107 (7). 

By this process we have three normal equations (5), (6) 
and (7), containing only three unknown quantities, and 
solving these by any of the algebraic methods, we find 

« = 2-4702, * = 3-5509, w=: 19157. 

If we substitute these values in equations (1), (2), (3) 
and (4) we shall find that they will not reduce to zero, but 
give the residuals 

t;, = - 0-2493, v, = - 0*0661, ^3 = 00945, t;,=.- 0*0704, 

the sum of whose squares is 00804. (Note, — the student 
should not forget the distinction between the errors or^, 
x^f etc. and the resid^ls v^, v^, etc. See Art. 24.) This 
quantity 0-0804 is less than the sum of the squares of the 
residuals resulting from any other values of s, t, and u. Let 
the reader test this by trying other values, for instance 

5=2f, <=3f, and tt=lf 

To solve observation equations of equal weight, we have 
then the following : For each of the unknown quantities f<yrm 
a normal equation by multiplying each observation equation 
by the coefficient of that unknown quantity in that equation 
(taJcen with its proper sign), otmd adding the results. Then there 
will be as mxiny normal equations as unknown quantities, and 
their solution will give the most probable values of those 
unknown quantities. In forming the normal equations it 
should be particularly noticed that the signs of the coefii- 
cients (+ or — ) are to be observed in performing the multi- 
plications ; and also that when the unknown quantity under 
consideration does not occur in an observation equation its 
coefficient is 0, For further illustration we give an addi- 
tional example, and a couple of problems as exercises for 
• the student. 
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If "we have given the observation equations 
8 = 14, 

t = 20, 

the coefficient of s in the first is 1, in the second — 1, and 
in the third 0. Hence multiplying each equation by these 
numbers and adding the results, we have 

25-^ = 7 

as the normal equation for s. Also the coefficient of t in the 
first equation is 0, in the second 1, and in the third 1, and 
in like manner the normal equation for t is 

2^-5=27. 
Solving the two normal equations, we find as the best values 

s = 13f, and « = 20J. 

Problems. !• Form the nonnal equations, and find the 
most probable values of s, t and u from the observation 
equations stated in Art 32. 

Ans, The normal equations are ^ 

3s — * — w = 5, 

-5-«4-2w = -ll, 
from which s = lOf , t = 17f , and w = 8 J feet. 

2. Find the most probable values of x, y, and z from 
the observation equations 

w= 5, * 

y + a?=ll, 
^ + y + a? = 13, 
«-f y — a; = 4. 

Ans. (c = 4§, etc. 

Adjustment of Independent Observations of equal weight 

35. The preceding principles furnish us with the follow- 
ing method of adjusting independent observations, either 
direct or indirect, upon several quantities. 
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1st. Represent each quantity to be determined by a 
symbol $, t, u, etc., and for each observation write an observa- 
tion eqvxition (Art. 32). 

2nd. From the observation equations form the normal 
eqiiations (Arts. 33, 34), which will be as many as there are 
unknown quantities. 

3rd. Solve the normal equations by any convenient 
algebraic method; the resultiog values of the unknown 
quantities will be their most probable values (Art. 33), that 
is, the best values which can be deduced from the given 
observations. 

The following example^ will illustrate the application 
of the method to cases arising in ordinary engineering 
practice. 

I. Adjustment of level lines. In the Report of the U, 8. 
Geological and Geographical Survey of the Territories for 
1873, Mr J. T. Gardner gives the following measurements 
as deduced from an examination of railroad profiles and 
coast survey levels. 

1. 8 above 0, 573*08 feet, by Coast Survey, and Canal 

levels, via Albany. 



2. 


T 


9> 


8, 2-60 


99 


)> 


Observations on surface 
of Lake Erie. 


8. 


T 


tf 


0, 575-27 


» 


99 


Coast Survey and R K. 
levels, via Albany. 


4. 


U 


99 


T, 16733 


9> 


» 


K. E. levels. 


6. 


X 


99 


U, 3-80 


99 


99 


99 99 


6. 


X 


99 


Z 170-28 


9> 


99 


„ via Alliance and 
Crestline. 


7. 


X 


99 


Y, 425-00 


» 


» 


» )f 


8. 


Y 


99 


0, 319-91 


99 


99 


E. E. and Coast Survey 
levels, via Phil. 


9. 


Y 


99 


0, 319-75 


99 


99 


E.E. levels, via Baltimore, 


In wh 


ich 













is the mean surface of the Atlantic Ocean,. 
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8 is the meaa surface of Lake Erie at Buffalo, 

T is Cleveland city datum plane, 

U is Depot track at Columbus, Ohio, 

X is Union Depot track at Pittsburg, 

Yis Depot track at Harrisburg. 

It is required to find the most probable elevations of 
each of these points above the datum 0, as given by these 
nine observations, supposed here to be of equal reliability or 
weight. 

We represent the unknown heights of 8, T, U, X and Y 
by the letters 5, t, u, x and y. Then the observations give 
us the equations 

« = 57308, 

< - 5 = 2-60, 

< = 575-27, 

n- * = 167-33, 

x-^u^ 3-80, 

a:- * = 170-28, 

a? - y = 42500, 

y = 319 91, 

y = 319-75. 

' All the coefficients of the unknown quantities are either 
+ 1 or — 1. Multiplying each equation in which b occurs by 
its coefficient in that equation, and adding the products, we 
form the normal equation for s ; then multiplyiDg each equa- 
tion in which t occurs by its coefficient, etc., gives us the 
second normal equation. Thus we have the fivo normal 
equations 

25- t =570-48=^, 

-5 + 4^- w- X =240-26 = 5, 

- < + 2u- X =163-53 = C, 

- i- w + 3a;- y=59908 = D, 

- a? + 3y = 214-66 = jB; 

containing only five unknown quantities. By w;hatever pro- 
cess these equations be solved, the values found for 5, f, u^ etc. 
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will be the same. We leave the reader then to choose his 
own process, simply remarking that the method of indeter- 
minate multipliers will prove the shortest. Representing the 
numerical terms by the letters A, B, 0, D and H, we find 

515=32^ + 135 + 11(7+ 92) + 3^= 29213-27, 

51t = 13^ + 265 + 22(7+ 18jD + 6iEr= 2933214, 

51m = 11^ + 225 + 60(7 + 27D + 9^« 3784470, 

I7aj= 3^+ 65+ 9(7+ 12D + 4JS?= 12672-37, 

17y= ^+ 25+ 3(7+ 4JD + 7^= 5440*53; 

hence we have 

5= 572-81, while Mr Gardnee gives 573-08, 

« = 575-14, 575-68, 

w = 742-05, 742-60, 

a = 745-43, 74600, 

y = 320-05, 319-91/ 

The discrepancy in the results is mainly due to the fact 
that Mr Gaedner has coDsidered the observations as of very 
unequal weight, taking s and t for instance, as given hy 1 
and 2 alone without reference to the other measurements 
(see Art^. 37 and 40). 

If these values be substituted in the observation equations, 
the residuals and their squares will be 

From our values. From Gabdneb's values. 



No. 


V. 


V*. 


1. 


0-27 


0-073 


s2. 


•27 


•073 


3. 


•18 


•017 


4. 


•42 


•176 


5. 


•42 


•176 


6. 


, -01 


•000 


7. 


•38 


•144 


8. 


•14 


•020 


9. 


•30 


•090 




v 


i;«=0-769 



Ifo. 


9. 


«a. 


1. 


0-00 


0-000 


2. 


•00 


•000 


S. 


f41 


•168 


4. 


•40 


•160 


5. 


•40 


•160 


6. 


•04 


•002 


7. 


109 


1-188 


8. 


•00 


•000 


9. 


!l6 


•026 




2 


t;*=l-704 
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The sum 0*769 of the squares of the residuals is the least 
arising from all the systems of values which can be attributed 
tp 8, t, u, etc. Mr Gardner's values, although exactly satisfy-* 
ing three of the equations, give the corresponding sum 1*704^ 

The above process, though simple, can be considerably 
abridged in the numerical operations by assuming approxi- 
mate values for the heights of S, T, etc., and regarding the 
unknown quantities as corrections to be applied to those 
assumed vsdues. We see at once from the observations that 
673 anji 675 feet would be approximate heights for 5 and T. 
If theii we place 

^=673 + 5', 

« = 575 + <', 
.. t^ = 742 + tt', . 

a? = 745 + a, 

y = 320+y, 

in which «', /, t*', etc. are corrections to be applied to the 
approximate elevations 673, 675, 742, eta, we may insert 
these expressions for 8, t, w, etc. in the observation equations 
above and obtain 

8^ 008, 

^-fi'= 0-60, 

t'= 0-27; 

t^'~f'= 0-33, 

af-^u'=^ 0-80, 

aj'-^= 0-28, 

;c'-y'= 0-00, 

y' = -0-09, 

y = -0-25. 

From these we form the normal equations by the same 
process as before, and have 

2«'~ If =-0-62=^', 

^ 8+U- w'- x' = 0-26 = -B', 

- «'+2w'- x' =-.0-47=0', 

- «'- tt+3x'- y'= 1-08=2}', 

- aj' + 3y = -0-34 = ^'. 
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axLd their solution gives 

51»' = 32^'+135' + 11(7'+ 9Z>' + 3^ = -9-73, 
61i' = UA' + 265' + 22(7' + 181/ + QE' = 706, 
etc., etc., 

from trhich 

»' = -0-19, «''=014, «' = 005, as' = 0-43, y = 005; 
and adding these to the assvuned approximate values, we have 
S = 573 -019 = 672-81, 
« = 574 + 014 = 57514, 
tt = 742 + 005 = 74205, 
« = 746 + 0-43 = 745-43, 
3jr = 320 + 005 = 32005, 
which are the same as obtained by the longer method. 

11. Adjustment of Angles taken at a point. At a station 
O there were measured the following horizontal angles, each 
i>eing the average of an equal number of readings. 

Fig. 6. ' 




^05 = 58»56' 42", 
AOD = 76 43 6, 
BOC^U 14 15 , 
BOD = 17 46 26, 
BOE^U 14 17, 
COD= 4 32 7, 
DOE^ie 27 54, 
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and it is required to find the most firobable yalaed 6t all the 
angles AOB, BOG, etc In the solution of sach a problem 
it will be found most conyenient to take as the unknown 
quantities the direction angles AOB = (, AOG '^x, A OD » y, 
and A OE = z. Then the ohserration equations will be 

«=58*56' 42", 
y = 76 43 6, 
«-l = 13 
y-< = 17 
«-<»34 
y— «= 4 
»-y = 16 

In order to avoid large numbers in the calculation* let us, 
as in the previous example, assume from the observations 
approximate values of the angles and desgnate bj f, ^, etc., 
the corrections to be apjdied to those approximate values. 
Thus if we place 

« = 58*56' 42' +«', 
«=72 10 57 +x', 
y = 76 43 6 H-y', 
«=:93 10 58 +z\ 
' Then, by substituting these in the observation eqoationd 
above, we get the simpler forms 

r =0, 

y=o, 

«'-«' = 0, 

y'-<'=2", 
«'- ^=1, 

y'-«'=-2, 
«'-y'=:2. 
in which the numerical terms include seconds only. From 
these we form the normal equations (Arts. 33, 34) 
4<'- 0^- j^- «' = -3", 
-^ + 2x'- y' = 2, 
-«'- a;' + 4y - «'=»-2, 
-f _ y' + 2a'- 3. 
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4»<1 by their solution find 

t'==-.0"l, aj' = rO, y' = 01, and «' = l-5; 
hence the most probable values of the direction angles are 

05 = 72 10 58 'O^AOC, 

y-76 43 6-9-^02>, 

« = 98 10 59 ^o^AOE, 

from which by simple subtraction we may find any required 
angle as BOCj GOD^ etc. (Strictly speaking, this problem is 
a case of conditioned observations, Chap. IV. By taking, 
however, a limited number of ui)kjiown quantities, the idea 
of condition need not enter the work until the final operation 
of Reducing any angle from tl^e direction angles by subtrac-^ 
tipn ; th^n, of pourse, the whol§ must equal the sum of its 
parts.) 

Problems. 1, To deten^ine the elevations of two points 
A and B above a datum 0, measuremients were made which 
give, A above = 12-3 feet, B above = 27 feet, and A 
telow B =« 14-1 feet What are the most probable elevations 
dfj.and£? Ans. 5 = 26-8 feet, etc. 

2. In the preceding example of the adjustment of 
XkngleSy what 9f e th^ best values if the sixth and seventh 
measwemeiJtQ wer§ not taken | 

Arts. AOC=-72P10'5T,eto. 



Observattona of Ufi4qual WeighU 

36. Suppose that four observations give each time the 
equa^on 

^nd that five other observations give each time 

The most probable values of 8 and t will be determined 
by writing the first equation four times, the second five 
times, and then by Art, 34 forming and solving the normal 
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equations. Considering for convenience only the unknown 
quantity 8, we may write the equations 

a^8 + ^j = 0, 4 times repeated, 

a^8 + ^j = 0, 6 times repeated. 
Forming from these the normal equations for « (Art. 33), 
we have 

4a, (a^ + N^ + ia^ (a^s + N^ = 0. 

If, instead of writing the observation equations 4 and 6 
times respectively, we should multiply the first by ^4 and 
the second by ,/5, giving 

»/Ea^8 + JEN^^0, 
and then considering each of these as repres^ting a single 
observation, form the tkor'mal equation for 8 by multiplying 
the first by >/4 a^, the second by ^6 a,, and adding the results, 
we would have 

4a,(a,^-fJV;)+5a,(a^ + -Kg=iO, 
which is the same as that found before. 

To solve observation equations of unequal weight, we 
have then the foUowing : — Multiply each observation equation 
hy the square root of its weight; then from tiiese form and 
solve the normal equatioriS asf beford* 

For illustration, let us suppose that four observations 
upon three unknown quantities have given the equations, 

os + y + z=^ 6*5 with the weight 3, 

aj= 0053 ....3, 

y« 0-003 :...3, 

«=-0043 1, 

that is, the last equation is the result of only one measure- 
ment, but each of the others of three. Multiplying each 
equation by the square root of its weight, we have 

JSx-hJSy-^JSz^ 5-5 J% 
JSw^ 0-053^3, 
JSy-^ O'OOSjS,, 
i5 = -0043. ^ 

Jigitized by Google 



UPON SEVEBAL QUANTITIES; 55 

Multiplying each equation in which x occurs by its co- 
efficient, viz. ,^3, we form by addition of the results the 
normal equation for a?; thus we have 

6aj + 3^+3^ = 16-659, 
3aj 4- 6y + 3« = 16-509, 
3i» + 3jr + 4« = 16-45'7, 
which being solved give 

a; = 0-9675, 3^ = 0-9175, ^» = 27005, 

which is the best system of values obtainable from, the 
observations. 

37. To adjust observations of unequal weight we have 
then, to write the observation equations, multiply each by 
the square root of its weight and from the reduced equations 
derive and solve the normal equations. The following 
examples will render the whole process clear.. 

I. Level Lines, The nine observations of Example 1, 
Art. 35, upon the elevations of the points, 8, T, etc. are of 
unequal weight. The least trustworthy is No. 9, because 
it is not known that mean tide at Baltimore is the same 
as the mean surface Of the ocean, and we call its weight 1. 
Nos. 3 to 8 inclusive are ordinary railroad levels and may 
with reference to No. 8 be given a weight of 4 Nos. 1 
and 2, being the result of carefully conducted government 
and canal levels extending over many years, are the most 
reliable of all, and we give them a weight of 25. The 
observation equations are the ^ame as before ; multiplying 
each by the square root of its weight, we have 

55 = 2865-40, 
6«-5i= 13-00, 

2^=1150-54, 
2w-2i= 334-66, 
2^-2w= 7-60, 
2x-2t=^ 340-56, 
2;B-2y= 85000, 

2y= 639-82, 
y= 319-75. ^ . 
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From tHese we form the normal equations "by tKe usual 
method, 

505-26^ « 14262-00, 

-255 + 37*-4u- 4ix = 1015-64, 

- 4t+8w- 4a? « 65412, 

- 4t-4w + 12a?-4y= 2396-32, 

- 4a!+9^ = -lQ0-61. 

Before solution these may be simplified by dividing the 
first by 25, and the third and fourth by 4, The resulting 
elevations are 

« = 572-98, 
^==575-48, 
w« 742-36, 
. ^« 745-72, 
y=: 320-25. 

A comparison of these values with those deduced in 
Art. 35 wm be interesting to the reader, as showing the 
influence of the weights and the closer agreement with 
Mr Gardner's values. (See also Art. 40.) 

II. Adhetment of Angles. In order to determine the 
angles AOn and £0V, Fig. 6, the instrument was set with 
its zero point in the direction OM and the following ob- 
servations were taken. 

Fig. 6. 




MOA « 46° 53' 29"-4, weight 4, 

MOB^ 83 14 36 -3, 16, 

JfOa=135 27 11-7, 9; 
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the first being the mean of 4 readings, the second of 16, and 
the third of 9. What are the most probable values of MO A, 
JLOJB&adBOCl 

Let X, y and « represent the required augleis, and to 
avoid using large numbers place 

«? = 46'53'30" + aj', 

ys^Se 21 7 +y; 

i5 = 52 12 35 -|-|j'; 

then the observation equations gxe 

a' « - 0"-6 with weight 4, 

y + a?'p=-0 7 16, 

^'+y' + »'=-0'3 ..,,,,.,. 9, 

Applying the weights, forming and solving tha normal 
equations, we find 

aj'«-0% y^O"-a* and «'»:-0"-54, 

and hence 

x^Wh2l %r^\^MOA, 

y^36 21 7-24=^05, 

a;c=52 12 34<46*£(?C; 

are the adjusted results. 

In the formation ai.nd solution of normal equations involv- 
ing large numbers of observations, the computer will derive 
much assistance from the notation of Gauss, which we have 
given in Part II., Arts. 37 and 58. In all cases the solution 
should be effected by independent niethods,in order to check 
the accuracy of the worJf , 

Problem, Levels are taken to determiae the elevation of 
three points A, B and C above the datum 0. The mean of 
five measurements show -4 to be 3*426 feet above 0, the 
xaeao of nine ebow 5 to be 10*828 above -4, the mean of four 
.give G 2*471 above B: and lastly, one measurement gives B 
I3'762 feet above OL What are the most probable elevations 
of tiie three points ? Ana. A ^ 34283, etc. 
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Probable Errors and Weights, 

38. The preceding methods are sufficient for the adjust^ 
ment of any common indirect observations. We now conie to 
^he question, ^hat degtee of confidence can we place in tho 
values of the quantities deduced from such measurements, 
or, in other words, ,to what degree of precision have we 
attained ? This will be shown by the probable errors of those 
quantities (Art. 16). * 

Let n be the mimber of the measurements, having, the 
relative weights ^r^, g^, ^j, etc.: the number of obsetvati^fi 
equations will also be n. Let q denote the number of im- 
known quantities 8, t, u, etc. whose values are to be deter- 
mined, and Og, Gi, etc, their respective weights. When the 
most probable values of 8,t, u, etc. have been found by the 
solution of the normal equations, let them be substituted in 
the observation equations {not in those equations after multii- 
plication by the square roots of their weights); they will not 
reduce these equations to zero, but leave the residuals v^, t?,, 
V,, etc. > 

Let S^rt?' denote the sum gjV* •\'gjo^ + g^^ + etc., that is, 
the quantity formed by multiplying the square of each residual 
by the weight of the corresponding equation, and taking the 
sum of the products^ Then, as proved in Part IL, the pro- 
bable error of an observation of the weight unity is : 

(87) r = 0-6745y^, 

hence (Art. 30) the jJrobable error of an observation of the 
weight g^ IS -j=:-. The probable errors of the adopted values 
of 8, tj u, etc. are also 

(81) S..^, a-^.etc 

to determine which we must compute r, and the weights 
0„ Gty etc. The weights g^, g^, etc. will be given by the con- 
ditions of the measurements. If, as /in Art. 34, the observa- 
tions are of equal weight, we have only to make all ^''s equal 
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to 1. Th^ weighlbs G^,> G^, of the determined quantities, are as 
yet unknown : to find them Gauss' method, as explained i^ 

the next article, may be: employed. - 

• 

39. As demonstrated in Part 11.,* the weights of the 
values of the unknown quantities may be thus found. After 
each observation equation has been multiplied by the square 
root of its weight, and thus all reduced to the same uiiit of 
weight (Art. 36), let the normal equations be formed (Art. 33). 
These will be of the /orm 

4,« + 5j* + G^i^ + eta =* ^, 
(80) u4,« + B,<+(7,tt + etc. = jB, 
-4^ + ^g* 4- C;,t^ + etc. = (7^ 
etc., etc., 

the first being the normal equation for «, the second for <, the 
third for u, etc. ; A^, A^... B^, B^, etc. being numerical co- 
efficients of the unknown quantities, and At B, C the absolute 
terms. The solution of these equations will give 

8 = aj4 -t- oi^B + a^G + etc., 
(92) « = )8i^ + /3^S+i8j(7+et€., 
u = ry^A + 7,5 + y^G+ etc., 
cftc., etc., 

that is, eacih unknown quantity will be given in terms of 
the absolute terms A, By G, etc. with the numerical coeffi- 
cients Oj, a,, )8j, )8„ etc. Then the weigbt of « is — , th6 

1 . ^ . 1 ' 

weight of Hs -5- , the weight of w is — , and so on. 

Stated in words the method is : Plcu^ instead of the 
absolute terms in the normal equations, the letters. A, B, C, etc., 
and solve the equations. Then the weight of any unhnown 
quantity ew s is the reciprocal of the coefficient of the absolute 
term A in the value of 8, the weight of t is the reciprocal of the 
coefficient of the absolute term B in the general value oft, etc, ; 
it being understood that A is the absolute term in the. normal 
equation for 5, J3 the. absolute term in the normal equation 
for t, and so on, 
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Thus in Art. 35, Example 1, wo have the normal equa- 
tions 

28^ t ^A, 

etc., 
whose solution gives 

* «gY -d +7^-8 + etc., 

^tc. ; 
and hence the weights of the values pf s, t, etc. are 

^ ^ etc 
32* 26' 

Problem. Given the observation equations 

2* + t = 7 with weight 3, 

e + Qt-^6 , 1, 

«- t^2 .,.e4, 

to find the values ^nd w^i^hts of 8 and t. 

Aniy tf = 3 with weight 16*43, 

f=l ,14-53, 

40. To illustrate in full the determination of weights 
and probable errors of indirectly observed quantities, let us 
take the example of Art. 37, iiX which the observation 
equations are 

No. 1. « = 573-OS with weight 25, 

,., % ir^$m 2-60 ,,,.-,... 25, 

... 3. ««:575 27 4, 

.,. 4. t*^< = 167-33 4, 

... 5. /c^u=: 3-80 4, 

,..6. a:- f = 170-28 ......,.,,.„., 4, 

7- OJ'-y 's 425*00 .,.^ ,.,^. 4, 

8. y=^ 319-91 4, 



• •f 



9. y = 319-75 .^ 1 , 
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and the Dormal equations are « 

60« - 25i - 1426^-00 = A, 

-255 + S7i-4li- ix ±s 1015-64:«jiS, 
«. 4«+8m- 4k = 65412= (7, 
- 4<-4M+12aj-4y=^ 2396-32 = JD, 
- 4a? + 9y = - 100-61 =JS; 
and it be required to find the probable error of each ob- 
servation and the weights and probable errors of the values of 
t and X. 

Placing the absolute terms equal to A, 5, C, etc. and 
solving the equations by any algebraic method (the method 
of indeterminate multipliers is the shortest), we find 

' 37^ + 745 + 640+542) + 24^ ^^_^ 
t^ — ^34j = 57o48, 

36^ + 725+ 1710+ 2702) + 120^ ^..»-^ 

^ 1788 ^*^*^2- 

By the rule of Art. 39, the weight of t is the reciprocal of 
the coefl&cient of the absolute term B in the normal equation 
for i, or 

^ 1 1341 -^.^ 

1341 

likewise the weight Of (2? is the reciprocal of the coefficient 
of 2), or 

(Note. If the object be merely to determine the weight of 
t, it is evident that it is unnecessary to retain Ay C, D and 
E in the algebraic work ; they may be placed equal to zera 
So in finding the weight of w, it is only necessary to retain 
2> in the computation.) 

We can now find the probable errors. The values of the 
quantities «, t, u, etc. as given by the normal equations are 

^ = 672-98, 

t = 675-48 with weight 18-12, 

w« 742-36, 

X = 745-72 with weight 6 62, 

y- 320-25. 
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Inserting these in the observation equations, the re- 
mainders or residuals v^, v^, etc. are plajced in the third 
column below, their squares in the fourth, and the product 
of each square by its corresponding weight in the fifth. 



No. 


0- 


V. 


«*. 


gvK 


1. 


25 


0-10 


0-010 


0-260 


2. 


25 


•11 


•012 


•800 


3. 




•20 


•040 


•160 


4. 




•44 


•194 


•776 


6. 




•43 


•185 


•730 


6. 




•02 


•000 


•002 


7. 




•48 


•210 


•840 


8. 




•34 


•116 


•464 


. 9. 




•50 


•260 


•260 


■ 






Sflri;2=3-762 



The sum of the products g^v^, gjo*, etc. is then 3'762. 
Since there are 9 observations and 5 unknown quantities, 
n = 9 and j = 6. Then (Art. 38) the probable error of an 
observation of weight unity, that is of No. 9, is 



r = 0-6746, 



' = 0-635 feet, 



and the probable error of observations 1 and 2 is 
0-635 



5 



= 0-107 feet, 



and of those from 3 to 8 inclusive is 



0-635 
2 



0-317 feet. 



Also the probable errors of the above values of t and x are 
0-635 



i2. = - 



0-153, and iJ. = -^S = 0-248, 
^6-62 



7i8-12 
and hence we may write 
t = elevation of Cleveland datum = 575-48 ± 0*153 feet, 
X = elevation of Pittsburg depot =745-72 ± 0-248, 
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It is then, as far as these observations show, an even 
wager that 675*48 feet does not differ from the true eleva- 
tion of the Clevelaud datjim by 0153, and aJso an even 
wager that 745*72 expresses the true height of the Pittsburg 
depot; within 0*248 feet. The first is therefore much the 
more accurately determined, 

Problem, Find the weights and probable errors of the 
values of ^, % and y in the preceding example. 

Am. s = 572-98 ± 0116 feet, etc. 

Other Applicatiom, 

By making j = 1 in the formulae of Art. 38, they reduce 
to those of Arts. 24 and 30. Direct observations upon a 
single quantity are then only a particular case of indirect 
ones, and the methods of the present chapter are suflBcient 
for their complete adjustment. Thus if there be only one 
unknown quantity z, and the measurements be made directly 
upon it, the observation equations are 

The normal equation formed from these will be 

and this gives at once the rule of the average. By Art. 39 
the weight of the value of z will be n, as in fact is implied 
in our .definition of weight (Art. 27), 

41. It is often the case that a quantity is measured 
in several parts, then of course its most probable value is 
the sum of the adjusted values of its parts. Thus, if we 
have measured by independent waya three parts of a base 
line, and find for them the values z^, «,, z^j the value to be 
taken for the whole line is 

Z=z^ + z^ + z^. 

If r^, r, and r, are the probable errors of z^j z^ and z^, the 
"e error of ^is 



probabh 



(100) B^Jr;' + r^''^r,\ 
In like manner, if a quantity Z is equal to any simple 
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function of the sum or diflterence of independently observed 
quantities, z^, «,, ^3, etc., that is, if 

-Z'=*i»i±^i±«,±etC, 

then having found the probable errors r^, r^ etc. of 0^, r^, etc^ 
the probable error of ^is given by the relation 

(100) fi' = r,* + r,* + r/+etc. 

As an example illustrating this, let us take a case of 
levelling. In order to determine the difference of level 
between two points A and B, a level was set up halfway 
between them, and 20 readings taken on rods held at those 
points with the following results : 

Rod at A. JW at B. 

7 readings gave 7'229 feet> 3 readings gave 9'806, 

8 7-230 ... 12 9-807, 

5 7-231 ... 6 9808. 

What is the most probable difference of level between 
the two points and the probable error of the determination ? 

The general mean (Art. 27) of the readings at A is 
7-2299, and of those at 5, 9-8071 feet. Hence the difference 
of level is 

9-8071 - 7-2299 « 2-5772 feet. 

To find its probable error, we have to find the probable 
errors of the two general means by Art. 30. For the mean 
7-2299, we find 

tg'(? = 00000118, w = S, G = 20, hence J2, = 0-00637, 

and for the mean 9*8071, we have 

Sflw* «: 0*0000078, n = 3, G^ := 20, hence R^ == 0-00031. 

Then from the above principle the probable error of 
the difference 2*5772 is 

70-00037' + 0-00031' = 000048. 
Hence the adjusted results are 

Reading at ^ = 7*2299 ± 000037, 
Reading at jB = 9-8071 ± 000031, 
Diff. of level = 2*5722 ± 000048. 
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Problem. The north declination of a star is 

S=:19'30'14"-8 

with a probable error of 0""8. The zenith distance of the 
same star is observed f = 21® 17' 20"'3 with a probable error 
of 2"-3. The latitude of the place of observation is ^ = S + f^ 
What is its probable §rror ? 

Ans. <t> = 40® 47' 35"-l ± 2"-44. 

42. If we have two quantities Z and z^ connected by 
the relation 

in which -4 is a constant, and if by measurement we find 
the value of z^ and its probable error r, the probable error of 
Z is given by 

(101) B^Ar^. 

Thus, the circumference of a circle is 3*1416 times its 
diameter. If we measure the latter and find its value 
z^ = 1000 ± 02, the value and probable error of the former 
will be Z= 3141-6 ±0-63. 

So in general if we have a quantity Z whose relation 
to the independent quantities z^, z^, z^, etc. is given by 

Z:=Az^ + Bz^ + Cz^ + etc., 

and if we find the values of z^, z^, etc. and their probable 
errors r^, r^, r^, etc. the probable error of Z is given by 

(102) J2' = A\^ + ffr^ + Cr^ + etc. 

We have here given merely the statement of these im- 
portant relations ; the proof is presented in full in the 
second part of this work. 

Miscellaneous Problems, 1. A chronometer is rated at 
a certain date and found to be 9™. 12"'3 fast with a probable 
error of 0''3. Ten days afterwards it is again. rated and 
found to be 9"*. 21'-4 fast with the same probable error. 
What is the probable error of the mean daily rate ? 



The rate in the whole interval is 

gm 21H - 9m 12-3 = 9«1, 
M. 
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with a probable error (Art. 41) of 



^0-3"+ 0-3' = 0-42. 
10 
= 0-042. 



9-1 
The mean daily rate is then -^^ = 0-91 and its probable 

error (Art. 42) is 

0-42 
10 

The clock gains then, daily, 9-1 + 0-042. 

2. Given the observation equations (all of equal weight) 

2a?- y+ ^= 3, 

32?+ 3y- « = 14, 

4a? + y + 4^ = 21, 

-5a? + 2y + 3«= 5, 

to find the best values of a?, y, and z, and their probable 
errors. 

Ans, 0?= 1-916 ±0026, 

y = 3-551 ± 0052, etc. 

3. A block of cast iron weighing 100 lbs. rests upon 
a horizontal table also of cast iron. A horizontal force is 
applied to the block and it is observed that it begins to 
move when the force is 16-5 lbs. If the probable error in 
the determination of this force is 0*5 lbs., what is the pro- 
bable error of the coefficient of friction fi ? 

Ans. fi = 0155 + 0005. 

4. The following levels were taken to determine the 
elevations of five points T, U, W, X and Y above the 
datum : 

T above = 115-52, X above W= 632-25, 

CT r= 6012, X r=21101, 

U 0=17704, Y Z7= 59612, 

W. r = 23412, Y Tr= 42718. 

W 17=17100, 

What are the adjusted elevations ? 

Ans. r= 115-61, CL^ 176-95, etc, 
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5. An angle is measured by a theodolite giving the 
value 37^ 16' 13" with a probable error of 5", and also by a 
sextant giving 37* 16' 10*' with a probable error of 3 "• What 
is the adjusted value and its probable error ? 

Further examples illustrating the adjustment of indirect 
observations which arise in physical investigations are given 
in Chapter V. 
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CONDITIONED OBSERVATIONS. 



43. All the observations thus far considered have been 
independent, so that an error committed in one measurement 
has had no connection whatever with those arising from the 

Fig. 7. 




others. Thus in measuring the angles A OB and BOG, whe- 
ther it be done directly or indirectly (Art. 21), each reading 
or observation has been entirely independent of those pre- 
ceding and following, and a variation in the value of A OB 
does not necessarily require a corresponding one in BOC. 
But if a third observation be made upon A OC, the values of 
A OB, BOO and AOC are no longer independent, for since 
the whole must equal the sum qf its parts, the relation 
AOG=AOB-h BOG exists, and it is impossible to suppose 
the value of BOG to vary without a corresponding variation 
in the value of AOG or AOB, 

We have then a new class of equations, viz. conditional 
equations, which must be exactly satisfied by the values 
adopted for the unknown quantities, since they are the 
expression of axioms and rigorous laws. The number of 
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these equations is always less than the number of unknown 
quantities, for if they were as many in number, the values 
of the latter would be determined by their solution without 
the necessity of measurement. From these conditional equa- 
tions and from the observation equations the values of the 
unknown quantities are to be found*. 



Observations of equal weight. 

44. The conditional equations being less than the 
number of unknown quantities may be satisfied in various 
ways. The observation equations cannot however be all 
exactly satisfied, and hence, as in the preceding chapter, the 
best system of values must be found. This is done (Arts. 16, 
33) by making the sum of the squares of the residuals a 
minimum. Hence we have to determine the values of the 
unknown quantities in such a way that they shall he the most 
prohaMe values for the observation equations and shall at the 
same time exactly satisfy the conditional equations. 

Thus in the case of the last Article suppose that the mea- 
surements give the values 

^05 = 36* 24' 30", 

BOG ^47 52 20, 

-40(7=84 17 10. 

If the value of AOC were exactly equal to the sum of the 
other two values, the quantities would need no correction. 
As they stand, however, the results are discordant and must 
be adjusted. Let the unknown values of AOB be x, ot BOG, 
y, and of AOG, z, then we have the observation equations 

a; = 36' 24' 30", 

y = 47 52 20, 

i5 = 84 17 10, 

♦ In most books npon tMs subject, the term *< equations of condition** 
is applied indiscriminately to both of these very distinct classes, and is a 
cause of some perplexity to the student. The excellent distinction of the 
Germans, Beobachtimgsgleichwng and Bedingtmgsgleichung, ought certainly 
to oome into use.; 
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and tlie rigorous conditional equation 

and we must determine the values of a?, y and z so that they 
will exactly satisfy the latter and, at the same time, be the 
most probable values for the former. If for ^, in the last 
observation equation, we place its value a? + y from the con- 
ditional equation we have 

a? = 36" 24' 30", 

y = 47 52 20, 
a? + y = 84 17 12, 
or simply three observation equations, each one of which is 
now independent of the other, and which can be solved in 
the usual way (Arts. 33, 34). Multiplying the first and 
third by 1 and adding^ them, we have the normal equation 
for Xy doing the same for the second and third we have the 
normal equation for y, or 

2a?+y = 120" 41' 40", 
aj + 2y = 132 9 32. 
The solution of these gives 

a; = 36' 24' 361"* y = 4r52'26r, 
and hence 

«=:a.+y=84« \r SJ". 

The process here exhibited has already been illustrated 
in some of the examples of the preceding chapter, in which, 
instead of stating the conditioned equations, we have simply 
written them as observation equations, omitting the addi- 
tional unknown quantity. (See Art. 35. Example 2.) We 
may therefore give the following as the statement of a 
method for adjusting conditioned observations of equal 
weight. 

1st. For every observation, whether direct or indirect, 
write an observation equation (Art. 32). Let n equal the 
number of such equations, and q the number of unknown 
quantities involved. 

2nd. For each rigorous condition, write a conditional 
equation (Art. 43), and let p equal their number.. ooIp 
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3rd. From the conditional equations find the values of 
p unknown quantities in terms of the others, and substitute 
them in the observation equations. There will then be n 
such equations containing o'—p unknown quantities, each of 
which will represent an independent observation. 

4th. From these observation equations, form and solve 
the normal equations (Arts. 33, 34); the resulting values 
will be the most probable. Then the values of the remain- 
ing unknown quantities may be directly found from the 
conditional equations. 

As a fuller illustration we choose the following. At the 
points A, B and 0, there are measured the angles 



« = 9r 27' 40", 

< = 43 52 60, 

tt = 44 39 60, 

sum = 180 20, 



y = 20^ 15' 10", 
« = 64 65 10, 



Fig. 8. 




subject to the geometrical conditions that at the point A 

and that in the triangle ABO 

« + < + w = 180^. 

To avoid the use of large numbers we adopt the method 
employed to some. extent in the last chapter of assuming 
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approximate values for the angles and regarding the cor- 
rections to be applied to those values as the unknown quan- 
tities J thus if we place 

8 = 91' 27'+ s', 

t = 43 52 + *', 

etc. 

and substitute the values in the above expressions, we have 

s' = 40", w'+y' = / + 60, 

<' = 50, «' + «' + «'« 120, 

u = 60, 

y' = 10, 

z' = 10, 

in which the numbers represent seconds only. From the 
two conditional equations we take the values of any two 
unknown quantities in terms of the others, for example, 

z' = u'+y' -60, 8'=120-u'-i, 

and substitute them in the observation equations, giving 

• «' + «' = 80", 

t' = 50, 

u = 50, 

2,' = 10, 

y' + w' = 70. 

From these we form the normal equations 

2<'+ u =130, 

^ + Su'+ y' = 200, 

tt' + 2y= 80, 
whose solution gives 

«' = 41"-25, m' = 47"-5, y' = I6"-25. 

Then from the conditional equations we have the remaia- 
ing quantities 

«'=3"-75, «' = 31"-23. 
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These are the corrections to be added to the appro^^imate 
values assumed, that is, the number of seconds to be taken* 
instead of the observed values. Hence the adjusted results 
are 

5= 91^27' 31"-25 

<«= 43 52 41 -25 w = 44' 39' 47"-50 

^= 44 39 47 '50 y = 20 15 16 '25 
sum = 180 «= 64 55 3 -75 

which exactly satisfy the two geometrical conditions. 

Problem, The three angles of a triangle were measured' 
as follows : 

ic =« 98*^ 17' 22", y = 70' 9' 56", and z ^ IV 32' 62". 
What are the adjusted values ? 

Ans. a; = 98' ir 18"-67, etc. 

45. The preceding method is general, and may be 
applied to any number of observations subject to any num- 
ber of conditions. Although the simplest in theory, it is not 
always so in practice, particularly when the number of con- 
ditional equations is large. The process in most common use 
is "Gauss' method of correlatives," which we will now proceed, 
to state and exemplify, referring the reader to Part II. for 
the proof and the fuller algorithm of the method. 

1st. Let the observations be adjusted by the methods 
of Chapter II., without reference to the conditional equations, 
and the results be called the measured values. 

2nd. For each rigorous condition write a conditional 
equation (Art. 43). If the measured values exactly satisfy 
these they will need no correction; if not, a further adjust- 
ment is necessary. 

3rd. Let /, t't u\ etc. be corrections to be applied to 
these values to make them satisfy the conditional equations. 
Let the conditional equations expressed in terms of these 
corregtions be of the form 

a/ + or/ + ffgU + etc. = jy', 
(104) i8/+y3/+y33w' + etc. = iV'', 
7i«'+7/+78«*' + etc. = j?^% 

in which a, ^y 7, etc. and N denote known constants. 
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4th, For each one of the unknown quantities s\ t\ u\ etc., 
write an equation of correlative containing as many new 
auxiliary unknown quantities as there are conditional equa- 
tions. Let them be of the form 

a,K,^fi,K, + 7A + etc. = /, 

(1 10) a,K, + I3,K, + y,K, + etc. = t\ 

a,K, + l3,K, + y,K, + etc. =w' 

etc. etc. 

in which K , K^^ K^, etc. are the new auxiliary unknown 
quantities, K^ having the same coefficients as in the firet con- 
ditional equation, K^ the same as in the second, eta; so that 
the first vertical row of coefficients corresponds to the first 
horizontal row in the conditional equations. There will be 
as many equations of correlative as there are unknown quan- 
tities «', t\ u\ etc. 

5th. From the equations of correlative let the normal 
eqtuitions (Art. 33) be formed, thus 

(«i* + «2* + ^^^') ^x + (« A + « A + e*c.) K^ + etc. 
(109) = a^s' + a/ + etc. = N\ 

(a,/8, + a^, + etc.) Z", + (/8,» + ^8,* + etc.) iT, + etc. 

= /8/ + i8/ + etc. = JV", 

which will be as many in number as there are conditional, 
equations, and whose absolute terras will be N\ N'\ etc. as 
given by the conditional equations. Solving these we find 
the values of K^, K^, etc. 

6th. Substituting these values of K^^ K^, etc. in the 
equations of correlative, we find the values of 8\ t\ u\ etc. 
which will exactly satisfy the conditional equations, and be 
the most probable corrections to the observed quantities. 

This process will be better understood, and its simplicity 
over that of the preceding Article be seen, by a consideration 
of some practical examples. We take first that of Fig. 9, in 
which there are only two equations of condition. 
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is 



At the points A, B and C there are measured five angles, 
each measurement having the same weight, 
5= 91*27' 40" 

t^ 43 52 50 tt = 44*39' 50" 

u^ 44 39 50 y = 20 15 10 

« = 64 55 10 



sum =180 20 



Fig. 9. 




A C 

The rigorous geometrical conditions Q,re 
w + y = «, 
5 + ^4- u= 180*, 
and since these are not satisfied by the measured angles, an 
adjustment is necessary. We assume the measured values 
as approximate, and take «', t\ u\ etc. as corrections to be 
applied to those values, then the conditional equations are 

u'+y'-V=^ 10" (1), 

«'+^' + w' = -20 (2). 

We next assume two unknown auxiliary quantities K 
and K^y and for each unknown quantity write an equation of 
correlative, thus 

^. = »' (3), 

K, = t' (4). 

K,-\-K,^^ (5). 

K, =y'...-. (6), 

-K, ^z' .^ (7), 
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the coefficients of K^ and K in (3) being the same as those 
of s' in (1) and (2), viz. and 1 : in (5) the same as those of v! 
in (1) and (2), viz. 1 and 1: in (7) the same as those of z in 
(1) and (2), viz. — 1 and 0. From these we form the normal 
equations (Art. 33) for K^^ and K^ ; viz, 

3ir, + ^, = tt'+y -«'= 10, 

^, + 3ir, = 5' + ^+u=-20; 

from which we have 

Z; = 6"-25 and ^, = - 8"-75, 

Substituting these values in (3), (4), (5), etc. we have 

/ = -.8"75, ^=-875, i^' = -2'50, ^ = 6-25, «'=-6-25. 

Applying these corrections to the measured angles^ they 
become 

8 = 9r 27' 3r-25 

t'= 43 52 41 25 t^' = 44^ 39' 4r''50 

^'= 44 39 47 -50 y = 20 15 1625 ' 
sum =180^ 0- / = 64 55 375 

or the same as deduced from the longer process of the pre- 
ceding article. 

46. As a second example we choose a case occurring in 
the common practice of every engineer, viz. the ' 

Adjustment of the angles of a Quadrilateral. In order to 
determine the distances WZ,ZY, etc (Fig. 10), the base line 
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TFiy was measured and nine angles of the quadrilateral ob- 
served as follows {X denoting the large angle at the comer 
JT, and X^ and X^ the small ones at the same comer, as shown 
in the sketch): 

F = 106° r 30" X = 66' 34' 9" 

. Xj= 36 34 21 r, = 49 17 23 

Z,= 37 18 12 Tr,= 64 8 34 



sum = 180 3 sum = 180 6 

Z = 84 7 18 

F,= 41 58 47 

F,= 53 53 50 

/ sum = 179 59 55 

In each of the four triangles which make up the figure, 
the sum of the measured angles differs from 180^ It will 
also be seen that at the comer TT, the sum of the two small 
angles does not equal the large one, thus 

Tf,= 4r 58' 47" 
F,= 64 8 34 



sum = 106 7 21 while W = 106° T 30" ; 

and the sum of the four quadrilateral angles is not 360^ 
The problem before us is, to adjust these angles so that in 
every triangle the sum of the three angles shall be exactly 
180^ so that at every corner the large angle shall equal the 
sum of the two small ones, so that the sum of the four qua- 
drilateral angles shall be exactly 360^ and further, so that 
the adjusted values shall differ from the measured values by 
the least possible amounts. 

We assume the above measured values as approximate 
ones, and take as unknown quantities the corrections to be 
applied to those values. Thus, at the corner W we designate 
the measured value of the large angle by W, the correction 
to be applied to it by w, and the required adjusted value by 
XWZ, so that 

XWZ = W+ w, ZWY=: W, + w^, WZX= Z^ + z^, etc., 
and the problem is to determine the corrections w, w^, z^, etc. 
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In order to avoid too many equations we select any corner 
as W, and take the three triangles WXZ, ZWY, and XFIT 
which meet at that point as the three triangles for correction : 
the angles of the fourth triangle XYZ will be found by 
simple addition and subtraction as soon as those of the other 
three are determined. 

The conditional equations for these three triangles are 
next to be stated. In the triangle WXZ the corrections to 
be applied to the measured angles are w, a;,, and z^* aod since 
the measured angles add up to 180® 0' 3", the algebraic Bum 
of the three corrections must equal — 3" ; and hence the first 
conditional equation, 

w + a?i + «, = -3 (1). 

Next, in the triangle ZWY the corrections are «, -i^j, and y^, 
and they must exactly balance the discrepancy between 180* 
and the sum of Z, TF^, and F,, or 

'^ + ^i + y, = 5 (2), 

Likewise in the triangle XYW, we must have 

» + yi + ^« = -6 (3> 

Also at the comer W the corrections w, w^ and w^ must 
exactly balance the discrepancy of 9" between w and the 
sum of w^ and w^, or we must have 

Wj + t^/g — w = 9 (i). 

These are the four geometrical conditions which the cor- 
rections must exactly satisfy. From these four equations we 
must find the values of the nine unknown quantities in such 
a way that they shall exactly satisfy the conditional equa- 
tion, and be the most probable system of corrections. This 
may be done either by the method of Art. 44, or by that of 
Art. 45. By the former we could write nine observati(m 
equations (viz. -m; = 0, 0?^ = 0, «, = 0, etc.), then from the five 
conditional equations find the values of five unknown quan- 
tities in terms of the other four, and substitute them in the 
observation equations, then from those nine equations deduce 
the four normal equations, and by their solution find the 
required corrections. We shall however follow the " method 
of correlatives," since it is by far the easier and shorter. 
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Our conditional equations are then the following : 

w +a;, +^, =-3 (1), 

^+^i + y«= 5 (2), 

X +yi +t^?, = -6 (3X 

w^ + w^ — w = 9 (4). 

We assume four auxiliary unknown quantities K^^K^^K^^ 
and K^t and for each unknown quantity write an equation of 
correlative, thus 

+ Z; ^K,^w (5), 

+ K, •\-K,^w, (6), 

+ ^, =«, 

the coefficients of K^ being the coefficients of the correspond- 
ing unknown quantities in equation (1), the coefficients of K^ 
being those of the unknown quantities in equation (3), and 
so on. From these equations we form the normal equations 
for K^y K^, etc., viz. 

SK^ - Z, = -3=«; + «, + 0„ 

3K^ + ^^ = + 5, = etc. 

3^3+ ^, = ^6, 

-K^-¥K, + K, + SK, = + 9, 

whose solution gives 

Z;=: 0-389, Jr, = 0-278, Jr3 = - 3*389, Z; = 4167. 

Substituting these in the equations of correlative (5), (6), 
etc. we have ^ i 
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w^K,-K, 


ss . 


- 3"-78, 


w, = K, + K, 


= 


4-44, 


w,= K, + K, 


= 


0-78, 


«=^. 


= - 


- 3-39, 


x, = K, 


s 


8-39, 


y.=^. 


= ■ 


- 3-39, 


y. = ^. 


= 


0-28, 


z^K, 


=s 


0-28, 


z,= K, 


= 


0-39. 



Applying these corrections to the meafiured angles, they 
become 



Tr+w, = 106' T 26"-22 

X, + a;, = 36 34 21 -39 

Z, + z^ =» 37 18 12-39 

sum = 180 

^ + « = 84" 7' 18"-28 

Tr, + «;^= 41 58 51 44 

r. + y, = 53 53 50 -28 

sum = 180 



X + oj = 66° 34' 5"-61 

r, + y, = 49 17 19 -61 

F, + M>,= 64 8 34 -78 

sum = 180 00 



F, + w, = 41° 58' 51"-44 
T^'. + to,= 64 8 34-78 
sum = 106 



7 26 -22 
= TT+w, 

The three triangles WXZ, ZWY, and XYW are hence 
exactly adjusted, as are also the angles at the comer W. 
The remaining angles are obtained by simple addition and 
subtraction of those already fouud, thus 

r4-y=(F, + yJ + (r.+y,) = 10rir 9"-89 
Z^^z^^{Z ^z) --(Z^ + z^i^ 46 49 5-89 
X,+ a;,= (X+a;) - (X, + a; J =_29_59_44j22 

sum = 180 -0 

and the angles of the quadrilateral also add up exactly to 
360^ 

The above adjustment is sufficient for cases arising in 
common engineering practice. In extensive triangulations, 
however, where the sides are many miles in length, a fifth 
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conditional equation is necessary. This arises from the fact 
that we have not considered the relation existing between 

Fig. 10. 




the sides and angles, so that if with the above adjusted 
angles we were to compute from the base WX the side ZY, 
through the triangles WZX and ZYX, and again through 
the triangles W YZ axid WZY, the two results would not 
exactly agree. In such triaujgulations, therefore, as are made 
in our Coast Survey, besides the conditional equations 
between the angles, others called side equations are intro- 
duced. We proceed to illustrate the development and ap- 
plication of this relation between the sides to the case of the 
above quadrilateral. 

Designating the adjusted angles by the same notation as 
before, we have the condition that in each of the three tri- 
angles TFZZ", WZY md WYZ, when adjusted, the sides 
must be proportional to the sines of their opposite angles, or 

WX_ BinWZX WZ yi nWYZ WYs mWX Y 
WZ^sm WXZ' WY^sux WZY' WX^sin WYX' 

multiplying these equations together, member by member, 

we have 

sin WZXsm WYZsmWXY _ ^ 
sin WXZ sin WZYsm WYX ' 

which is the fifth equation of condition. To adapt it to 
numerical computation, we write it . 

sin TTXr sin WYZ sin WZX = sin WXZ sin WYZ sin WZY. 

M. & 
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Now WXY, WYZy etc. are the adjusted values of the 
angles, or 

TFXr=X+a?, TrrZ= r,+y„ etc.; 

hence the equation is 

sin (X+ x) sin ( F, + y^ sin (Z^ + z^ , 

= sin (Xj + ajj sin (F, +y^ sin {Z+ z), 

in which x, y^, &c. are the unknown corrections. Applying 
logarithms, the equation becomes i 

log sin (X+ a?) + log sin (r, + yj + log sin {Z^ + -sr^) 

= log sin (Xj + x^ + log sin (Fj 4- yj + log sin (Z+ ^). 

Now the corrections x, y,, ;8rj, &c. are very small, each 
being only a few seconds of arc. Therefore very nearly 

log sin [X-\-x) = log sin X + 0? log. difif. 1", 
log sin(F,+y8) = log sin Fj+y^log. diff. 1", etc., 

in which log. diflf. 1" denotes the tabular logarithmic differ- 
ence for 1" corresponding to the measured values X, F^, etc. 
Hence, inserting these values and transposing, our equation 
becomes , 

a;, diff. l"+y,diffi 1" +;?,diff l"-a;,diff I'-y^diff. l"-^,diffi 1" 

= log sin Xj + log sin F^4- log sin Z— log sin X— log sin F^ 

— log sin Z^, 

The logarithmic sines of the measured angles and their 
tabular differences may then be taken from a table and 
arranged as below : 



Angle 


log. sin 


diff. r 


Angle 


log. sin 


diff. 1" 


f. 


9-77612924 
9-87967913 
9-99771028 


0-00000284 
0-00000181 
0-00000021 


X 


9-96262639 
9-90739060 
9-78249744 


0-00000091 
0-00000154 
0-00000277 


sum 


=29-66261865 


sum 


=29-65261343 
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Substituting these values the equation is 

■ 0-00000091aj + 000000154.y, +,etc. 

= 29-65251865 - '29-65251343 ; 

or after multiplying by 10000000 to avoid decimals, 

91a? + 154y, + 277-^, - 284«r^ - ISly, - 21^? = 522 (5) 

which is the conditional equation, that any side computed 
from the measured base shall always be of the same length, 
whatever set of triangles be employed. 

The four conditions between the angles are the same as 
before ; and we have to determine the nine unknown cor- 
rections so that they shall be the best values, and at the 
3ame time exactly satisfy the five conditions, • 

w +aJi + ^a = -3 • (1)' 

2 +Wi + y2= 5 (2), 

a? +yi + «^2=-6 (3), 

10^ + 10^-10=^ 9 (4), 

91^ + 154y, + 277;^, - 2Six, - ISly^ - 21z = 522. ..(5). 

This we perform as before by the method of correlatives. 
Assuming five auxiliary unknown quantities, K^,K^, etc., 
the equations of correlatives are 

^ +K, + K, =t^„ 

+ Z3 + 91^,= ^, 

-fZ, -284ir,= a?,, 

+ K, -181^, = y,, 

+ ^^ +154ir, = y„ 

+ K^ + 277 K,^z^; 

the coeflBcients of K^ being those of the corresponding un- 
known quantities in the conditional equation (1), the coeflS- 
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cients of K^ being those of the corresponding unknown quan- 
tities in equation (5), etc. From these equations the normal 
equations for K^, K^j etc. are formed, viz. 

SJT, +^,+ 133ir,=: 5, 

sjf.+jr,- 90ir,=-6, 

- ir,+ K^^ ir,+3Z, - 9, 

-7^i + 133ir,- 90^3 +222584Z,= 522. 

Solving thetfe equations by any algebraic method, we 
have 

K, = 0-393, Z,-0-238, Z, = - 3-366, r,=4174,'ir = 0000855. 

Substituting these in the equations of correlative we get 
the values of the required corrections, viz. 

«; =Z;- iT, = -3"-78, 

w^^K^+ ir,= 4-41, 

l^?, = ir,+ K,^ 0-81, 

X =^8+ 91irB = -3 -29, 

fl?j = Z;-284ir,= 15, 

y, = Z3-181ir, = -3-52, 

y. = J5r.+ 154S;= 0-37, 

z =J5r,- 21^,= -22, 

«, = ir,+ 277Z,= '63. 

These values exactly satisfy the conditional equations 
(1), (2), (3), and (4), and almost satisfy (5). Applying them 
to the measured angles we have the adjusted values, which 
differ by only one or two-tenths of a second from those 
determined without the use of the side equation, but which 
conform equally well to the geometrical requirements of the 
angles, and which in delicate computations will give no dis» 
cordant results in the lengths of the sides. 

The student who verifies the above results will have 
acquired some idea of the processes involved in the adjust- 
ment of extensive triangulations, for the method here given 
is that employed by the computers of the U. S. Coast Survey 
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Office. As the number of triangles increases the number of 
conditional equations and unknown quantities rapidly mul- 
tiplies, and their proper solution requires great skiU and 
patience. In the U. S. Coast Survey Report for 1854* men- 
tion is made of the adjustment of a triangulation involving 
sixty-five conditional equations. In such cases, we ought to 
say, the work is reduced to a routine by means of the algo- 
rithm exhibited in Part II., Arts. 37 and 60, and the results 
are obtained independently by different computers. In fact, 
it would be well-nigh impossible to deal with such a large 
number of equations, without performing the numerical 
operations by a systematic routine, and employing frequent 
checks to test the accuracy of the work. 

As an exercise in operating with side equations we give 
the following problem, which requires no angle equations 
and only one side equation, which will be similar in form to 
that of the preceding example. The answers to the pro- 
blem may be seen in Sie K IS. Coast Survey Report for 1854^ 
Appendix, p. 82. 

Problem. In the quadrilateral WXYZ, Fig. 10, five 
angles were measured as follows : 




X=66'34.' 2" -745, 
Y, = 49 17 21 -310, 
r, = 53 53 54 -075, 
Z, = 46 49 5 -167, 
Z. = 37 18 11 -542. 
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What are the adjusted values, so that in computing any 
side from the measured base WX, it shall always be of the 
same length whatever set of triangles be employed ? 

Observations of Unequal Weight. 

47. A general method of adjusting conditioned obser- 
vations of diflferent weights follows directly from the pro- 
cesses of Arts. 37 and 44. It may be thus stated : 

1st. For each of the n observations write an observation 
equation (Art. 32) in terms of the unknown quantities. 

2nd. For each rigorous requirement write a conditUmaZ 
equation (Art. 43). 

3rd. From the conditional equations find the values of 
p unknown quantities in terms of the others, and substitute 
them in the n observation equations, each of which will then 
be independent (Art. 44). 

• 4th. Multiply each observation equation by the square 
root of its weight (Art. 37). 

5th. Form and solve the normal emudiona (Arts. 33, 34), 
thus obtaining the most probable values of the unknown: 
quantities. 

• For example, suppose that at the point 0, Fig. 11, four 
angles are measured as follows : 

Fig. 11. 
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AOB = 40' 52' 37" weight 16 
jB0(7= 92 25 41 „ 4 
(702)= 80 6 15 ^ 3 
D0A=-U6 35 20 „ 1 

sum = 359 59 63 

the angle AOB being the mean of 16 observations, BOG ot 
4, etc. 

We assume the measured angles as approximate, and 
take w, X, y, and z as corrections to be applied to those 
values. Then the observation equations are 

«; = with weight 16, 

a? = „ „ 4, 

y = o „ „ 3, 

« = „ „ 1, 

and the rigorous conditional equation is 
io + x + y+z^T. 

Taking from this the value of «, inserting it in the obser- 
vation equations, and multiplying each by the square root of 
its weight, we have 

4tt? = 0, 

2a? = 0, 

N/3y = 0, 

tr + ic + y = 0. 

From these the normal equations are 

I7w+ a;+ y = 7, 
w + 5itt+ y = 7, 

w+ aj + 4y = 7, 
from which we find 

i£; = 0"-27, a? = r-06, ,y = r-42, and hence « = 4"-25. 
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Adding these to the measured values, we have 

BOG = 92 25 42 -06 
COD =80 6 16 -42 
DOA = 146 35 24 25 



sum = 360 



Problem. To determine the angles of the triangle ABC 
the following measurements were taien. 



Fig. 12. 








A 




a 












ACB= 


= 91" 


27' 35" 


weight 


4, 










CBA = 


= 43 


52 42 


» 


9. 










MAC = 


= 64 


55 5 


» 


2, 










MAB = 


= 20 


15 15 


n 


2. 








What 


ftre the adjustec 


[ results ? 
















Ans. 


AGB= 


91" 


27' 33"-71, 












CBA = 


43 


62 


41-43, 


etc. 



48. As in the case of equal weighted observations, 
Gauss' method of correlatives may be advantageously em- 
ployed when the number of unknown quantities or condi- 
tional equations is greater than two or three. Although the 
process of Art. 47 is perfectly general and applicable to all 
cases, it is not in practice of so easy application as the 
method of correlatives, which is hence generally employed. 
The following is the statement of the process to be fol- 
lowed : 
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1st. Combine the observations by the methods of Chap- 
ters IT. and ILL, and find their most probable values, and 
their weights. 

2nd. Let s', ^, u', etc. be the best system of corrections 
to be applied to those values to make them conform to the 
conditional equations, whick expressed in terms of those 
corrections are of the form, 

ri04^ a/ + a.«' + a3Z.'+etc. = iV^', 

^^"^^^ 13,8' + /3/ + ^,u'+ etc. = JV^". 

3rd. Let g^, g^ etc. be the weights of the values to 
which d, t\ etc. are corrections, and for each of these cor- 
rections write an equation of correlativej thus 

. |^(«A + A^a + yx^3 + etc.)=a', 
(110) 1 ^a,K, + fi^, + y,K, + etc.) = t\ 

i- ^a^K, + fiJS:^ + y^, + etc^^u\ etc., 

in which K,, K^, etc. are new auxiliary quantities, K, 
having the same coeflScients as in the first conditional equa- 
tions, K^ the same as in the second, etc. 

4th. From these form the normal equations j for iT, Zl, 

11 

etc. (not regarding — , — , etc., as coefficients of jff^, K^, etc.). 

' 9i 9t 
Then the second terms of these normal equations will be 
N\ N'\ etc., as in the given conditional equations. 

5th. Solve the normal equations, and then substitute 
the resulting values of K,, K^, etc. in the correlative equa- 
tions, thus finding the best system of corrections. 

Although somewhat tedious to state, this method is of 
very rapid application. As a first example we choose a case 
involving but three conditional equations. 

At the points -4, B, 0, and i). Fig. 13, there are mea- 
sured the following seven angles. 
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Fig. 13. 




^* — ^^^ . 

\^ y B 

n 

BAG= So'' 2' 10" weight 2 
G^J5^= 40 12 17 „ 3 

A OB =104i 45 39 „ 5 

sum = 180 6 

i?a4 = 109' 52' T weights 
CAD^ 45 29 „ t) 

^2)a= 25 7 20 „ 1 

sum = 179 59 5ft 

Da5 = 145* 22' 15'' weight 7 
^CB = 104 45 39 „ 5 

DC4 =109 52 7 ^4 

sum = 360 1 

the angles of the. triangle ABC being 6" greater than 180^ 
those oi AGD being 4" less than ISO* and those around the 
point being 1" greater than 360 ^ 

We assume the measured angles as approximate, and 
take a, f, u^ etc. as corrections to be applied to BAG, GBA, 
etc. (see Fig. 13). Then the conditional equations are 

« + «+w = -6" (1), 

x+y +-?= 4 (2), 

U'hw-\-x=-l (3). 
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The weights of s, t, u, w, x, y and z are 2, 3, 5, 7, % 6 
and 1, corresponding to the weights of the measured angles 
to which they are correctious. The equations of correlative 
are then 

i( +K, + K,) = x. 

the coefficients of K^ being those of the corresponding 
unknown quantities in equation (1), those of K from those 
of (2) and so on. From these we form the normal equations 

(^l + l + l^jK^ + lK, ^x + y + z= 4. 

, |^. + |^. + g + i + ^)^.= « + - + «' = -l. 

or after reduction 

31 Z, + 6Z, = -180, 

34.^,+ 6J5r,= 96, 
28 Ki + 35,^, + 83 Z, = - 140, 
from which we have by solution 

Z, = - 5-601, £, = + 301, £. = -1058. 
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Substituting these values in the correlative equations, 
we have 

s' = -2"-80, « = -l-87, t* = -l-33, w=^-0'15, 

X = 0-49, y = 0-50, z = 3-01. 

Adding these to the measured values we obtain 

BAC=^ S5' 2' 7"-20 DCA^lOd' 52' 7"-49 



CBA= 40 


12 15 -13 


CAD= 45 





29 -50 


^C5=104 


45 37 -67 
0-0 


ADO'^ 25 


7 


23 -01 


sum = 180 


sum = 180 





-0 




Z)GB = 145"' 


22' 14"-85 








^CB=104 


45 37 -67 








DCA ^109 


52 7 -47 







sum = 360 0-0 
Which is the best system of adjusted angles? 

Problems. 1. In a quadrilateral WXYZ the four angles 
are measured as follows : 

Tr= 65* ir 64" with weight 5, 

X= 66 24 15-5 „ „ 10, 

r= 87 2 24-7 „ „ 12, 

Z =141 21 20- „ „ 1. 

What are the adjusted angles? 

Ans. F= 65* 11' 54"-84, eta 

2. In a spherical triangle XYZ the three angles were 
measured 

X=93* 48' 15"-22 with weight 30, 

r= 51 55 -18 „ „ 19, 

Z=34 16 49 72 „ „ 13. 

The spherical excess is 4 "•054. What are the adjusted 
angles? 
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The conditional equation is 

x+r-fZ= 180^ 0' ro64, 

or taking m, y^ and z as corrections^ it is 
a;-fy-f««-l"066. 
The correlative equations then are 

From which the normal equation is 

and hence the corrections are 

a; = -0-217, y = - 0-347, « = - 0*503, 

which applied to the measured angles will adjust them so as 
to agree with the geometrical condition for a spherical 
triangle. 

3. To adjust the percentages of the following chemical 
analysis of pig iron, so that the sum shall be 100. 

Carbon 2'30 per cent. — weight of determination = 1 

Silicon 1*36 „ „ „ „ 5 

Phosphorus 1*02 „ „ „ w 5 

Sulphur 0-61 „ « „ „ 6 



Iron 94-93 



« " " 10 

Manganese 0'81 „ „ „ „ 1 

Ana. Sulphur = 0*594. Iron = 94*11 eta 

4. A base line AG is found by one measurement to be 
1599 feet; by another measurement it is chained in two parts, 
and AB found to be 791, and BC 806 feet. Considering the 
first as of weight 2 and each of the others as of weight 5, 
what is the adjusted length? ^ i 
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Probable Errors. 

49. The determination of the weights and probable 
errors of conditioned observations follows directly from the 
method of Art 38, since all such observations can be reduced 
to independent ones; ■ 

Let n denote the number of observations, q the num- 
ber of imknown quantities involved, and p the number of 
conditional equations. Let ff^, g^, ffs*'^^* ^® ^® weights of 
the several measurements, and j/^, v^, i/g, eic. the residuals 
arising from subtracting tjie observed and adjusted values, and 
Xgv^ the sum oig^v^-\-g^^ + g^v^ + etc. Then the probable 
error of an observation whose weight is unity is 



(111) r= 0-6745 A /-^^^. 

and the probable error of an observation of the wiaight g^ is 
~=r (Art. 30). Also if G„ G^, etc. be the weights of the 

adjusted values of s, t, etc. their probable errors are 
(81) Jl..-^_. B,=^^,e... 

to determine which we must first find r, and the weights 
Gg, Ot, etc. If as in Arts. 44 and 45 all the observations 
are of equal, weight, we have oiily to make g, g^, g^, etc. 
unity or simply omit them from the formulae. To &id the 
weights Og, Ot, etc. we have only to form the independent 
observation equations by elimination of imknown quantities 
from the conditional equations (Art. 44), and then having 
multiplied each by the square root of its weight, proceed 
with the normal equations by the method of Art. 39. A 
single example will render the whole operation clear to the 
reader. 

We take a case of levelling involving but one conditional 
equation. There are three points -4, B, and C, situated at 
nearly equal distances apart, but upon diflFerent levels. In 
order to ascertain with accuracy their relative heights, a 
levelling instrument was set up between A^ and, 5, and 
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readings taken upon a rod held at those points, with the 
results, 

On Eod at -4, 8*7342 feet — ^mean of 12 readings, 

„ „ B, 2-3671 „ — „ 9 „ 

The instrument was then moved to a point between B 
and C, and the observations taken 

On Eod at B, 5*0247 feet — mean of 7 readings, 
„ „ (7, 11-2069 „ — „ 4 , 

Lastly, the level was set up between G and A, and the 
rods observed 

On Rod at G, 0*4672 feet — ^mean of 5 readings, 
„ A, 0-6510 „ — „ 3 „ 

It is required to find the adjusted values of these read- 
ings and their probable errors, also the most probable dif- 
ferences of level between the points and their probable 
errors. 

First let us atrange these measurements as they would 
be written in an engineer's level book, and assuming the 
elevation of J. as 0*0, find the heights of the other points. 



Sta. 


Back Sight. 


Fore Sight. 


Height of 
Instrument. 


Elevation 
above A, 




8-7342 
6-0247 
0-4672 


2-3671 

11-2069 

0-6510 


8-7342 

11-3918 

0-6521 


0-0 

6-3671 
0-1849 
0-0011 



The number of readings or the weight of each sight is placed 

in the first column preceding and following the name of the 

station ; thus ^B^ denotes that the back sight on B has a 

weight of 7, and the fore sight one of 9. Regarding the 

elevation of A as 0, that of B comes out 6*3671 feet, that 

of 0, 0*1849 feet, and on returning to the starting-point we 

find that A is 0*0011 feet instead of as it ought to be.T 
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If we represent the back sights upon Ay B, and O by 
Sy T, and W, and the fore sights upon B, G, and A by X, 
T, and Z, the rigorous conditional equation is 

Since this condition is not exactly fulfilled by the ob- 
served values, they must be corrected so as to cause the dis- 
crepancy of 0*0011 to disappear. To avoid the use of large 
numbers let us take the measured values as approximate, 
and represent by 8, tf w, etc. the connections to be applied to 
8y Ty Wy q\^ Thcu our conditional equation is 

«+« + w-a?-y-« == - 00011, 

The weights of «, ty w, x, y, and z are 12, 7, 6, 9, 4, and 
3, corresponding to the weights of the observations to which 
they are corrections. The most probable values of the un- 
known quantities which will exactly satisfy the conditional 
may then be found, either by the method of Art. 47 or by 
that of Art. 48. The latter is much the shorter and simpler. 
Leaving the solution as an exercise for the student, we give 
merely the results, viz. 



t(; =r - 0-00020, 
z = 000033. 



ff = - 0-00008, t = - 0-00014, 

»= 0-00011, y^ 0-00024, 

Applying these to the observed values we have the ad- 
justed results, viz. 



Sta. 


Back Sight 


Fore Sight. 


Elevation 
above A, 


A 
B 
C 
A 


8-73412 
6-02456 
0-46700 


2-36721 

11-20714 

0-65133 


0-0 

6-36691 
018433 
0-0 



To find the probable errors of these values we must first 
determine the weights of the adjusted back and fore sights 
by the method of Art. 39. We give the process in full for 
those necessary to determine the elevation of JB, that is for 
SajidX 
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Remembering that 8, t^ etc. are corrections to the ob- 
served values, the observation equations are 

5 = weight 12, a? = weight 9, 

< = „ 7, y = „ 4, 

t(? = „ 5, « = „ 3, 

and the conditional equation is as before 

« + ^ + w;-fl?-y-0 = - 0-0011. 

Taking from this the value of y in terms of the other un- 
known quantities, inserting it in the observation equations, 
and multiplying each by the square root of its weight (Art. 
36), we have 

^/^2J? = 0, 
s/7e = 0, 

3a? = 0, 

^/3« = 0, 

2* + 2* + 2t(; - 2aj - 2i8 «- 0-0022. 

From these we have the normal equations (Art. 34), 

16«+ 4e+4K;- 4aj - 4« = - 0-0044 « -4, 

4tf + ll^ + 4u?- 4aj -4« = - 0-0044 = jB, 

4«-f Uf^- 9m7 - 4aj - 4^ = - 00044 = C, 

-4«- 4«-4M>+13a?+4;5« 00044 = 2), 

-4d- 4«-4m;+ 4tx+7si^ 0*0044 = ^, 

the first being the normal equation for a, the second for t, 
the third for «(;, the fourth for x, and the fifth for z. 

To determine the weight of «, we place the absolute 
terms in the normal equations, equal to A, B, C, i), and E, 
and then solve the equations by any algebraic method. 
Then the weight of 8 will be the reciprocal of the coeffi- 
cient of -4 in the value of «, and the weight of x will be the 
reciprocal of the coefficient of i) in the value of x (Art. 39). 
The solution gives 
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3921 



8 = 



50832 



A + terms. in B, (7, A and E=- 0*00008, 



Q-IO 

Of = ^^D -f terms in A, B, G, and E = 0-00011. 
Hence the weight of the value of j? is 

and the weight of the value of aj is 

a = ?-^^ - = 9-99 
' 212 ^^^• 

The original weights of the observed values of 8 and X 
were 12 and 9. The adjustment has then increased each of 
these by nearly unity. (In finding the weight of 8, if we 
do not at the same time wish to find the value of «, the 
terms B, C, JD, and E may be made zero, and the numerical 
work thus shortened.) 

We can now find the probable errors. The residuals 
t?j, t?,, Vg, etc. are in this case our corrections 5, ty w, etc. : 
squaring these and multiplying each square by the weight of 
its corresponding observation, we obtain 

2^v'=.0'000001079. 

Since there are six observation equations, six unknown 
quantities and one conditional equation, we have n — 5+ p = 1. 
The probable error of an observation whose weight is unity, 
that is of a single reading of the rod, is then 



r = 0-6745 ^0-000001079 = 0-000701, 

and the probable error of any observation is this value 
divided by the square root of its weight. Hence the pro- 
bable error of the adjusted reading /S= 8*73412, whose 
weight we have found to be 1296, viz. 

r ^OjOOOTOl^^, 

JO, s/12-96 
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and the probable error of the adjusted reading X whose 
weight we have found to be 9*99, is 

5«=-7L: = 0-00022. 
^/lO 

The elevation of B is the difference of the readings S 
and X or 

B = 873412 - 2-36721 = 6-36691. 
Hence (Art. 41) its probable error is 



ii = Vi^/ + ie/ = 000029. 

The adjusted values may then be written 

Back sight on A^ 873412 ± O'OOOIO, 

Fore sight on 5 = 2*36721 ± 000022, 

Height of B above A = ^^OeQl ± 0-00029. 

In the same way the probable errors of the other adjusted 
measurements and the probable error of the elevation of G 
may be deduced. 

Problems* 1. In the preceding example what is the 
probable error of the observed readings b*7342 and 2*3671, 
and the probable error of the elevation 63671 1 

2.. The three angles of a triangle were measured as 
follows, all having the same weight, 

X = 98°. 17' 22", y ^ 70^ 9' 56", z « 11* 32' 52". 

What are the adjusted angles and their probable errors? 

Am. «= 98° 17' 18"-677 ± V 007. 
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CHAPTER V. 

THE DISCUSSION OF PHYSICAL OBSERVATIONS. 

50. In the preceding pages we have given the methods 
of adjust! ug and comparing such simple observations as 
arise in the common practice of the civil engineer. It is, 
indeed, in the delicate measurements of higher Qeodesy and 
Astronomy that these methods of combination find their 
most extended application, but it would interfere with the 
plan of this book to introduce examples of such cases. One 
or two applications, which we have not yet exemplified, 
are however of such great use, not only to engineers, but to 
all who cultivate the physical sciences, that we shall devote 
a few pages to their development and illustration. 

The Deduction of Empirical FormvlcB. 

51. Observations are frequently made to discover the 
laws which govern phenomena. It is one of the most 
happy applications of the method of least squares, that it is 
often able to determine firom such observations the relati(Mi& 
between the observed quantities, and to deduce formulsB 
expressing these relations in a convenient form for use. 
When such formulaa merely represent the particular obser- 
vations from which they are deduced, they are called empi- 
rical, but when they are applicable to all phenomena of the 
same class they become physical laws. The discussion of 
observations to deduce empirical relations thus often leads 
to the discovery of important laws by which our physical 
theories are extended and improved. 

To illustrate, let us suppose that the law of falling 
bodies is unknown, and that in order to discover the relation 
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between the time of falling and the space passed over, we 
construct an apparatus by which a body can be allowed to 
fall certain given distances^ and the times of its descent be 
accurately registered. Suppose that with this apparatus the 
following observations are made at New York, 

1. A body starts from rest and fall)3 10 feet in 0*788 seconds. 

2. „ „ „ 20 „ 1115 „ 

o 

**• » » » 

4 

^' » » >9 

*'• >> » » 

Now what relation or law exists between the space s and 
the time t ? First we observe that as the space increases so 
does the time, but the first very much faster than the second. 
The relation between a and t cannot therefore be expressed 
by such an equation aa s^^At where ^ is a constant, but a 
must depend upon higher powers of t In like manner we 
might try in succession the hypotheses that s=^At\8=At-{-Af, 
etc. Of course in this simple case the reader will at once 
recognize the relation which exists. In most cases, however, 
like those that will be hereafter given, the relation cannot be 
determined by inspection or even by the trial of hypotheses. 
In our ignorance, then, of whether $ depends upon t, f, or f, 
or upon a combination of these, it is best to write 

(112) 8^At + Bf+Ce, 

in which A, B, and G are constants to be determined, and as 
yet unknown. Then our first observation gives « = 10, 
t = 0*788, the second gives 8 = 20 and t = 1-115, etc. Substi- 
tuting then in the assumed formula, we have 

1. 10 = 0*788 A + 0*788« B + 0*788' (7, 

2. 2o=rii54-f rii5«5H-rii5»a 

etc. etc. 

There will be as many observations as there are obser- 
vations each containing the unknown constants A, B, and C. 
They are then the observation equations (Art,^32) from 
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ivhich we may derive a normal equation (Art. 34) for each 
of the unknown quantities, and by their solution find the 
values of J., B, and C. In the case before us this process 
will give 

-4 = 0, JB = 16-08, and (7 = 0, 

and hence the assumed formula becomes 

5=1608^. 

With reference to the five observations this is art empi- 
rical formula, but as it is found to represent all such obser- 
vations made at New York, it is hence also an expression of 
the law of falling bodies at that place. 

52. A very large class of physical phenomena may be 
represented by the general equation 

(112) y = -^+^«+C'^ + -D^ + etc. 

in which y and x are two related quantities of different kinds ' 
and Ay B, (7, etc. constants. As such we may indicate many 
of the relations between space and time, space and velocity, 
volume and temperature, pressure and density, etc. It in- 
cludes hence many of the empirical formulae of hydraulics, 
of the theory of heat, and of other branches of physics. No 
general rules, however, can be laid down to show when and 
where it is applicable. 

Temperature of the Earth at depths behw the surface*. 
It has long been known that as we descend below the surface 
of the earth the temperature increases. Geologists formerly 
supposed that the increase was nearly uniform, and at the 
rate of about 1 degree Centigrade to every 30 metres in 
depth. The observation at the artesian well at Grenelle near 
Paris, showed however that the increase was not at so rapid 
a rate. The observations there taken were the followingj: 

♦ From Von Frebden's Prctxis der Methode der Kletnsten Quadrate. 
Braunschweig, 1863. 

t If F=Fahrenheit degrees and C = Centigrade degrees, F=82+ JO. 

1 meter=3-28087 EnglUh feet=3-28071 U. S. feet. t 

o 



66 


if 


= 12 -90 


173 


» 


= 16 -40 


248 


»> 


= 20 00 


298 


»> 


= 22 -20 


400 


» 


= 23 -75 


503 


» 


= 26 -46 


548 


99 


= 27 -70 
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Mean yearly temperature of the surface = 10*'60 Centigrade. 

1. Temperature at a depth of 28 meters = 11 -71 
o 

^' ft • M >> 

**• » 99 » 

^' »> 99 99 

*^' it » » 

^* » >» >» 

7 

"• » >» » 

Let us represent the mean temperature of the surface by 
T, and the temperature at any depth by t : also let x be any 
depth. We then assume that the temperature increases with 
the depth according to the approximate law, 

t^T + Ax + Ba^, 

and we have to find A and B from the above observation. 
The first one gives t = ll"-7l and x ~ 28; or 

ll-71 = 10-6 + 28ul+28»5. 

The second one gives 

12-90 = 10-6 + 66 ^ + 66«jff. 

Transposing the constant term 10°*6 to the first members, 
the observation equations may be written, 

ril= 28^+ 28»J?, 

2-30= 66^+ 66' J?, 

5 -80 = 173 -4 + 173' JS, 

9 -40 =248 ^ + 248' 5, 

ll-60 = 298^ + 298«JS, 

13 -15 = 400 ^ + 400' 5, 

15 -83 = 505 ^ + 505" 5, 

17 -10 = 548 ^ + 548' 5. 

We now proceed to form the normal equations ; multi- 
plying the first equation by 28, the second by 66, the third by 
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173, etc., and adding tke products, we have the normal equa- 
tion for A, viz. 

29599-23 = 900706 A + 404557966 B. 

Also multiplying the first equation by 28*, the second by 
66', etc., and adding the products, >we find 

13068985-39 = 404557966 A + 1934100018145 

as the normal equation for B. Combining these two normal 
equations we have the values 

A = 0042096, jff = - 0000020558, 

and hence the empirical formula is 

t = 10*-6 + 0'-042096 x - 0^-000020558 x'. 

. Thus by this formula the temperature t at any depth w 
may be found, x being in meters and t in Centigrade degrees. 
The negative sign of the term containing x* shows that the 
increase of temperature is not so rapid as the depth. If in 
this formula we substitute a? = 28, a; = 66, etc., and compute 
the corresponding values of t, the accuracy of the formula 
may be seen by comparing the computed and the observed 
temperatures. The following table is such a comparison. 



Depth. 


Observed 


Calculated 


temperature. 


temperature. 


Cm. 


10-60 C. 


10-60 C. 


28 


11'71 


11-76 


66 


12-90 


13-29 


173 


16-40 


17-26 


248 


' 2000 


19-78 


298 


22-20 


21-32 


400 


23-76 


24-16 


505 


26-43 


26-61 


548 


27-70 


27-49 


1000 




3213 



The close agreement between experience and calculation 
show that as far as the above observations are concerned, the 
formula is sufl&ciently accurate. A more recent boring, at 
Speremberg near Berlin, which reached a depth of about 
1200 meters, seems to indicate a somewhat faster rate of 
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increase. At a depth of 1000 metres the observed tempe- 
rature was about 43** C, while that computed by our formula 
is only 32^*13, and it can hence be regarded only as empi- 
rical, and as representing the Grenelle observations. A dis- 
cussion of the results from this deep boring at Speremberg 
may be seen in Nature, Oct. 21, 1875, where it is claimed 
that after a depth of about 5000 feet the temperature is 
uniform, and probably not more than 49° Centigrade or 120° 
of Fahrenheit's scale. 

Problem. Volume of water, at different temperatures. In 
MuLLEB's Lehrbiich der Physik, Vol. II. p. 505, the following 
are the relative volumes of water at temperatures from 4° C. 
to 23° as derived from experiments by Despretz : 

At 4° Ceutigrade the volume is 1* 

6 „ „ 10000309, 

8 



10 
15 
20 
25 



1-0001216, 
1-0002684, 
1-0008751, 
1-0017900, 
10029300. 



Required a formula to represent these experiments. 

Let F represent the volume at any temperature t Then 
by the same process as that of the preceding example, we 
find 

7= 1 - 0000061045 1 + 00000077183 <•- 0-00000003734 f. 

Although agreeing tolerably well with the observations 
this formula does not hold for temperatures much higher than 
25° C. From 25° to 50° Mctller gives 

F= 1 - 0-000065415 1 + 00000077587 ^-0-00000003541 f. 

It appears then that the law connecting the volume and 
temperature of water is not yet understood, and hence that 
the above formulae are merely empirical 

63. Another large class of phenomena may be repre- 
sented by the general equation 
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/iijx >! . n • 360 . „ 360 
(114) y=^A + B. sin a?+jB,cos x 

+ C, sm 2a? + C, cos 2a?+ek 

in which as a; increases y passes thotigh repeating cycles. 
As such may be mentioned the variation of temperature 
throughout the year, the changes of the barometer, the ebb 
and flow of the tides, the distribution of heat on the surface 
of the earth depending on latitude, aiid in fact all pheno- 
mena which repeat themselves like the oscillations of a pen- 
dulum. The letters A, B^, B^, etc. represent constants wnicb 
are to be found from the observations by the method of least 
squares, while m is the number of equal parts into which the 
whole cycle is divided, and must be taken in terms of the 
same unit as x. If the several cycles are similar and regular, 
only the first three terms are required to represent the 
variation, if not, additional terms are necessary. 

Declination of the magnetic needle. The magnetic needle 
has always in the New England States since the earliest 
records pointed to the west of true north. Its declination 
has moreover been slowly changing. Numerous observations 
from all parts of the country have been collected and care- 
fully discussed by Mr Schott of the U. S. Coast Survey*, 
from which it has been established that the declination 
passes through a recurring cycle of about 250 or 300 years. 
If we divide a horizontal line into equal parts representing 
years, and erect verticals giving the declination at various 
times, we shall have a curve, such as is roughly represented 

Fig. 14. 




i900 



noo 



JBOO 



J900 



in the annexed figure. Thus in Connecticut the westerly 
declination Vas in the year 1600 gradually increasing; it 

* See U. 8. Coast Survey Reports, 1865, 1858 micL 1859^ 
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attained a maximum of 11° or 12° about 1675, then again 
decreased until it reached a minimum of 5° about 1800, and 
has since been steadily but not uniformly increasing. 

If then we have observations enough to determine the 
equation of this curve, we shall not only have an empirical 
formula representing those observations, but an expression of 
the law of the change. At any place where observations 
extending over fifty years exist, an approximate formula may 
be found. 

At Hartford, Conn., we have the following records : 

In 1786 the declination was 5°25' W= + o°-42, 

„ 1810 „ „ 4 46 „ = 4 -77, 

„ 1824 „ „ 5 45 „ = 5 75, 

„ 1828 „ „ 6 3 „ = 6 05, 

„ 1829 „ „ 6 3 „ = 6 -05, 

„ 1859 „ „ 8 4 „ = 8 -07. 

We assume the expression 

n ^ . D • 360\, ^ 360\ 
D=^A-{'B,sm t + B^ cos t, 

as a sufficient approximation to the law of the declination, 
in which D is the value of the declination at any time t, t 
being the number of years counted from an assumed epoch, 
say the year 1830, m the number of years in the cycle, and 
-4, jBj and B^ constants to be determined. The value of m 
is not exactly known, and must be adopted from observations 
at other places, or be determined by trial. Mr ScHOTT finds 
its most probable value to be 288 years for Hartford. Then 
our equation is 

D^A-^-B^ sin 1-25 t-^B^ cos 1*25 t 

The first observation is 2) = 5'**42 for i = — 44, the second 
D = 4^77 for « = ^ 20, the last D = 8°-07 for « = + 29-6, 
t being the number of years counted from 1830 as an 
origin, and hence negative for the years preceding that epoch 
and positive for those following. The first observation equa- 
tion is then 

5H2 = A- 0-819 B, + 0-574 B^^^^ .^Google 
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and the five others are similar. (The sign of the sine and 
cosine of a negative arc should be carefully regarded.) 

From the six observation equations we form and solv6 
the normal equations, and find 

A = 8«-60, B, = 2-54, 5, « - 2-54, 
hence our formula is 

i) = 8'-60+ 2^-54 sin r-25^-2«-54 cos r-25<, 

which agrees very closely with the observations. This equa- 
tion may be discussed like that of any algebraic curve, and 
the times and values of the maximum and minimum decli- 
nation found. Thus, according to the formula, the last mini- 
mum at Hartford was 5®'01 for ^ = — 34, or in the year 1796 ; 
the next maximum will be 12'*19 for < = 110, or in the year 
1940. The formula therefore extends the law of the varia- 
tion with a fair degree of accuracy to times considerably 
preceding and following the observations. As indicated in 
Part II. it may be placed under the forms 

D = 8^-60 + 3°-59 sin (l'*-25 1 - 45°), 

D = 8°-60 - 3^-59 cos (l°-25 1 + 45°), 

either of which is more convenient for discussion than that 
given above ; and the latter of which is the form given by 

SCHOTT*. 

54. Other general empirical and theoretical fonmilae 
such as 

y = -4 + Bx^ + Cb* + etc. 

y = ^ + 5 sin mcc + C sin* mx + etc., 

are sometimes used in discussing physical phenomena. Rarely 
also coeflBcients such as A^, ff, etc. are introduced, but their 
determination is generally laborious (see Part II. Art. 59), 
Exactly what formula will apply to a given set of obser- 
vations, so as to agree well with them, and at the same time 
be of use in other similar cases, can only be determined from 
theoretical considerations or by trial. The computer mu8t> 
first, from his knowledge of physical laws, assume such au 

* U. 8, Coast Survey Report, 1859, p. 298, or Am, Journal of Science^ 
1860, p. 335, where decimation fonnrilsd for many other places are given. 
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expression as seems most plausible, then deduce the con- 
stants by the method of least squares, and test them by 
comparing the observed with the calculated results. If the 
agreement is not sufficiently close, he must assume another 
expression and try again. In this way hypotheses may be 
tested, and often important physical laws be discovered. 

The precision of such empirical formulae may be ascer- 
tained not only by their agreement with the observations, 
but by computing the probable errors of the constants, 
which the methods of least squares has deduced (Art. 38). 
The probable errors of the calculated results may also be 
found!, but for such computations we must refer the reader 
to the larger and more complete treatises upon the subject, 
a list of which will be found m the Appendix. 

Problem. The following observations have been made 
in different parts of the earth, upon the length of the 
second's pendulum. 



Plact. 


LaUtode. 


Length of Pendulum 
in meters. 


Peru 

Petit Goave . 

Paris . 

St Petersburg. 

Lapland 


00 
180 27' 
48 24 
6S 16 
67 4 


0-990664 
0-991150 
0'993867 
0-994589 
0-995326 



Find a formula representing these observations, observing 
that according to the theory of the variation of gravity it 
must be of the form 

? = ^+£sin*^, 

4 

in which <^ is any latitude, and I the corresponding length 
of the second's pendulum. 

Ans. I = 0"*-990555 + 0~-005679 sin'^. 

BowDiTCH, in hia Commentary on Laplace's Mica- 
nique Celeste, finds a formula for the length of the pen- 
dulum by using fifty-two observations. The student who 
has solved the above problem will have some idea of the 
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patience and labour required in such processes, and will be 
interested in looking over Bowditch's admirable discussion, 
which is in Vol. ii. p. 481 of the above-cited work. The 
formula which he adopts as the best is 

I = 3901307 + 0-20644 sin* ^, 

in which the constants and I are in inches, while those of our 
problem are in meters*. Ivory's interesting papers in the 
Jjondon Philosophical Magazine for 1826 may also be referred 
to in this connection. 

Prdbahility of Errors, 

55. In Arts. 12 and 13 we have stated a property of the 
probability curve, and given a table of the probabilities of 
errors, which has many interesting applications in various 
branches of science. That table gives the areas of the 
curve on both sides of the axis of F, corresponding to suc- 
cessive numerical values of hx, h being the constant measure 
of precision and x any error, or, more strictly, it gives the 
number of eiTors comprised between the limits — x and + a:, 
when the whole number of errors is unity. As the fre- 
quency of an error is proportional to its probability (Art. 5), 
the numbers in the columns P' denote also the probabilities 
that an error will be comprised between the assumed limits. 
Thus if in the figure representing the probability curve, 
we lay off the positive error OM equal to the negative error 
OMj and draw the ordinates MG and if (7, the area, MCACM 
is (if the total area be unity) a fraction expressing the 
probability that an error will lie between the limits— OM and 
+ OM, or be numerically less than OM. So also if the area. 
PBAPB be 0-5, OP is the probable error (Art. 16), or the 
error such that the probability of an error less or greater 
than it is \, If we were to lay off an abscissa positively 
and negatively equal to three times OP, the corresponding area 
would be 957, and the probability that an error taken at 

• 1 meter= 89-37044 English inclie8= 39 '36862 CT.Sf. inches. I do not 
know whether the above formula is in English or United States inohM. 
This annoying difference between onr own measures and those of England 
is one of the strongest reasons why we should throw both of them oyerboard 
and adopt a reasonable system. 
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Fig. 1. 



Ill 




g5'7 

random would be less than 3. OP is tkkf\ > and that it would 



1000' 



43 



be greater ^7^. ; so that it would be a wager of 957 to 43, 

or of nearly 22 to 1, that any error is less than three times 
the probable error* 

In discussing the probabilities of errors and the accuracy 
of sets of observations, it is hence convenient to have a 
table giving the numbers P' in terms of r instead of L 
The following is such a table, which we have taken from 
Chatjvenet's Treatise on Least Squares^ and which gives 
these areas or probabilities, P corresponding to successive 

numerical values of - , a? being any error and r the probable 

error. To illustrate its use let us consider some practical 
examples. 

1. Suppose that we have measured an angle with two 
instruments A and B, and find for the results 

with A 37" 42' 13"-9 ± O'H, 

with B 37' 42' 18"-8 ± 0"-2, 

the probable error of the first being twice that of the second, 
and hence its precision one-half as great (see Fig. 2, Art 16). 
The probability that the result A does not differ Jirom the 

y Google 
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Prohahility of Errors, (116). F 






e-^dt, t=fer=:0-4769- 



X 

r 


F 


X 

r 


F 


X 


F 


X ■ 

r 


P 


0-00 


0-00000 


1-20 


0-68171 


2-40 


0-89460 


3-60 


0-98482 


0-04 


•02152 


1-24 


•59705 


2-44 


•90019 


8-64 


-98688 


0-08 


•03228 


1-28 


•61206 


2-48 


•90562 


8-68 


•98691 


012 


•06451 


1-32 


•62679 


2-62 


•91082 


3-72 


•98786 


0-16 


•08594 


1-86 


•64102 


2-66 


•91678 


8-76 


•98872 


0-20 


0-10731 


1-40 


0-66498 


2-60 


0-92061 


8-80 


0-98962 


0-24 


•12860 


1-44 


•66858 


2-64 


•92608 






0*28 


•14980 


1-48 


•68184 


2-68 


•92934 






0*82 


•17088 


1-62 


•69474 


2-72 


•93344 






0-86 


•19186 


1-66 


•70729 


2-76 


•93784 






0-40 


0-21268 


1-60 


0-71949 


2-80 


0-94106 


4-00 


0-99302 


0-44 


•23336 


1-64 


•78134 


2-84 


•94468 






0-48 


•25388 


1-68 


•74286 


2-88 


•94793 






0-62 


•27421 


1-72 


•75400 


2-92 


•96111 






0-66 


•29436 


1-76 


•76481 


2-96 


•96412 






0-60 


0^31430 


1-80 


077628 


8-00 


0-95698 


6-00 


0^99926 


0-64 


•83402 


1-84 


•78642 


8-04 


•96968 






0-68 


•35362 


1-88 


•79622 


3-08 


•96224 






0-72 


•37277 


1-92 


•80469 


812 


•96466 






0-76 


•89178 


1-96 


•81883 


8-16 


•96694 






0-80 


0^41062 


2-00 


0^82266 


8-20 


0-96910 






0*84 


•42899 


2-04 


•83117 


824 


•97114 






0-88 


•44719 


208 


•83936 


8-28 


•97306 






0*92 


•46509 


2-12 


•84726 


8-82 


•97486 






0-^6 


•48270 


2-16 


•85486 


8-86 


•97573 






100 


0-60000 


2-20 


0-86216 


8^40 


0-97817 






1*04 


•61699 


2-24 


•86917 


844 


•97967 






1-08 


•63366 


2-28 


•87691 


8-48 


•98109 






112 


•66001 


2-32 


•88237 


862 


•98237 






1-16 


•66602 


2-86 


•88867 


8-66 


•98360 


00 


1^00000 
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true value of the angle by 0"*4 is ^, and the probability 

that the result JB does not diflFer by 0"'2 is also ^ . Now 

what are the respective probabilities that these results are 
correct within 0"*1 ? We take then x = 0"'l, and have 

for A, -=:,^ = 0-25, and for 5,^=^ =0-5; .. 

r ()4 r 0*2 ' v^ 

. hence from the table 

for - = 0-25, P' == 0-134., and for - = 05, P' = 0264. 
r r 

The probability that the result A is within 0"*1 of the truth 

is ■ , and that B is within the same limit txt^v or nearly 

twice as much. Hence for A we could wager 134 to 866 or 

1 to 5 that the result was within 0"*1 of the truth, and for 
B 264 to 736 or 1 to 2*8 that such was the case. Again, 
what are the respective probabilities that these results are 
within 0"-6 of the truth? We have 

« . a? 0-6 - ^ , « ^#0? 0*6 ^ ^ 
for -d, - = p-r = 1*0, and for ^, - = tt-tt = 3*0, 
r 0*4 ' r 0-2 

and then from the table 

for ^ =; 1-5, P' = 0-688, and for - = 3, P' = 0-957. 
r T 

The probabilities are then ^kf^ and j^ respectively ; 

and we could afford to lay a wager of 688 to 312 or of about 

2 to 1 that A is within 0"'6 of the truth, and one of 957 to 
43 or of more than 22 to 1 that B is within the same limit. 

2. A series of 64 observations upon an angle gives for 
the probable error of a single observation 1"*5. What is the 
probability that the mean is correct within 0"'76? 

The probable error of the mean will be -j== = 0"'185. 
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Then - = ^^ = 4-05, for which P^ 0-994 

r 0185 

994 
The required probability is hence ^t^f^ , and it is a wager of 

994 to 6 or of 166 to 1 that the mean is within 0"'75 of the 
true value. 

3. An angle is measured with an instrument graduated 
to 1'. The error which is liable to occur in a single measure- 
ment (that is, the probable error r) is 45". How many obser- 
vations are necessary in order that it shall be a wager of 9 to 

1 that the mean is within 5" of the truth ? 

9 
A wager of 9 to 1 corresponds to a probability of ^ . For 

P' =: 0'9 we find from the table - = 2*44, in which r, being the 

r 

45" 
probable error of the average, is equal to -j= . Then we 

have 

2-44 = - = ^ = ^, or ^ = 21-96, and n = 482, 
r 45 45 ^ 

and hence 482 observations are necessary. 

Problems, 4. The average of several observations gives 
for the value of an angle 33' IT 30"-8 ± 0"-3. What wager 
can we afiford to lay that 33" 17' 30"-8 is within 1" of the true 
value? Ans. 39 to 1. 

6. A line is measured 500 times. If the probable error 
of each observation is 0*6 centimeters, how many errors will 
be less than 1 centimeter and greater than 0*4 centimeters? 

6. An angle is measured by an instrument graduated to 
quarter degrees, the probable error of a single reading being 
12 minutes: how many observations are necessary that it may 
be a wager of 5 to 1 that the mean is within one minute of 
the truth? Ans. 22 observations. 

7. The length of a line was estimated by 16 persons of 
equal skill. If the probable error of each is 4 centimeters, 
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what is the probability that the mean of their guesses is 
within 2 centimeters of the truth? 

56. The preceding examples illustrate the great value of 
large numbers of observations even when made with poor 
instruments, provided only that no constant cause of error 
exists (Art. 2). The preceding table is thus of the greatest 
assistance in discussing statistics of social science, even when 
those statistics are confessedly inaccurate, provided only that 
they are numerous and taken by unprejudiced observers. It 
is however not the place here to give examples of such appli- 
cations, and we would refer the reader desirous of looking up 
such investigations to the admirable work of Quetelet (No. 
17 in our list of literature in Art. 64 of the Appendix), in 
which numerous interesting examples of such discussions may 
be seen. 

The principles of probability of error here set forth are, it 
must be borne in mind, only applicable after all constant 
errors (Art. 2), and all mistakes (Art. 3) have been elimi- 
nated from the numerical results. If a single cause of con- 
stant error exists it may sometimes be detected by a com- 
parison of the results with those obtained by a more accurate 
instrument. The following illustration may in this connection 
be interesting. 

Suppose that an angle is laid out with very accurate in- 
struments and tested in many ways so that its true value may 
be regarded as exactly 90®. Let 25 observations be taken 
upon it with a transit whose accuracy we wish to test, and 
let the mean of those measurements be 89° 59' 57" + 0"'8. 
Then we see it is extremely probable that a constant error of 
about — 3" exists in the instrument. To find the numerical 
expression of this probability let us suppose that the time 
value of the angle was unknown, and let us ask the proba- 
bility that the mean is within 2" of the truth. Then for 

ic 2 

- = ^7^ = 2-5 we find F = 0'908, so that it is a wager of 908 

to 92 or of almost 10 to 1 that the mean is between the 
limits 89* 59' 55" and 89° 59' 59". Hence since the angle is 
known to be 90°, it must be the same probability and the 
same wager that there is a constant error lying between the 

8-§ 
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limits — 1" and — 5". So also if we take ^ = 3" we can show 
that it is a wager of 39 to 1 that there is a constant error 
between 0" and - 6". 

Problem, The capacity of a certain large vessel is un- 
known: 1600 persons guess at the number of gallons of water 
which it will hold and the average of their guesses is 289 
gallons. The vessel was then measured by a committee and 
found to hold 297 gallons. If we regard the probable error 
of a single guess as 50 gallons and also consider it impossible 
that there can be any constant source of error in guessing, 
what is the probability that the committee made an error in 
their measurement of between 3 and 13 gallons? 

993 
Am. frjTxT)* or a wager of 142 to 1, 



The JRejection of Doubtful Observations. 

67. It not unfrequently happens that in a set of measure- 
ments there are certain values which seem to be so much at 
variance with the majority, that the observer rejects them 
in adjusting the results. This procedure however, unless 
governed by proper rules, is a dangerous one and not to be re- 
commended. A conscientious observer having conscientiously 
made several series of measurements will give each its proper 
weight (Art. 28) and deduce the most probable result, and 
give to it the confidence which its probable error shows that 
it deserves. The too common practice of taking twenty 
measurements (for instance, readings of a levelling rod) and 
then throwing away all except two or three which happen to 
agree, is one which cannot be strongly enough condemned. 
After making such measurements and eliminating all known 
constant errors (Art. 2), no results except those which are 
unquestionable mistakes (Art. 3) should be rejected. All 
remaining discrepancies will then fall under the class of 
irregular errors (Art. 3), and the adjustment of such observa- 
tions should be made in accordance with the principles 
governing them, the methods for which we have presented in 
the preceding Chapters. ^^^^^^ ^yGoogle 
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We must mention, however, that for very delicate and 
precise observations such as sometimes arise in Astronomy, 
the principles of probability itself furnish a means of deter- 
mining whether or not a given observation may be rejected. 
The discovery of this important criterion is due to Prof. 
Peirce of Harvard University, and presented in the (7am- 
hridge Astronomical Journal for 1856. The student wishing 
to inform himself concerning its theory and application, may 
consult the original paper in that Journal or the works Nos. 
28 and 34 quoted in the list of literature given in Art. 64 of 
our Appendix. The use of this criterion is not necessary ex- 
cept in some very accurate astronomical investigations. 



Concluding Remarks, 

58. The student who has carefully read the foregoing 
pages and has verified the examples and solved the problems 
presented, will have acquired a fair knowledge of the princi- 
ple of Least Squares and of its simpler applications to engineer- 
ing practice, and will be prepared to study the more complete 
theory of the subject with interest and profit. The second 
part of this book will afford him an introduction to the theory, 
which has been prepared with especial reference to the needs 
of beginning students. He will also be prepared to take up 
the adjustment of Geodetic and Astronomical observations 
involving large numbers of observations, and as excellent 
books for study and consultation in this connection, we would 
refer him particularly to the works of Chaxjvenet, Von 
Freeden, Dienger and Helmert, whose titles will be found 
in the list of literature appended in Art. 64. The brief his- 
tory of the origin and development of the science given 
in the Appendix will also prove suggestive in directing his 
further studies. 

We take occasion to again call attention to the fact that 
the formulae given for probable errors are all based upon the 
supposition that the number of measurements is large. 
Hence in using different formulae a perfect agreement in the 
results is not to be expected, unless suflSicient observations 
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have been taken to exhibit the several errors in proportion to 
their respective probabilities, and this would require a very 
large number. We would also again mention that the whole 
theory is based upon the supposition that only accidental or 
irregular errors (Art. 4) affect the measurements. 

Although we have only given methods for the adjustment 
of observations involving equations of the first degree they 
are sufficient for the discussion of all cases ; for, by the pro- 
cess of Art. 59 (Appendix), equations of higher degrees can 
always be reduced to linear equations, and the observations 
thus be brought under the rules which we have presented. 
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PART 11. 

THE THEORT OF LEAST SQUARES AND PROBABLE 
ERRORS, 

CHAPTER VI. 

DEDUCTION OF THE FUNDAMENTAL PRINCIPLES. 

1. In Part I. we have presented the rules and methods 
for the adjustment and comparison of ordinary observations 
and have illustrated their application by numerous practical 
examples and problems. We shall now take up the subject 
from a more mathematical point of view, and give demon- 
strations of all the formulae which have been there employed 
without proof, and also discuss the met^hod in a more com- 
plete and algebraic manner. While the previous chapters have 
been written with more especial reference to the wants of 
practical computers, the following pages will be designed 
to meet the requirements of students who wish to acquire 
a tolerably thorough knowledge of the theory qf the subject. 
Although each Part may be read independently, each is 
in fact a supplement to the other, and hence in order to 
render reference easy, we use the same numbers to mark the 
corresponding paragraphs and formulae. 

Since all measurements even when made with the utmost 
precision give discordant results, all of which cannot be 
true, it is evident that we can never be sure that we have 
found the absolutely true value of a quantity, which has 
been the object of measurement. In combining then such 
measures or observations, we seek a method which shall 
furnish us with the most advantageous or most probable 
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result, that is, a result which (as far as our observations go) 
we can regard as the nearest approximation to the true 
value. Such is the method of Least Squares, the theory 
of which we are to develope. Further, as we cannot regard 
our adjusted result as absolutely true, we must also establish 
what degree of confidence it is entitled to, and this involves 
the theory of probable errors. The first treats then of the 
adjustment, the second of the comparison of observations. 

When we measure several times a quantity and obtain 
discordant results, we recognise that each measurement is 
probably incorrect. The difference between the true and 
an observed value is called an error; it is taken as positive 
if the true value exceeds, and negative if it is less than 
the observed one. Since the true value cannot be exactly- 
determined, these errors can never be definitely known, 
nevertheless they can be made the subject of mathematical 
investigation, 

2. Constant errors which always affect our observations 
by the same amount and whose causes are understood, are 
no longer errors, as they may be always eliminated from 
our numerical results. They are not the errors which we 
are to discuss. 

3. Mistakes, whose distinguishing feature is that they 
"admit of conjectural correction," are also not included 
among the errors of which we are to treat. 

4. Accidental or irregular errors are those discrepancies 
which remain after all constant errors and mistakes have 
been eliminated, and are hence those produced by irregular 
and varying causes whose degree or manner of action cannot 
be estimated. The word error, as used in this book, means 
therefore "discordance, of which the cause is unknown." 
These errors being produced by many causes, all unknown 
as to their laws, are governed by the principles of pro- 
bability. 

Probability. 

5. We give therefore, by way of introduction, the de- 
finition and some of the first principles of probability. 

o 
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The word prohdbility as used in mathematical reasoning 
means a number less than unity, which is the ratio of the 
number of ways in which an event may happen or fail to 
the total number of possible ways. Thus if an event may 
happen in a ways and fail in h ways, and each of these ways 
is equally likely to occur, the probability of its happening is 

7 , and the probability of its failing is v. Thus 

probability is always expressed by an abstract fraction, 
and is a numerical measure of the degree of confidence 
which we have in the happening or failing of an event. 
If the fraction is it denotes impossibility, if \ it denotes 
that the chances are equal for and against its happening, 
and if it is 1 the event is certain to happen. 

6. Hence unity is the mathematical symbol for certainty. 
And since an event must either happen or not happen, the 
sum of the probabilities of happening and failing is unity. 
Thus if P be the probability that an event will happen, 
1 — P is the probability of its failing. 

7. If an event may happen in a ways and also in a 
ways and fail in b ways, the probability of its happening 

Cb "t* Cb 

is, by Art. 5. T-r\ \ and since this is the sum of the 

probability of happening in a ways and of that of happening 
m a' wsijs, it follows that if an event may happen in different 
independent ways the probability of its happening is the 
sum of the separate probabilities. 

8. Let us now ask the probability of the concurrence 
of two independent events. Let the first be able to happen 
in ttj ways and fail in 6 ways, and the second happen in 
a, and fail in 6, ways. Then there are for the first event 
®i "+■ K possible cases, and for the second a^ + b^ : and each 
case out of the a^ + b^ cases may be associated with each 
case out of the a^ + b^ cases, and hence there are for the two 
events {a. + b^) (a, + oj compound cases each of which is 
equally iitely to occur. In a^a^ of these cases both events 
happen, in bfi^^ both fail, in afi^ the first happens and the 
second fails, and in a A the first fails and the second 
happens. Thus we have for two independent events 
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Probability that both happen 



Probability that l)oth fail 



JA 



Probability that the first happens and the second fails 

Probability tliat the fii;st f^.and the second happens 

and the sum of these is unity since one of the four events 
is certain to occur. Now considering each event alone the 

probability of the first ^happening is ^, and of the 

second ^ , and since 

a, + b^ 



we have established the important principle, that the pro- 
bability of the concurrence of several independent eveuts, 
is equal to the product of the separate probabilities. 

Thus if there be four events and P^, P^, P^ and P^ be 
the respective probabilities of happening, the probability 
that all the events will happen is F^F^P P^y and the pro- 
bability that all wiU fail is (1 - PJ (1 - PJ (1 - PJ (1 - PJ. 
The probability that the first happens and the other three 
failisP,(l-P5^(l-P3)(l-Pj; and so on. 

9. If there be several events whose separate proba- 
bilities are known, what is the case most likely to occur? 
Let P be the probability of the happening of an event in. 
one trial and Q the probability of its failing so that 
P+ Q = l: and let there be n such events. Then by 
the preceding Art. the probability that all will happen is 
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P*; the probability that one assigned event will fail and 
n — 1 happen is P^'^Q, and since this may occur in n ways 
the probability that one will fail and w — 1 happen is nP'''^Q. 
Sinularly the probability of two assigned events failing 
and w — 2 happening is P""* Q, and since this may be done 

in —^-5 ways*, the probability that 2 out of the whole 

number will fiiil and w-2 happen is ^ — ^P""*Q*. 

If then (P + Q)" be expanded by the binomial formula, 
thus, 

(P + (2r = P* + nP«-^Q + '^J^P*^Q« + ... 

the first term is the probability that all will happen, the 
second that n — 1 will happen and 1 fail, and the m + 1^ 
term is the probability that n — m will happen and m 
fail. 

To determine then the most probable case we have only 
to find the term in this series which is the greatest. If we 
consider that n coins be thrown, P= Q = \^ and the series 
becomes 

«(«-!)(«- 2). ..(w-m+l) . , 
"*■ 1.2.3...m W +••• 

in which if n is even the middle term is the greatest, and if 
n is odd there are two equal middle terms greater than any 
other. Thus if » = 9, the series is 

1 9 36 84 126 
612 "^ 612 "*■ 612 ■*" 512 "^ 512 

126 84 36 9 1 
"^ 512 "^ 512 "^ 612 "•" 512 "•" 612' 

* See ToDHimTEB's Algebra for Schools and Colleget. London, 1870, 
p. 288, p. 466. i,i,„ed by Google 
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and lience if 9 coins be thrown — ^ is ^^ probability that 

9 

all will be heads, ^^o ^^^ probability that 8 will be heads and 

1 tail, and so on. Since one of these 10 groups must happen 
the sum of the series is unity. The most probable group is 
that whose probability is the greatest, and this is the one 

corresponding to ^y^, viz. 5 heads and 4 tails, or 4 heads and 

tails. 

If then we have several events either simple or com- 
pound we recognise that the most probable case is that whose 
mathematical probability is the greatest Hence if in given 
sets of observations various adjusted values of the measured 
quantity have different probabilities, the one to be chosen is 
that which has the maximum probability (Arts, 11, 15), 



Law of the Probability of Error. 

10. Although it would seem at first sight that acci- 
dental errors of observation could hardly be made the subject 
of mathematical reasoning, yet the very fact of their irregu- 
larity brings them under the laws of probability. Moreover 
we recognise that they must be subject to the following 
fundamental laws of arrangement: 1st, Small errors are 
more frequent than large ones ; 2nd, Positive and negative 
errors (that is, measurements greater and less than the true 
value) are equally probable, and hence in a large number of 
observations are equally frequent; 3rd, Very large errors do 
not occur, so that in every set of observations there is a limit 
I, such that all the positive errors are included between and 
+ 1 and all the negative ones between and — l. 

These are the three fundamental axioms which, in connec- 
tion with the principles of probability, are the foundation of 
aU our following reasoning. 

11. Therefore, different errors are not equally probable ; 
for in a large number of observations a small error occurs 
more frequently than a large one, and hence has a greater 
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probability, and an error greater than the limit I has a pro- 
bability of or is impossible. Hence the probability of an 
error is a function of that error, so that calling x any error 
and y its probability, the law of probability of error is repre- 
sented by the equation 

(1) y=/{^l 

and will be determined if we can find the form of /(a?). , 

Although practically there is a limit in the graduation 
and use of instruments by which x can have only definite 
numerical values (thus if an observer reads a theodolite to 
10", the values of x can only differ by 10" or some multiple 
of 10"), we must in our mathematical treatment regard x as 
a continuous variable. This is evidently perfectly allowable, 
since as the precision of our observations increases the suc- 
cessive values of x are separated by smaller and smaller 
intervals. Taking y then as a continuous function of x the 
equation (1) represents a curve whose form we are to deter- 
mine. The axioms of Art. 10 show that its general form 




must be that of Fig. 1, for by the first small values of x 
must have the largest probabilities y, by the second f{x) 
must equal /(— x) or the curve be symmetrical to the axis of 
Y, and by the third y must be zero for all values of x 
greater than ± I. This last requirement is one which ex- 
tremely embarrasses the mathematical treatment of the sub- 
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jeot, since it is impossible to determine a continuous function 
of X which shall become zero tor x = ±l and also be zero for 
all values of x from ±lto ±00. But since this limit I can 
never be accurately assigned we shall extend our limits to 
± 00 , and shall determine our curve in such a way that the 
value of y, although not zero for large values of x, will be so 
very small as to be practically inappreciable. 

For determining the equation of this curve we give two 
demonstrations, the first for the case of direct observations 
(Art 18), and the second for the more general case of indirect 
ones (Art. 19), both of which lead to the same result. These 
demonstrations are due to Gauss (see Art. 65). 

11a. Let one and the same quantity, for example the 
angle M, be measured n times with equal care, giving the 
values M^y M^y -3f« ... M^. If z is the true value of the angle 
the errors {z—Mj), (^ — JbTJ, .... (z — MJ have been com- 
mitted, and their respective probabilities are 

(1) y,=f{z-M,), y,=f{z-M:) ... y.=/(^-ig, 

that is, if 71 is the whole number of errors, and the error 

{z — Ifj) occurs Wj times, y^=^— —f{z — M^ ; and we suppose 

enough observations to be taken to exhibit the several errors 
in proportion to their respective probabihties. Now the pro- 
bability of committing all these errors is by Art. 8 the 
product of the separate probabilities, or if P denote this 
quantity , 

(2) P=y,y,y,...y,=f{z-M,)f{z-M^...fiz-M,). 

This probability P depends upon z the true value of the 
quantity, which is unknown. If we give to z various values 
we shall have corresponding values of P; and in the impossi- 
bility of finding the true value of z we can only find its moA 
probable value as given by the n observations, and the most 
probable value of z is that for which P is a maximum 
(Art. 9). To find the value of z which makes P a maximum 
we take the natural logarithm of each member of (2), giving 

(3) logP = log/(^-lfJ + log/(«-ig + ... 
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and by the usual rules find the value of z for which log P is 
a maximum. Diflferentiating (3) and placing it equal to zero, 
"we have 

6.P_ df{z-M^ Af{z-M ^ a^iz-iQ 

or since 



(5) 



df(z^M^)^<f>{z^M,)f{z--M:)d[z^M,), 

(in which j>{z — M^, ^{z— JfJ, etc: denote new fimctions of 
the errors) we have 

(4) ^^<i>{z'-M^dz + i^{z^M^dz + ... 

+ <^ (fe- M^ dz = 0. 

Hence the value of z which makes (2) a maximum is that 
which satisfies the equation 

(6) i^(«-Jf,)+0(«-^^+0(«-J/8)+... 

And this equation will furnish us with the most probable 
value of Zy provided that we can determine the form of the 
function <^. 

Now it is universally accepted as an axiom* that in direct 
observations made upon one quantity, the average or arith- 
metical mean furnishes the most probable result, that is, the 
most probable value of z in the case under consideration is 
found by taking the sum of M^, M^.,.M^^ and dividing it by 
n the number of measurements, or 

n 
This equation may be written 

which by transposition becomes 

(7) {z - M,) + (^- J/J + (^-1/3) + ... + («- ilfj = 0, 

• See Art. 66, p. 195. Digitized by GoOglc '. 
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that is to say, the arithmetical mean requiries that the 
algebraic sum of the differences or errors shall be zejr o. 

Comparing now (6) and (7) we see that the symbol 
means merely the multiplication by a constant, since the 
value of z must be the same in both ; hence 

(8) ^ (;5- IfJ + ^ (2r- ig +etc. 

= A; (« - JlfJ + i (^- JQ + etc., 

where h is any constant. Inserting in this the values of 
^ (2J — JlfJ etc. from (5), we have 

= i (^j - JfJ + i {z-M^ + etc. 

And since this is true for any number of observations, it 
must be true for one, or two, or three. Hence the corre- 
sponding terms in the two members must be equal If then 
X be any error, and y its corresponding probability, so that 

(1) y=fip\ 

we have from (8*) 

Integrating this we obtain 

(10) logy = ^ + ^, 

in which &' is the constant of integration, and the logarithm 
is taken in the Napierian system. Passing from logarithms 
to numbers we have 

(11) • y = e**^e*'=ce^,' 

in which e is the Napierian base, and the constant e*' is placed 
equal to c. Now y is to be positive, and is to decrease as x 
increases either positively or negatively, hence the constant 
h is negative. Placing then JA; = — A^ we have 

(12) y = ce-^'^ 

as the equatioh of the curve represented in Fig. l.)OgIe 



DEDUCTION OF THE FUNDAMENTAL PRINCIPLES. J25 

This equation satisfies the conditions imposed at the 
beginning of our investigation, for y is a maximum when 
a? = 0, it is symmetrical with respect to the axis of F, since 
equal positive and negative values of x give equal values for 
y ; and when' x becomes very large y is very small. The 
constants c and h will be particularly considered hereafter. 

116. Let us now consider the more general case of 
indirect observations in which the quantities 8yt.,,z are to 
be determined by measurements of a related quantity 
M=f{8,t,...z), Let n observations be made giving the 
values M^, M^.^. M^, and it be required to find from these 
results the most probable values of s,t f,, z. The measure- 
ments being imperfect the results Jfefj, M^, etc. cannot be 
perfectly accurate. Comparing each with the corresponding 
true value the diflferences will be errors (Art. 4) which we 
represent by x^, sc^y a;^...x^, each »of which is a function of 
(8,t...z\ Then from our equation (1) we have for the 
several errors 

(1) y.'/Ca'J, 2',=/K).-y«=/(-^«). 

And by Art.. 8 the probability of committing the given 
system of errors is 

(2) P = y,y,y,...y.=/k)/(a^,).../(aTJ. 

Applying logarithms to this expression it becomes 

(3) logP = log/(a?J + log/(a?J + ....+log/(^J. 

Now the most probable values of the unknown quantities 
8yt ... z are those which render P a maximum (Art. 9), and 
hence the derivative of P with respect to each of these 
variables must be equal to zero. Indicating the diflferentia- 
tion we have then the following equations : 

^P _ df(x,) dfja,;) dfix,) 

Pds -f{x,) ds ^/(x;) ds'^- "*"/(^J di - "' 
(*) 

Fdi -fix,) dt V W dr * • • "^ / W d< " "' 

etc, etc. 

M. Digitized by GoO^IC 
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the first being differentiated with reference to 8, the second 
with reference to t, the last to z, and so on, K in these 
we place 

(5) df{x^^4>{x^f{x^dx^, df{x;)^<l>{x;)f{x,)dx,,etc, 
they become 

(4) 

etc. etc. 

and being as many in number as there are unknown quan- 
tities they will determine the values of those unknown 
quantities as soon as we know the form of the function <f). 

Since these equations are general and applicable to any 
number of unknown quantities the form of the function 
may be determined from any special but known case. 
Such is that in which there is but one unknown quantity 
and the observations (8) taken directly upon that quantity. 
Thus if there be only the quantity s and the measurements 
give for it the values M^, M^.,, M^, the errors are 

from which 

dx^ __ dx^ _ _ dx^ _ 
ds ds '" ds "" ' 

and the first equation in (4) becomes 

(6) <^ K) + ii>{x^ + ^ {x^i + ... +1^ (o^J = 0. 

In this case also the arithmetical mean is the most 
probable value, and the algebraic sum of the errors will 
be zero (Art. 11a), or 

(7) a?, + aj, + a?3+...+a;„=-0, 
and equations (6) and (7) can only agree when 

(8) ^ (^i) +<A (^2) + ... + <^ W =^kx^ + kx^ + ... + A?(arJ, 

Jigitized by VjOOQ I 
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in which k is any constant. Eeplacing from (5) the values of ' . 
4> (a?J, <^ {x^, etc., it becomes 

and since this is true whatever be the number of observa- ; 

tions, the con-esponding terms in the two members are equal. 

Hence if ^ be any error and y =/(a;), we have "> 

(9) |^=^ = fc,. 
f[x)dx ydx 

Multiplying both members by dx and integrating, we obtain 

(10) logy = ^+&', 

which by passing from logarithms to numbers becomes 

(11) y^e\^^. 

As shown before the constant h must be essentially 
negative ; replacing it then by — SA'*, and also placing ^ ^ c, 
we have 

(12) y=ce-^, 

as the equation of the probability curve, x being the 
abscissa, y the ordinate, and c and h constants depending 
upon the precision of the observations. Considering c and 
h as unity the values of y corresponding to a few values 
of X are given in Part I. from which Fig. 1 has been con- 
structed, the vertical scale being double the horizontal in 
order to exhibit more clearly the form of the curve. 

The ordinate of the curve for x=0 is y = c, and hence 
c is the probability of the occurrence of the error x. The 
probability of committing any given error » is from the 
equation 2/ = ce"*'^, and this probability will be smaller as 
h is larger. Hence A is a measure of precision of the mea- 
surements. The more accurate the observations, the greater 
is A. 

12. For all kinds of observations, then, the law of the! 
probability of error is expressed by the equation 

(12) y = ce-'^A \ 

Digitized bv V^OQQ IC ^.^^ 
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which corresponds to the curve of Fig. 1. In discussing the 
properties and the consequences of the law, it is convenient 
to know the values of the constants c and K We proceed 
to find c (for A, see Arts. 25, 80). 

Let 00^, fl?j, x^;..x^ be a series of errors, fl?^ being the 
smallest, x^ the next following, and x^ the last, the £flfer- 
ences between the successive values being equal (thus if 
the measurement be made by a rule graduated to millimeters, 
x^ is V^ greater than x^, x^ is 1™™. greater than x^, and so 
on). Then by Art. 7, the probability of committing one 
of these errors, that is, the probability of committing an 
error lying between x^ and x^ is the sum of the separate 
probabilities ce~**^*, ce"'^^, etc., or if P' denote this sum 

(13) F = c {e-^^' + e-*^«' + e-*«V + . . . + e'^^'n^, 

which may be written 

(13) P' = c2^e-^, 

which denotes the sum of the probabilities of all the errors 
from x^ to x^ inclusive. P' denotes then the probability 
that an error lies between the limits x^ and x^. Now if 
i denote the small interval between the successive values 
of X, and if our observations are accurate enough so that 
X may be regarded as a continuous variable, i will be veiy 
small and equal to dx. Then 

(14) it;;e-^=re'^'^dx, 

* J Xi 

from which by comparison with (13) we have 

(15) • F=rcX:^e'^^^J\-^dx, 

which expresses the probability that an error will lie between 
the limits x^ and x^. Now it is certain that the error will 
lie between — oo and + oo , and as unity is the symbol for 
certainty (Art. 5), we have 

(16) l = |J^V*^tfe, 

Digitized by LjOOQIC 
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from which the value of the constant c will be known as 
soon as we find the value of the integral between these 
limits. The integral calculus gives * 



(17) 






>^is=^. 



* The most conyenient method of determining this integral is the geo- 
metrical one dne to Poisson. I give the process nearly as presented hy 
Stvbm in his Cows d^ Analyse, Paris, 1857, Vol. ii. p. 16. 

Since y—e~^^^ is the equation of a carve of the same form as shown 

in Fig. 1, the integral jydx^jer^^^ expreeses the area between that curve 

and the axis of x ; and since the curve is symmetrical to the axis of Y, that 
integrtd between the limits - oo and + oo will be equal to doable the integral 
betjireen the limits and + <x> * Placing also Ax= t, we have 



C'-*»^^=lj\-^dt, 



and we have to determine the integral in the second member. 

Fig. 15. 




If we take three co-ordinate rectangular axes OT, OU, and OT, and 
change t into «, we have 

A = Je't^dt = area between curve VtT and axes, 
4=je~«^(itt= area between curve Fu [7 and axes, 
and A^= r r (Td^v^dtdu. ^ t ^ 

•'O •'<> Digitized by COOgle ^j|| 
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and hence from (16), 

/io\ t ^s/'"' ^* 

inserting this in (12), we have for the equation of the pro- 
bability curve, or the law of the probability of error, 

(19) y^hiTT'^e'^^, 

in which x is any error, y its probability, h the measure of 

precision (a quantity so that t is a concrete number of the 

same kind as x), and i is the smallest graduated division 
in the instrument of measurement^ and hence a quantity 

of the same kind as a? or t^ . The probs^bility y is then an ab- 
stract number^^ as of course it ought to be. 

Inserting also the value of c in equation (15), we have 
h f*- 



(20) F=^-^re'^'^'dx, 



which expresses the probability that an error will fall between 
the limits x^ and a?^. Also since the integral between the 

Now v=c~<^ is the equation of the ctirve YtT and t;=e"«* is the equation 
of VuU^ and if either of these curves revolves about the axis of V it gene- 
rates a surface whose equation is i;=«-t3-u«. Hence the double integral A* 
is one-fourth of the volume included between that surface and the hori- 
zontal plane. If we suppose a series of oyUnders concentric with the axis V 
to form the volume, the area of the ring included between two whose radii 
are rand r+dr is 2'nrdr, and the corresponding height is v=e^^u=e~f^. 
Hence one-fourth of the volume is 

which, since /i~«*2rdr = - «"**, is equal to ^. Therefore we have 

and hence, finally, 

(17) f^e-l^^dx = I re-i^dt = ^. 

Digitizfed by VjOOQ iC 
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limits — X and + a? is twice the integral from — a: to or from 
to + a?, we have 

(21) F = A f ^""e"*'^^ = ^ \\-^dx, 

as the probability that an error taken at random is between 
the limits — x and + x, or is numerically less than x. 

Now (19) is the equation of the probability curve, and the 
area between the curve and the axis of a? is 



jydx == -i je-^'^dx. 



J' 

Hence if (16) be multiplied by i it will be the total area of 
the curve, and if (21) be multiplied by i it will be the area 
between the limits — x and + x. Hence expressions (16), 
(20) and (21) are proportional to the areas of the probability 
curve corresponding to those limits, and if we regard the 
total area as unity, the partial area between the limits — x 
and + a? will be a fraction given by P in equation (21). 
Further, since errors are committed in proportion to their 
probabilities, these integrals and their corresponding areas 
are proportional to the number of errors which we should 
expect to find between those limits. If then we compute 
values of P' corresponding to successive numerical values of a? 
in equation (20) they will be fractions proportional to the 
number of errors numerically less than a?, and at the same 
time express the probabilities of committing an error less 
than X. As however the constant h depends upon the pre- 
cision of the measurements, and hence varies in different sets 
of observations, we write equation (21) under the form 

(22) ' Fr.^re-^dM~\^e-m, 

and compute the values of P corresponding to successive 
numerical values of hx or t, by the usual methods of the 
integral calqulus*. A table of these values is given in 
Part I. 

• Developing ert^ into a series by MoLaubin's theorem, multiplying by 
dt and integrating, we get r^^^^r^T^ 

Digitized by VjOOQ IC 
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13. Thus from our table we see that the probability of an 
error corresponding to Aa? = oo is 1, that is, it is certain that 

all errors will be less than x = y.. For hx = 1*24 we have 

P'ss 0'9205, that is, the probability that an error will be 

1*24 9205 
committed less than x = -r— is ^ ;^,x^\. Or in other words, if 

fl lUUUU 

we have 10000 observations we sliould expect that in 9205 

of them the errors would be less and the remaining 795 

1*24 
greater than -y- • 

The Principle of Least Squares, 

14. The law of the probability of error represented by 
the equation of the probability curve leads directly to import- 
ant results. The equation is 

(12) yr^ce-^, 

in which y is the probability of committing the error a?, and 
c and h are constants independent of x. Now, considering 
the general case of independent indirect observations (Art. 
19), let w equally good observations be made upon a quantity 

^^^ (:' •- 3 + O- 6 - 17273 • 7+^S • 

which is convenient for small values of t. For large values we integrate by 
parts, thus 

1 « 1 -« . 3 r^-^ J. 1 

"=-2l'"'^ + 2^' «»+25/-ir^« = et°-' 
e~^dt = ^^ as shown in the preceding footnote, we have 

From these two series the valnes of P' can be found to any reqtxired 
degree of accuracy for all values of t or fcc. r^ ^ ^ ^ 

Jigitized by VjOOy I 
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M—f{Sy t.,.z) for. the purpose of determining the magnitude 
oi8yt,.,Zy and let the observed results be if ^, if ^ ... MJ each 
of which is a certain function oi 8,t... z. Subtracting each 
observed value from the corresponding true value there 
result the errors x^yX^,,,x^y having the respective proba- 
bilities 

(12) y,^ce'^^\ y^^ce-^^\ y=ce-^\ 

Now, by Art. 8, the probability of committing all these 
errors either simultaneously or successively is the product 
of the separate probabilities, or 

(23) P = ce-^^^ ce'^^ ce-'^' = c^^-'^^', 

in which 2ic* denotes the sum of the squares of the n 
errors, or 

(24) 2a? = a;,V<+a?3»+...+a:,». 

Now each of these errors is a function of the quantities 
s,t,..z, which we are to determine. Their true values we 
can never be sure of having obtained however accurate be 
the measurements. Moreover, in equation (23) the proba- 
bility P will vary with these quantities, or as %a? takes 
different values, and hence out of the many systems of values 
which many be assigned to s,t ... z, we must take the most 
probable as approaching nearest to the true system, that is 
the system for which P is a maximum (Art. 9), and P is a 
maximum when 2a?' is a minimum. Hence the most probable 
system of values of observed quantities is that which renders 
the sum of the squares of the errors a minimum. 

This is the fundamental principle of least squares for ob- 
servations of equal precision. If they are not equally good, 
the constants c and h are different in each observation or set 
of observations, and we have 

(12) y^=cj&-^K y,= c,e-'^«''«' y»=c«e-'^^' 

as the respective probabilities of the errors ar^ , a? ^ , . . . a?^. Then 
as before the probability P of the total system of errors is 

(25) P=CM...c^e-<^^'-'^'^^--'^^2,^Coog\e 
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and the most probable system is that for which P is a maxi- 
mum, or that for which 

(26) 2AV = h^x^ + h^x^ + . . . + h,^x^ = a minimum. 

Hence in measurements of unequal precision, ih^ square of 
each error must he multiplied by the square of its measure 
of precision and the sum of the products be made a minimum. 
If h is the same for all* errors this evidently reduces to the 
rule previously given. Thus arises the term "Least Squares." 

15. Thus if all the observations are equally good and 
made directly upon one and the same quantity M whose true 
value is z, the errors are 

ajj= z — M^, x^ — z — M^,.,x^=^z — M^ 

and the expression (24) is to be made a minimum, or 

(24) 2 (« — MY = a minimum. 

By the usual process for determining minima this gives 

2 (^ — M) = 0, from which z = , 

n 

which is merely the expression of the law of the arithmetical 
mean, or of the equation (7), which we took as the founda- 
tion of our reasoning. 

The Measure of Precision and the Probable Err(yr. 

16. Let us consider two series of observations, one having 
the measure of precision \ and the other h^. The probability 
that a single error in the first series will lie between the 
limits —a;^ and '\-x^ will be expressed by the integral 

2K 



(21) P;A^e-^'-'dx, 



and the probability that an error in the second series will fall 
between the limits — x^ and + x^ is also 

(21) ^;=^ e-'^'^dx. 

Digitized by VjOOQ IC 
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These two integrals are equal when \x^ = h^x^. If the first 
series is three times as precise as the second, we have \=^ 3A,, 
and hence the integrals will be equal when x^ =» 3^7^, that is, 
the probability of committing an error less than x^ in the 
first series is the same as that of committing three times as 
large an error in the second series. Hence the accuracy of 
different sets of observations is directly proportional to their 
measures of precision. 

Owing however to the circumstance that A is a quantity 

of the same kind as - , and hence expressed in terms of an in- 

% X 

convenient unit, it is usual to employ other constants for the 
comparison of the accuracy of sets of measurements. The 
one m most common use is called the probable error, which 

is an error of such a value that the probability (22) is ^ : and 

is hence an error such that it is an even wager that an error 
taken at random will be greater or less than it. The proba- 
ble error is then the value of x given by the equation 



1 2 r** 
(27) i^4= e-^-'dJix. 

^ JttJ 



By interpolation from the table in Part I. Art. 13, we 
find 

for A^= 0-4769, P' = 0-5; 

hence denoting this value of x by r, we have for the probable 
error 

(28) hr = 0-4769 or r = ^J— . 

Hence if in Fig. 1 we lay off the abscissae — OP and + OP 
equal to the probable error r, and draw the ordinates PB and 
PB, the area PBABP will be one-half of the total area of 
the curve, and in any series of observations we may expect 
that one-half of the errors will be less numerically than r, 
and the other half greater than r. 

17. If then we have two sets of observations whose 
measures of precision are h^ and h^ and probable errors r^ 
and r,, we have 

(28) V. = 0-4769. V, = 0-476|^^^^^^^CoogIe - 
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and hence h^r^ = \r^. Now if the precision of the first series 
is three times that of the second A. = 3Aj, and hence r, = Sr^, 
that is, the probable error of the first series is one-third that 
of the second. Hence the probable error serves to compare 
the accuracy of measurements equally as well as measures of 
precision. The smaller the probable error, the better are the 
observations. Thus if two sets of observations give for the 
length of a line in centimeters 

i,= 427-32 ±004 and X,= 427-31 ± 016, 

in which 004 and 0*16 are the respective'probable errors, the 
meaning is that it is an even wager that the first is within 
0*04 of the truth, and also an even wager that the second 
is within 016 of the true value ; and the precision of the 
measurements in the first set is four times that of those in 
the second. 

We have now given the fundamental theory of least 
squares and probable errors, and shall proceed in the next 
Chapter to develope its practical features. 
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CHAPTER VII. 

DEVELOPMENT OP PRACTICAL METHODS AND FORMULA. 

18. We distinguish the following kinds of observationfi : 
Direct observations upon a single quantity, in which the 
measurements are made directly upon that quantity, 

19. Indirect observations upon one or more quantities 
s,t ... zhj the measurement of functions of those quantities 
(Arts. 11& and 14). 

20. Conditioned observations, which considered singly 
are independent, but which collectively are subject to vv^ot- 
ous requirements or conditions. They may be either direct 
or indirect. , 

21. Independent observations, which are also either 
direct or indirect, but between which there exists no con- 
ditional requirements. 

Thus if the sides and angles of a field are measured, each 
observation taken alone is direct. If we find its area from 
the sides and angles the measurement of that area is in* 
direct. Further, any two sides considered are independent 
of each other, but if we consider all the sides and anglea 
they must fulfil the condition that when plotted they shall 
form a closed figure. 

Direct Observations upon one Quantity, 

22. We take up first the case of direct observations 
of equal precision upon one and the same quantity. Tiiese 
are combined by the use of the principle of the a^ crage 
or the arithmetical mean, whose use is limited to this single 
class. 
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23. As stated in Art. 11a, the most probable value of 
a quantity which is measured n times with the results if^ 
Jfj... M^ is the arithmetical mean, or 

^'^ ^ ^ n n 

From the second fundamental axiom of Art. 10, which 
asserts that in a great number of observations positive and 
negative errors are equally probable, it would seem that 
the average was also the true value of the measured quan- 
tity. As far as our observations show, this is the case, for 
it is the most accurate value deducible from them, and must 
be used as if it were the true value. For an infinite number 
of measurements the average would be the absolute true 
value; for a limited number it can only be regarded as 
the most probable value, that is, as the nearest approxima- 
tion we are able to make to the true value. 

24, Probable errors. Having taken the average of n 
equally precise measurements upon a single quantity, we 
next proceed to investigate the accuracy or precision of 
the result. 

Let the observations give the values Jf^, M^,..M^ whose 
mean is z^y and whose corresponding true value is z. Let 
the errors committed in the several observations be 

ajj, a?jj ... a?„, 
so that 

(30) Xj^^z-M^y ajj = «-ilfj,...fl?„ = 2r- Jf^. 

Also let v.yV^.,.v^ be differences resulting from sub- 
tracting each observed value from the arithmetical mean, or 

(31) v, = z,--M,, v, = ^,-Jlf,...v„ = ^,-Jf^. 

If Zq were the true value of the quantity, or Zq=^z^, then 
the errors x would be the same as the differences v. But 
as we can never be sure that z^ represents the true value, 
we can never determine the errors a?^, x^,.,x^. The values 
Vj, Vjj ... v^ which are readily found m any particular case, 
we call residuals; they should be carefully distinguished 
from errors. ^ 
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Each of our measurements Ifj, if, ... Jf^ is probably 
incorrect, a« likewise the arithmetical mean z^. The degree of 
confidence which we can place in each of these results will 
be shown by their probable errors (Art. 16). Considering 
h as the measure of precision of a single observation and 
r as its probable error, the relation between them is given by 

and hence to determine r we have only to find h. There 
is, however, no known method of finding the exact value 
of h : the best that we can do is to determine an approxi- 
mate value, which moreover shall be the Tnost j>robable one 
(Art 9). 

25*. The probability of the occurrence of any error 
xia 

(19) y = Ai7r-»e-*^, 

and the probability of the occurrence of the system of errors 
iPj, iCj ... a?^ is by Art. 8 the product of the separate pro- 
babilities. If then h is the same for each of the n observa- 
tions, we have as in Art 14, 

(32) P-y,y....y. = AV7r-i«e-'^, 
where Saj' denotes the sum 

Kow in this expression h is unknown, and further, we have 
no means of findihg its exact value. But whatever be its 
value, P must be a maximum in order to give the most 
probable value of the measured quantity z. We are there- 
fore led to conclude that for a given system of errors the 
most probable value of h is that which renders P a maxi- 
mum. Diflferentiating (32) with reference to h, and putting 
the first difierential coefficient equal to zero, we have 

(33) ^ = nA""S^7r-i»e-'«*'-2A2a?'e-^^AV7r-»^ = 0; 

* Simplified and considerably altered from the demonstration given by 
l>iBKO£B in his Avsgleichung der Beobactungsfehler^ Braunschweig, 1857, 
P* 59. 
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dividing this equation by A*"^, i*, tt"** and e"^*^, we obtain • 

(33) n-2A*2a?' = 0, 

from which 

(54) 2a? = ^.andA«=2^. 

We have thus the value of h in terms of Xa^, which, 
however, we have no means of obtaining, since the errors 
ajj, ajj ... a?^ depend upon the unknown true value z. If 
the number of observations were infinite 'S,a^ would equal 
Xv* (Art. 23), and as the latter value is determinate A would 
be known. In a large number of observations, therefore, 
the equation 

A'= ^ 



2St;» 

will always give a close approximation to the value of 
h. But as the sum Sv' is always less than 2a^ (since from 
the principle of least squares, Art. 14, the first is the mini- 
mum value of the second) we may place 

(35) 2a? = 2v' + A;', 

in which ^ is a constant to be determined, and then our 
value of h will be correctly given by inserting in (34) the 
value of 2a?' from (35). As however 2^?^ cannot be exactly 
found, we cannot hope to find the exact value of A*, but 
must be content with determining ap approximate one. 

Now the probability of committing the system of errors 
a7j, iPj ... a?^ is 

(32) P = c*e-'^«^, 

or, inserting for 2«i' its value from (35), 

(36) P^^c"" e-^^"^-^^^ = c* e-^'^'^e-^*^. 
Placing in this the constant terms equal to c, it is 

(36) P = ce'^. 

Hence the law of the probability of any value k is the same 
as that of an error x, as shown by (12). We may regard then 
(36) as the equation of a curve of the same form as Fig. 1, 
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and, as in Art. 12, we may show that c = hiir^^, where i is a 
small constant of the same kind as k Hence the equation 



(37) P = Ai7r- e-^'^ 

shows the probability P of a value k. Now in this equation 
both h and k are unknown, and, as we have said before, we 
can only expect to determine their most probable values. 
The value of A' although unknown is fixed and definite, and 
hence we conclude as before that the most probable value of 
h will be that which makes P a maximum (Art. 9). Differ- 
entiating then the equation (37) with reference to h and 
placing the first differential coefficient equal to zero, we have 

(38) ^ = iir-i e-^ - 2hk'e'^ hiTr'^ = 0. 



Dividing this by i, tt"' and e-*'**, we have 


(38) 


l-27t'A'=0, 


from which 




(39) 


*^ 2A'' 


Therefore for the equation (35), we have 


(40) 


2«^=2.« + i 



as the nearest possible approximation. Since from (33) the 

Til 

value of Sa;' is gp, we have accordingly^ 

(^). 21^ = ^'^ + ^- 

from which we find the most probable value of h in terms of 
the known sum of the squares of the residuals (31), or 

(41) ^^^VW- 
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Inserting tliis value in (29) we have as the probable error of 
a single observation 

(42) r = 0-4769y^=0-6745y^. 

The above is only one of the many demoastrations of the 
formula for the probable error r. None of them are, from 
the nature of the case, entirely satisfactory, since it is im- 
possible to find the exact value of h. In Art 30 we give 
another demonstration which may be readily applied to the 
case here considered by making ^r = 1, or regarding the 
measurements as of equal precision. 

26. The formula (42) just deduced gives the probable 
error of a single measurement M. We next inquire, what is 
the probable error of the arithmetical mean of the n measure- 
ments ? 

The probability of the arithmetical mean is the proba- 
bility of committing the system of errors v^, i?^ . . . t;„ or 

(23) P,= c"e-'^'^^, 

and the probability that the true value of z is, z^-^ x^ is the 
probability of the system of errors v^ + x\ v^ -f a?', etc., or 

(43) P,^ = c"e-^(-+»)». 

Since S (v + xf = tv* + 2x tv + nx*^ and tv = 0, this be- 
comes 

(43) p^, = c«^-ft2(Zr»+«r'2,^ 

Hence we have 

(44) P, : P^ :: e-^"^"^ : 6-'^«(^«=-»-««'^ :: 1 : e-«^^, 

that is, the probability of the error in the arithmetical 
mean is to that of the error £c' as 1 is to e-"'^'*. For a single 
observation whose error is x, we have however from (2) 

(45) • J/o' !/.::! : e'^"^, 

or the probability of the error in a single observation is to 
that of the error x' as 1 is to e*"'^^^'^. Hence comparing (44) 
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and (45) we see that if h is the measure of precision of a 
single ohservation, sjnh must be the measure of precision of 
the arithmetical mean z^. Therefore, h^ denoting the measure 
of precision of 5?o, we have ^ 

(46) h^^h'/ii. 

Denoting by r^ the probable error of z^, 

(28) Vo =0-4769 and Ar=^ 0-4769. 

Inserting from these the values of h^ and h in (46) we have 

Vq = -.:=-, that is, the probable error of the arithmetical mean is 

wn 
equal to the probable error of a single observation divided by the 
square root of the number of observations. Hence from (42) 



(4t) r. = ^ = 0-6745y^^|^^-^. 

which is the expression used without proof in Part I. 

27. We take up next the adjustment of direct ob- 
servations' of unequal precision. If the several observations 
give the results M^,M^.,. M^ whose measures of precision are 
h^,h^...h^, the principle of least squares requires that the 
quantity 

(26) A, V + KW + — + KW = a minimum. 
If j? be the true value of the quantity the errors are 
a?i = 2? — JI/j, x^ = Z''M^ ... OD^^^z — M^, 
and the expression to be made a minimum is 

(26) h,^{z-MJ + h/{z-M,r + ...+h:iz-MJ. 

By differentiation, we find that the value of z which makes 
this a minimum is 
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1 



Owing however to the fact already alluded to in Art. 16, 
that h is expressed in terms of an inconvenient unit it is 
usual to employ numbers having the same ratios; thus if 

(49) g,:g,'. 9, :: K ■ K • K. 

the expression (48) becomes 

The numbers g^^g^ ..»gn are called the weights of the ob- 
servations M^yM^,..M^y and are merely relative numbers pro- 
portional to the absolute quantities AAA/ ... h^. If, as in 
Art. 23, all the observations are of equal precision, we take g 
as 1, and the formula (60) agrees with the law of the arithme- 
tical mean (27). The weights of measurements are then 
numbers proportional to the number of single observations 
to which each is equivalent. Thus if M^ is the equivalent of 
g^ .single observations of the weight unity its weight is gr^, 
and as the arithmetical mean z^ is the equivalent of n single 
observations its weight is w. The weight of Z, which we call 
the general meariy is then 5^i + 5^, + ... + jr^. 

28. If we have different sets of measurements of unequal 
precision upon one quantity, the first giving the average z^ 
from rij measurements, the second giving z^ from n^ measure- 
ments, the adjusted value is furnished by 



(50) 



2^' 



and the weights gi, g^'^gn i^^st be found from the propor- 
tion (49). The measures of precision of the several averages 
being h^, \, etc., we have from equations (41) and (46) 

in which Sv* denotes the sum of the squares of the residuals 
in the first set, 2v"* in the second, and so on. Inserting these 
in (49) and omitting the common factor 2, we have 

f52^ a ' a -a - ^>^^ • ^«K^^) • ^K^^^ 
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from wliicli we may find the relative weights, and then the 
general mean by (50). The weight of the general mean is 
of course 2^* 

29. If the probable errors of the averages of different 
sets of measurements have been found by (47), their relative 
weights are easily determined by (49). For if the measures 
of precision be h^,h^,..h^, and the corresponding probable 
errors be r^, r^ ... r^, we have 

/OQN 7^. 0-4769 , 0-4769 , 

(28) K' = —r- > K=' —zr- ^ etc.; 

'i 'a 

inserting these in (49) and omitting the common factor 
0*4769, we have 

/RQ^ 111 

that is, the weights of observations are inversely proportional 
to the squares of their probable errors. 

30. We now proceed to find the probable error of the 
general mean. Let Z be the general mean, O its weight, 
and R its probable error, also let r^ and r^ be the probable: 
errors of observations of the weights g^y g^. Then from the 
above principle 

(53) G I g,: g,y.-j^:^:^,, 
from which we find 

(54) i?=-?^L=-?^ = etc. 

Also since ^g^+g^ + ... + fl^„, we have from (53) 



1=1+1+ +1 

M r. r^ r^ 



(55) pa""xra'^r^"^'"''^iri* 



Hence, having by (46) found the probable errors r,, 
^2 ••• ^n> the probable error of Z may be found by (55), or, 
having by (52) found the weights /7i, 5^, ...fl^n> whose sum 
is G, we may by (54) also find the probable error^ Z^ I 
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As however the computation of all the weights or all 
the probable errors is sometimes laborious, we shall deduce 
another formula. Let, as before, O be the weight of the 
general mean, and R its probable error, and let r be the 

5robable error of an observation whose weight is unity, 
'hen we have 

(53) ^:i::^:l, 

from which 

C5«) B.^. 

Hence having found the probable error r of an observation 
whose weight is unity ^ the probable error of a result whose 
weight is G is found by dividing r by the square root of G. 
We proceed to develope a method for finding r. 

Let n be the number of observations or sets of measure- 
ments, A„ A, ... A» their measures of precision, and jr,, g^^.^.g^ 
their relative weights. Also let h be the measure of pre- 
cision of an observation whose weight is unity. Then 
from (49) 

and hence 

(57) K=ff,K V =ffJ<'* '•' K'=ff^*- 

The law of the probability of error gives for the corre- 
sponding errors x^,x^...x^, 

(19) y, = V>"*«"'^'*^'> y^^K^ir-W^"^" e^n 
or from (57), _ 

yi = A J 9, \ ^-* ^"''''^^^^ y, = '^ J 9, \ ^-* e-^^'"' 

Hence in general if x be any error, g the weight of its 
corresponding observation, and h the measure of precision 
of an observation qfthe weight unity ^ we have 

(58) y^hjgiir-^e-^^, 

as the probability y of any error a?, whose measure of pre- 
cision is hjg, or as the probability of any error xjg whose 
measure of precision is A, r^ t 
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Now n observations being made, it is evident that the 
quantity 

^ ^ n n 

will have a certain definite value; the probability of the 
occurrence of g^x^ will be y^ of ff^os^ will be y^, etc. : and 
each term in the numerator will occur a number of times 
proportional to its probability, provided that n is a very 
large number. Hence g^x^ occurs ny^ times, g^^x^ occurs 
ny^ times, and the quantity (59) becomes 

(59) g,x^\ + g^x^^y^ +...+ g^x,\ = %^y. 

But whan » is a very large number, the errors will be 
distributed according to the law of the probability curve 
from —00 to +oo, and if the measurements are accurate, 
i in (58) will be dx. Henee we have, by inserting for y its 
value, 

(60) ^ = ^gx'y « V^ f ^ V e-^'''' dx. 

n ^ir J -00 

Taking in this hxjg=^t as the unit variable, it may 
be written 

(60) ?2^*== i.r>^-^^^, 

and as the value of the integral is -^- * we have 
* From the foot-note to equation (17) we have 



■p 



,^ 



*-# 



Placing t^tijly this becomes 



I 



>.dt=^.. 



JO w» 

Differentiatiiig this egnation mth reference to $, and regarding t as con- 
stant, we have 
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which gives the value of h in terms of the sum Xgo^. 

Now let Vj, v^.,,v^ denote the residuals or differences 
between the general mean Z and each observation. Then 
the sum Sjri? is greater than 2grt?', since the second is 
the minimum of the first. If we place then 

(62) 2^-S^«;* = ifc^ 

and suppose S^raj* to have all possible values greater than 
S^ri;*, and each to be repeated a number of times proportional 
to its probability, we may consider the mean of all the values 
thus found for i* as the best approximation a,ttainable to 
its value. The law of the probability of these values of k 
is as in (36) given by the equation. 

and if n be the number of possible values of h^ the value 
k^ will occur nY. times, k^ will occur nF, times, and hence 
the mean of all the possible values will be 

(63) tk'Y= ^r^k'e'^'^dk = ^,. 

If then we place in (63) for Xga? its value from (62) and 
for Ar* the value just found, we have 

(63) 2T' = 2^^ + i' 

from which we find for h 



m A=,/| 



-1 



Hence from (29) we have the probable error of an ob- 
servation of the weight unity 



/^rx 0-4769 rt^^.r /2flr«* 

(65) r = — ^=0-6745 y'^. 

Diyiding this hy -ds and making «=1, we obtain 



{. 



I c-^ t^dt = y^ = one-half of the integral in (60). 
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Whicli being computed, the probable error of the general 
mean is found at once by (56). 

If in (64) and (65) we place ff = l, they reduce to the 
expressions (41) and (42) obtained in Art. 25 by an entirely 
different method. Moreover in that case G becomes n, and 
the expression (56) coincides with (47). Both methods of 
reasoning thus lead to the same results. 

31. The above includes the whole theory of direct 
^measurements upon a single quantity. Observations upon 
several quantities will be investigated in the next section; 
and it is interesting to observe that all the preceding methods 
are but particular cases of the more general theory of in- 
direct observations, as we shall show in the following articles. 

Independent ohservations upon several quantities. 

32. Independent observations are those which are sub- 
ject to no conditions except those imposed upon them by 
the measurements themselves, so that, before taking the 
observations, all systems of values are equally probable. 
The manner of measurement may be either direct (Art. 18) 
or indirect (Art. 19), but for convenience we shall consider 
only the latter, of which the first is a special case. 

Indirect observations being made upon functions of the 
quantities to be determined, require in general the statement 
of equations between the measured quantities and those re- 
quired. Thus if in order to determine the magnitude of the 
quantities s,t...a we make observations upon the related 
quantities Jif, Jf, which are connected with the first by the 
relations 

M=f(s,t...z), 

3f=f(s,t...z), 

and find the values M^, M^, Jf/, M^\ etc., each observation 
furnishes us with an equation, whicn we call an observation 
equation. The number of these equations is the same as the 
number of observations, and generally greater than the num- 
ber of unknown quantities 5, ^ . . • -e^ which we are to determine. 
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Hence in general no system of values can be found for s, ^ ... 2 
which will exactly satisfy the observation equations. They may 
however be approximately satisfied by many sytems of values, 
and we propose then the problem to find out of these systems, 
all equally possible, the. one which is the most probable, and 
hence the best. 

The equations between the observed and unknown quanti- 
ties may be either linear or non-linear, exponential or trans- 
cendental ; but we shall treat only rf linear equations, to 
which all the others can always be reduced by the methods 
of Art 69. 



S3. Taking up first the case of observations of ( 
weighty let the equations between the observed and the 
measured quantities be of the form 

as + bt-h .,. + lz==M, 

in which Syt...z are the unknown quantities to be deter- 
mined, a, J...?, constants given by theory and absolutely 
known, and M the measured quantity. For each observation 
we shall have a similar equation, and in all the following n 
equations 

a^s + bJ;-^...'^l^^M^ 
(66) ii^s-{'b,t'i^... + l,z^M^ 



the first of which arises from the first observation, the second 
from the second, and the last from the n^« 

Now as the number of these observation equations is 
greater than that of the unknown quantities, they will not 
be exactly satisfied for any system of values we may find. 
Hence if 8,t .,, z denote the true values of the quantities to 
be determined, the general form 

as + bt-^-.-. + U-M^s^O 

is only approximately correct, for the equations do not re^ 
duce exactly to zero. Let us designate then by a?., f, ... a?., 

jdgle 
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the errors whicli thus arise when (or 8,t ...z in. (66) we place 
their true values. Then we have strictly 

a^s + &i^+ ... + ?!«- J^i = iri, 
ajs + 5J + ... + Lz — -Jf = a?j, 

(67) 

a^8 + 6^* + .•. + l^z — M^= x^. 

Now in the impossibility of finding the true values of 
8yt...z from these equations, we must determine their most 
probable values as the nearest attainable approximation to 
the truth. The most probable system of values is, by the 
fundamental principle of Art 14, thatwhich makes the sum 
of thie squares of the errors a minimum, that is which makes 
ajj* + a?/ + ... + x^ a minimum. By the use of this principle 
we have in Part I. deduced a method of finding the most 
probable values. We give here a more general proof which 
follows directly from our demonstration of Art. 11 J. 

From equations (8) and (9) of that Article we have 

<j> (a?,) = Icx^, <f> (a?J = hx^ ... ^ (a?J = kx^, 

in which k is any constant. Substituting these in the differ- 
ential equations (4) and dividing each by k, we obtain 

dx. , dx^ . ^^8 . i ^^n A 

(68) ar,-^-' + a-,-^«+ar.-^«+... + a..-^=0. 



Now by differentiating equations (66) with reference to each 
variable, we obtain 

dx^ _ dx^ _ dx^ _ 



dx^ _^ J dx^ __ , dx^ ^ J 
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which substituted in (68) give 



(70) 



^1^1 + k^% + K^z + ••• + i^x^ = 0, 



which are the conditions for determining the most probable 
values of 5, t...z, since they render equation (2) a minimum, 
and which are as many in number as the number of those un- 
known quantities. K in these we substitute ioi x^fX^\..x^ 
their values from (67), we have then the equations from which 
8,t...z can be determined. These final equations we call 
normal equations, and we see by (70) that the first is formed 
by multiplying each observation equation by the coefficient 
of 8 in that equation and adding the results. Inserting in 
the first of (70) the values o{ x^x^.,. x^irom (67), we have as- 
the first normal equation, or the equation for 8, 

(71) (a^^ + a^^+... + a^*)8+(aJ>,+ aJ)^+... + aJfJt + ... 

+ {a,l, + a,l^+.,.+aJ^)z 
- (a^Jf, + a,;»f, + ... + a„JlfJ = ; 

and in like manner we form a normal equation for each of 
the other quantities t,u ... z. To abbreviate the expression 
of these equations let us place 

2 a» = a^" + a,* + 03' + . . . + a^', 
tab = afi^ + ajb^ + ajt, + . . . + aj)^, 

(72) Xal = aj^ + aj^ -^aj^^^ ... + aj^, 

2aJf = a^M^-h a^M^ + a^M^ + ... + a.il/^, 
etc. etc. etc. 

and then the normal equations (71) may be written , 

. oogle 
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2a* . 5 + SoJ . ^ + ... + 2ai. « = 2a3/, 

(73) ^ac.8 + 'Zbc.t + ...'\-XcLz=XcM, 



Xal.8 '+ 2W. ^ + ... + 2P . -^ = tlM. 

The coefficients of the unknown quantities, it will be re- 
marked,, present a curious symmetry; thus the coefficients of 
the first horizontal row are the same its those of the first 
vertical row, those of the second horizontal row the same as 
those of the second vertical row, and so on. The period placed 
between each coefficient and its unknown quantity shows that 
the sign 2 extends only to the former. 

34. Thus if we have n observations for determining q 
unknown quantities, the most probable values of the un- 
known quantities are obtained by writing n observation 
equations as in (66), then forming the q normal equations as 
in (73), whose coefficients 2a*, 2aJ, etc. are given by (72) ; 
then the solution of these normal equations will furnish the 
most probable values of 8,t ... z, -In the most common 
cases the coefficients in the observation equations (66) are 
+ 1, — 1 or 0, and in the formation of the sums (71) the signs 
must be carefully regarded. Stated in words, the process for 
forming the above normal equations is the same as given by 
the rule deduced in Part I, 

35. To adjust indirect independent observations of equal 
weight, we have then only to form and solve the normal 
equations, thus obtaining the best system of values for the 
unknown quantities. 

The solution of the normal equations may be effected by 
any algebraic method. When there are only two or three 
equations the usual methods of substitution or addition are 
perhaps the quickest, but for many equations they are tedious. 
A process of solution by substitution, assisted by the notation 
of (71), which is due to Gauss, is here valuable to the com- 
puter (see Art. 60). The method of indeterminate multipliers 
is likewise often of quick application. As this method is not 
presented as fully as it ought to be in the common text-books 
in Algebra, we think it worth while to give an example illus^ 
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trating its use, particularly as we shall have occasion to refer 
to it again in another connection. 

We take the five normal equations given in Part L Art. 
35, viz. 

28-1 ^A, 

— * — u + Sx— y =«J9, 

and to solve them by the method of indeterminate multipliers, 
we multiply the first by a number yS, which is as yet un- 
known, also the second by a number /S,, the third by yS,, the 
fourth by /8^, and the fifth by /8,. We then add the resulting 
equations and place together the terms containing like un- 
known quantities, thus 

Now if we wish to find the value of y we have only to require 
such relations to exist between these multipliers that all of 
the terms of the first member shall disappear except that 
containing y, that is, we must have 

2/3, -i8, =0, 

-A + 4/3,-/33-/S, =0, 

-/3, +2^3-/3, =0, 

^d then we find 

_ ^,A + fi,B + ^ ^C+l3,D + fi,E 

Now /Sj, /8j, etc. may have any values which will satisfy the 
four imposed conditions. If then we take ^ = 1, the first of 
these conditions gives y3,= 2^ and the other three become 
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7-/S3-/9, «0, 

-2+2/3,-/9, =0. 
-2-/9. +3/9,-/5. = 0. 

Solving these equations by the common algebraic method of 
addition we find /Sg = 3, ^^ = 4, and fi^ = 7. Hence our value 
of V is 

^+25+3(7+40+7^ 
!/ 17 ' 

In like manner if we wish to find the value of s, we re- 
quire that all the terms shall disappear except the first, or 
that 

-/9, + 4/3,-/9,-/9. =0, 

-/S, +2/3,-/3, =0, 

- /3, - /3. +3/3,-/8, = 0, 

-/3.+3/3,= 0, 
{tad then we have 

•- 2/3,-/8, 

and as the multipliers may have any values which will satisfy 

the four conditions, we take ^^^1 from which ^^ = 3, and 

11 13 32 

hence from the other three ^, = — , /8j, = -^ , and ^^ = -q- . 

Then 

32^+ 135 + 11(7+92) + 3^ 



s = 



51 



Thus by the operation of indeterminate multipliers we 
reduce in this case the five given equations to three much 
simpler equations whose solution is readily effected by the 
common methods. The values of the other unknown qu^.nti- 
ties can now be either found directly from the normal equa- 
tions by inserting for s and y their values, or by imposing 
new conditions and finding new sets of values for the indeter- 
minate multipliers. r^ t 
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36. Observations oftmequal weight next claim our atten- 
tion. As before, let the observations be represented by the 
equations (66), the measure of precision of the first being A^, 
of the second A,, etc. and their corresponding weights ^j, y,, etc. 
Then the errors being a?j, x^...x^, we may as in (67) write 
them 

aj5 + 5j< + . . . + l^z '-M^ = x^ with weight ff^, 

a^s -h h J ■\-*..'\' I jS'r'M^=oci^Yfith weight ff^. 

Now by the principle of least squares (Art. 16), the most 
probable values of the unknown quantities s, t, z are those 
that make 

(26) h^x^ + h^xl 4- . . . + h^x^ = a minimum. 

But if in this we place 

(57) A/ = A.V». K^'h'9.-K = 'h'3^ 

in which h is the measure of precision of an observation whose 
weight is unity, it becomes 

(75) h^ia^x^ + g^^ -{-g^^ + ... +5^»0 = a minimum, 
so that we have to render a minimum the quantity 

(76) ^gx''=g,x,' + g,x,' + ... +g^x^\ 

Remembering that x^, x^ are functions of the variables 
Sji...z Bs given by equations (74), we must to determine th^ 
minimum differentiate '^ga? with reference to each of those 
variables and place the several differential coeflficients equal 
to zero; thus, after dividing by 2, 

dx. . dx^ , , dx^ ^ 

,^^-v dx. . dx^ , , dx^ ^ 

(77) ff^'^.-df + ff^.-^+"-+9nO'n'-Jt=0> 



dx. . dx^ , , dx ^ 
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Next dififerentiating (74) with reference to each variable, 
we have the coefficients -r-* = a^, -r-^ = ^g, etc. exactly as in 
(69), and, inserting these in (77). we have the conditions 
9x^x^1 + 9J^^^ + .t: + 5^ A^« = 0, 
(7B) gh^^ + gi>^^ + . . . + ^A^H = ^> 



which will be as many as there are unknown quantities 
B,t,..z, If in these we place for x.,x^,,,x^ their values 
from (74), we have the final normal equations which de- 
termine the most probable values of the unknown quantities. 
As in Art. 33, we abbreviate the expressions of these equa- 
tions by placing 

(79) 2gab =ff,aj>, +^M + ••• +5'A*»» 
^gaM^g^a^M, +gfi^M^ + ... +g,%M^, 

etc. etc. etc. 

And thus have the normal equations 

Xga'.8-\-^gab,t+ ... +'Zgal,z = 2gaM, 

(80) Xgab.8 + XgV .t+ ... +Xgbl.jg = XgbM, 



l^gaLs'h'^gbl.t + ... -{-^gP, z = XglM^ 

by whose solution we find the values of «, ^ ... «, 

If all our observations are of equal weight, we may place 
g = l, »nd then our equations (77), (78), (79) and (80) 
reduce to (68), (70), (72) and (73), as determined for that 
case by another method. 

37. The notation above exhibited is very useful as a 
guide in dealing with large numbers of observations, aii(} 
in reducing the numerical operations to a routine for com- 
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puters. To illustrate its uses, suppose we have five ob- 
servations giving the following equations : 

2« - 2^ + 3w = 5 with weight 3, 

«+2*+ 1*= 9 2, 

3^- t+ w= 12 5, 

it-Su- 5 1, 

-«- ^+2m = -7 4, 

from which we wish to find the most probable values of 
s, t and lb. 

Comparing these with (74), we have 

a,= 2, i, = -2, c,= 3, M^^ 5, <7, = 3, 

a,= l, i,= 2, c,= 1, i¥,== 9, 5r, = 2, 

«. = 3> i8=-l' ^»= 1> ^^= 12, 5^3 = 5, 

a, = 0, t,= 4, c, = -3, J/,= 5, (7, = 1, 

a,= l, J, = -l, c.- 2, if,= -7, 5r, = 4. 

Then we have 

a^' = 2, a^Jj = - 4, i^c^ = - 6, ajb^ = 0, etc., 
and hence from (79) we form the sums 

^ga^^gfl^ +g^a,' -rgji^ +g,< +gfi,' = 63, 

Igab =i7A&i+5^M +9z<^A '^9.<^a\ +g,^A = - 27, 
Igac = 5r^a,c^ + g^a^c^ + fir^v, + ^r.a/, +g,a,c, = 43, 
2<7ailf = 200, SjjrJ' = 45, %Jc = - 39, 
2;9'Jif=-6, 27C» = 59, 2^cJf« 62, 
and inserting these in (80), we have the normal equations; 
63«-27^ + 43i* = 200, 
-27« + 45«-39m = -6, 
435-39^ + 59u= 62, 
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whose solution will fumish the most probable values of s, t 
and u. For Gauss' method of still further employing the 
notation in the solution of the normal equations, see Art. 60 
of the Appendix. 

The student will observe that the above process is iden- 
tical with that of multiplying each observation equation by 
the square root of its. weight, and then forming the normal 
equations by the method used in Part I. 

38. To determine the probable errors of the determined 
quantities s, t... z, we let G^, (?, .-. G, denote their weights, 
and Eg., Eg,.. U^ their probable errors. Then if r be the 
probable error of an observation whose vmght is unity y we 
have 

(53) (?.;«.:!:: -1:^,:1, 

from which we find 

(81) R. = ^, 5.= ^, etc. 

Hence in order to find the probable errors of «, ^ . . . «, we 
have only to determine the probable error of an observation 
of the weight 1, and the weights G^, G^,,. 0^. And in 
general, if r^ be the probable error of a measurement whose 
weight is g^y r^ is equal to r divided by the square root of 
.7, (Art. 30). 

To find the probable error of an observation whose 
weight is unity we give the following reasoning. 

Suppose that we have formed and solved the normal 
equations (80), and found the most probable values of the 
unknown quantities. Let those most probable values be 
represented by s,t ... z, and the corresponding true values 
by « + &, t-^Bt... z + 8z, in which S5, Bt ... Bz are small 
unknown corrections. Now if in (74) we substitute the 
values Syt...Zj they will not reduce to zero, but leave the 
residuals v^, n^...v^. Thus 
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a^s + &i« + . . . + l^z - Jfj = Vj with weight g^, 

(82) a,s+i,^+... + i^-3f, = Va gtf 

a^s + 6,< + . . . + Iji - M^ = v^ with weight (7^, 

while if the corresponding true values are inserted, we shall 
have the errors x^, x^..,x^. Thus 

a, (8 + ha) + ii (« + SO + • • • - ^1 = ^1 > 

(83) . a,(« + &) + J,(« + SO + ...-il/, = aj„ 



a.(8 + &)+6„(< + 80 + ...-J*^n = ^n- 

Now each one of the latter equations may be written in 
the form 

a8-\'lt 4- ... + ?« — -Jf, 

+ ahs + 68^ + ... + ZS^ = a?, 

and as (82) are of the form 

a«+ W + ... + Zi^— if = V, 
the equations (83) will reduce to 

Vj + a^8 + h^t + . . . + l^z = ajj weight g^, 
(84) r, + a,8« + i,S^ + ... + l^hz^x^ flr„ 



v» + «n^« + K^t + . . . + /«?« = a?„ weight g^. 

Now the principle of least squares requires that the 
expression 

(76) ff^^' + g,^^+'- + ffn< = ^9^ 

shall be made a minimum to give the most probable values 
of 8, t,..Zf and by the solution of the normal equations 
we have for its minimum value the sum S^rt;*. From the 
residual equations (84), we may find a relation connecting 
the two sums S^rv' and ^ga? by squaring both members of 
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each of those equations, multiplying each by its corresponding 
weight, and then adding the results. Without actually per- 
forming these operations, we see that if the squares and 
products of Ssy Bt ...Bz be neglected as small in comparison 
with SsySt.,. Sz, the result will be of the form 

(85) Sflrv* + k^Ss + kJ5t + ...-]- k^U = Xgsi?, 

in which A? , k^,,,kq are coefficients of the unknown correc- 
tions, and dependent only upon the known constants a^, ij, 
etc. If the number of unknown quantities is j there will 
be q of these tenns. Placing 

kjSs = k^\ k^U='kt^ etc., 

we write the relation 

(85) 25rt;« + Ar/+A;,'+...+A'.» = 2(7af». 

Now the probability of the occurrence of the error x^ 
whose measure of precision is A^, and whose weight is g^^ 
is by (19) and (57), 

(58) y,-=hJg,f,ir'ie-^'^^\ 

in which h is the measure of precision of an observation of 
the weight 1. And hence by exactly the same reasoning as 
in Art. 30 we may show that when w is a large number, 

(61) S^^^i'- 

Further, if we suppose that all the q unknown quantities 
Sf t ... z except one are zero, the equation (85) being true for 
any number q, will hold good for g = 1, and will contain then 

only one k^, whose value as shown in Art. 30 is ^ . Hence 

since each of these A* is a constant, the value of each must 

be ^ whatever be their number. Hence the equation 

(85) is 
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and since there are q terms whose value is ^ . ^ , it is 

from which we find 

Therefore from the constant relation between h and r 
(28) the probable error of an observation of the weight 
unity is 



(87) r= 0-6745./^. 

V n — q 



In the above we have followed closely the reasoning of 
Art. 30. The student can also readily apply the method of 
Art. 24 to produce the same result, by finding from (58) the 
product P of the simultaneous existence of the errors 
x^y x^ ... x^y and then determining h by dififerentiation. 

If in (87) we make g = 1, or consider only one measured 
quantity, it reduces to formula (65), and if we make ^r = 1, 
or regard the observations as of equal weight, it becomes (42). 
Thus measurements upon one quantity are but a special case 
of the more general one of indirect observations. 

39. The probable errors of the values of Syt ... z can 
now be found from (81) as soon as the weights ©,, (?,, etc., 
are known. We now proceed to determine these. 

The observations (74) upon the quantities Jf^, J!f«... J/^ 
furnish the normal equations (80). The solution of these 
equations gives the values oi s^t ... zin terms of -Jf^, M^... 
M^ and coefficients independent of those quantities. Suppose 
the general solution to give 

(88) t = T^M^ + T^M^ + T^M^ + ... + T^M^y 



' = ^^' + ^^« + ^^. + -ot£*oogle 
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in which the coefficients o-, t ... ? depend only upon the con- 
stants a, J ... Z and the weights g in (80). Then if J?, is the 
probable error of a and T^yr^...r^ the probable errors of 
M^yM ... M^y we have, by a principle to be proved in Art 42, 
since the measurements are independent, 

(102) iZ/ = o-,V,« + <r>,' + . . . + ay: = 2<r»r«. 

Now O^ being the weight of s, and ffi, fff-ffn *^® weights of 
the measurements J/^, M^...M^, we have from (53) and 
(54) 

and by substitution in equation 

(90) ij.*=0.i?;(^'+^V...+^') = G.i?.'2^, 

from which we find 

(91) ^- = -T- 

in which the usual notation for 6ums is followed. In like 
manner we may show that the weight of t^ is the reciprocal 

of 2 — , and that the weight of ^ is the reciprocal of 

2-^. 
9 

Owing however to the labour of finding the coefficients 
cr, T...f it is better to deduce these expressions under a 
different form. Let us suppose the normal equations (80) 
to be solved, giving 

8 = a^XgaM-\- a^^ghM-\- ... + a^^glM, 

(92) t^^^lgaM+^^'S.gbM'h... +fi,tglM, 



-•\lgaM-^ \lgbM+ ... +\%^-^GooQle 
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in which oe, /9 . . . X are coefficients independent oi M^,M^.., 
M^. Then the respective weights of 8, t ... z will be 

- , 3- ... — (« being the last unknown quantity and j the 

number of unknown quantities). In order to prove this let 
us find the weight of t By comparison of (88) and (92) we 
have 

-^i = Px9x «i + Pt9j>x + • • • + P,9xh^ 
(93) T, = A5r,a, + /9,5r.6, + ... + p,g,l,. 



Squaring each of these equations^ dividing each by its g and 
adding the results, we have 

(94) +/8.(/3,2flra6+ fi,l.gV+ ... +fi,l.gbl) + .., 

+ 13, {fi.lgal + fi.Xgbl + ... + /9.2^P). 

Now if we were to solve the normal equations (80) by the 
method of indeterminate multipliers (Art 35), we might 
multiply the first by a number 4,, the second by /8,, the q^ 
by /8^, and add the resulte ; then if upon these multipliers we 
impose the conditions 

fi^lga^-^-fi.tgab + ... + I3,lg(d = 0, 

(95) l3,lgab+fi,Xgb''h...-\'fi,Xgbl^l, 

^^Xgal + ^,l.gbl+ ... ^-fi.tgV = 0, 

all the terms except those involving t will reduce to zero, 
and the value of t will be the same as given by the second of 
equations (92). * Comparing then (94) and (95), we have 

(96) 2^ = ^,. or^ = | = (?.. 

9 
Hence the weight of < is ^ , which was to he proved. 
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Although we have here supposed the sohitlon to be made 
by the method of iudeterminate multiplers, it is evident that 
the same result will be reached whatever be the method 
employed. Therefore to find the weights of the values of 
the unknown quantities, we have only to solve the normal 
equations preserving the absolute terms in literal form; then 
the weight of 8 is the reciprocal of the coefficient of the 
absolute term in the normal equation for s occurring in the 
general value of s, the weight of t is the reciprocal of the 
corresponding absolute term in the general value of t, and 
so on. 

40. Thus in the example of Ai-t. 35 we have found the 

17 

general values of y and s. The weight of y is then -=- , the 

reciprocal of the coefficient of E in the value of y ; and the 

. 51 . 

weight of 5 is ^ , the reciprocal of the coefficient of A in the 

value of 3. 

If there be but one unknown quantity, the observation 
equations will be 

and if the weights of these be ^r^, jr, ... gr^, the single normal 
equation is Xga* . z = ^gaM. If a^^a^ = 1, the observations 

are made directly upon JIf, and the equation gives z = -^ — , 

which agrees with the general mean (50). By Art. 39 the 

* 1 

weight of z must be the reciprocal of =- , or 

as shown otherwise in Art. 27. Further, if g = 1 the normal 
equation becomes nz = 23/., agreeing with the law of the 
arithmetical mean (27). Thus is the correctness of our 
methods verified by these mutual checks. 



41. The probable errors of functions of independently 
observed quantities will next be investigated. Let us take 
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first the most simple function of two independently measured 
quantities, viz. 

in which z^ and z^ are the measured quantities ; then it is 
evident that if z^ and z^ are the most probable values found, 
the most probable value of Z is the sum or difference of z^ 
and jz^f as the case may be. Let the errors arising in the 
measurements be 

for^^; a-/, .r/', a?/", etc, ; 
for z^ ; w,\ x^\ x^'\ etc. ; 
then the errors of Z will be 

X^ = ^i' ± «^./, -3^2 = «?/' ± «j", X3 = w^" ± x^\ etc. 
Squaring and adding these errors, we have 

XI + Z/ + eta = (aj/ ± </ + « ± <T + etc., 
or developing and adding 

(97) SJT = 2x,» + 22a:,aj, + Xxl 

In a great number of observations there will probably 
be as many positive as negative products of the form x^, x^, 
and hence we shall have Xx^x^ = 0. Hence (97) becomes 

(98) XX^ = Xx^'+%x^\ 

Denoting the measures of precision of Z, z^ and z^ by 
H, Aj and A^, this becomes by (34) 

' Denoting the probable errors by R, r^ and r^, we have 
(28) iriJ = V, = V2 = ^'^769, 

and hence this relation in connection with (98) gives us 

(99) J? = r,' + r,». 
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In like manner, if we have a function of the sum or 
diflferenoe of several independent observed quantities, viz. 

the probable error of Z will be given by the relation 

(100) B^^T^ + r/ + r,« + ... + r,«. 

42. Next let Z be a function connected with an ob- 
served quantity z^y by the constant relation 

Then if the probable error of z^ is rj, it is evident, since an 
error x^ in the measurement of z^ gives to Z an error Ax^, 
that the probable error of Z is 

(101) R^Ar^, or JK» = ^\«. 

Hence combining this with the principle deduced in 
the previous Article, if 

Z= Az^ + Bz^ + Cz^ + etc. 

and if z^, z^, etc. are independently observed with the pro- 
bable errors r^, r^, r^, etc., the probable error of Zis given by 

(102) R = A\^ + Rr^ + (7V3*+ etc. 

Thus if ^1, ^2, 2?3 ... z,^ are several observed values of the 
sarrw quantity, the probable error of their sum is 

-2= n/^' +< + ...+ rj=^jm^, 
and by (101) the probable error of -th of this sum is 
_ Jnr' _ r 

which is the probable error of the arithmetical mean, as 
has been otherwise shown in Art. 26. 
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Conditioned Observations. 

43. In all that precedes we have supposed that the 
quantities to be determined by observation were independent 
of each other. Although they have been related to each 
other through the observation equations, and have been 
required to satisfy approximately those equations, they have 
been so far independent, that any one unknown quantity 
might be supposed to vary without affecting the values of 
the others. The methods above developed show how, out 
of the many equal possible systems of values, we can deter- 
mine the most probable. 

We now come to a second class of observations in which 
all systems of values are not equally possible owing to the 
existence of conditions which must be exactly satisfied. 
Thus having measured two angles of a triangle, the adjusted 
value of one is entirely independent of that of the other, 
but if the third angle be measured, the three angles are 
subject to the rigorous geometrical condition that their 
sum must be exactly 180^ We have then in conditioned 
observations two classes of equations, observation equations 
and conditional equations, the number of the first being 
generally greater than the number of unknown quantities, 
and that of the latter always less. 

44. The number of observation equations we designate 
as before by n, the number of unknown quantities by q, 
and the number of conditional equations by p. If no con- 
ditional equations existed, the principle of least squares 
(Art. 14) would require that the adjusted system of values 
should be the most probable for the n independent observa- 
tion equations. But here these n equations are conditioned 
by p conditional equations. The p conditional equations 
being less in number than the q unknown quantities, may be 
satisfied in various ways, and further, the final adjusted 
system of values must exactly satisfy them. Hence we 
conclude that of all ilie systems of values which epcactly satisfy 
the p conditional equations, that one is to be chosen ck the best 
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which in the n ohservation equations makes the sum of the 
squares of the residuals a minimum. 

We may then reduce the problem of conditioned obser- 
vations to that of independent ones, by finding from the 
p conditional equations the values of q unknown quantities 
in terms of the remaining q — p quantities and substituting 
them in the n observation equations. There will thus result 
n observation equations, containing, however, only q — p in- 
stead of q unknown quantities, and each of these equations 
will represent an independent observation. From these 
equations we proceed to form the normal equations (73), 
whose solution will give us the most probable values of the 
q — p unknown quantities. Substituting these values in the 
p conditional equations, we find the values of the remaining 
q unknown quantities. Thus the system of values will 
exactly satisfy the conditional equations, and at the same 
time be the most probable system for the observation equa- 
tions. This, therefore, is a general solution of the problem, 

45. Although this is perfectly general and simple in 
theory, it gives rise in practice to tedious computations, for 
we have &:st from the p conditional equations to eliminate 
q unknown quantities, and again solve the normal equations 
to find the values of the remaining q — p quantities. The 
process generally used by computers is Gauss' Method of 
Correlatives, which in Art. 45 of Part I. is fully illustrated by 
examples of conditioned observations of equal weight, and 
which in Art. 48 we shall proceed to develope for the more 
general case of unequal weights. To apply the proof to 
those of equal precision or weight, we have only to place 
^ = 1, that is, omit g from the formulae. 

47. The method of Art. 44 is perfectly general, whether 
the observations be of equal or unequal weight. We have 
in the latter case to find from the p conditional equations 
the values of any q unknown quantities in terms of the 
remaining j— p quantities, and substitute their values in 
the n observation equations (74), each of which is then in- 
dependent ; then, applying the weights we form the normal 
equations (80), whose solution gives us the best system of 
values for the j — p quantities. The remaining p quantities 
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are then directly found firom the p conditional ecfnations. 
The method of correlative© is however greatly to be preferred 
for simplicity in the numerical operations ; and this we now 
proceed to develope for the ordinary case of linear con- 
ditional equations, to which all others may always be 
reduced. 

48,* Let n observations be made to determine the 
values of q unknown quantities, which axe subject to p 
rigorous conditions. Whether the measurements be direct 
or indirect, let them first be supposed independent, and let 
them be adjusted by the methods of Chapters il. and in., 
and let the resulting values be fif, T ... Z, having the weights 
9gf 9i '" 9z- 1^6* *he most probable system of values which 
we are to find be «, ^...a, and the conditional equations 
which they are to exactly satisfy be 

fltj^ + a^j* + ... + a^« = iVj, 

(103) ^,s + i8,«+_+i8^^=iV„ 



Xj5 + X2^+ ... +X,2?=JVp, 

Then if the values S, T ...Z exactly satisfy these equations 
no further adjustment is necessaiy. If not, let 8\i ... z be 
a system of corrections which applied to S, T ,,. Z will make 
them equal to «, t ... z^ so that 

8^S-V»\ t = T-\-t\..z=^Z-\-z. 

Then substituting these in (103) they reduce to 

(104) 

a,8' + a,t' + ... + a,z' ==N,^{a,8 + aj+ ...) = N\ 

* Essentially the demonstration given by Chauvenet in his excellent 
Treatise on Least Squares, the Appendix to his Spherical and Practical 
Astronomy, Philadelphia, ld67» 
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in which 8*, t\.. z^ are the unknown quantities, and N*, N'\ 
etc. constants depending upon the theoretical constants 
N^, N^, etc. and the observed values 8, T... Z. The number 
of these equations is p. 

Now the n approximate observation equations are 

which we may write 

(105) »-iS=5'=0, f-r«^=:0...^-Z=«' = 0, 

whose respective weights are g^ffff^fft] ^^^ the values 
which we are to determine for 8yi ,..z must not only 
exactly satisfy the conditions (104), but also be the most 
probable set of values for (105). Since «', i ,.. z' are them- 
selves the residuals in the observation equations, this latter 
requirement is by our fundamental principle (26) satisfied 
when the quantity 

(76) 2^v* = g,8'^ + 5^^^ + . . . + ^,«'* = a minimum. 

Putting then the differential of this quantity equal to zero, 
we have 

(106) g.sds' +g,fdt' + . . . + g.z'dz = 0. 

Now if «', { ...z were independent of each other, the differ- 
ential coefficient of Igv^ (76) with respect to each of the 
variables would necessarily be zero (as in Art. 36) ; and then 
8,t\..z being each zero, the most probable values of 
s, t ... z would be their adjusted or observed values S, T ...Z. 
But this expression (106) is conditioned by the equations 
(104), and no values oi 8,1 ...z can be admitted which do not 
exactly satisfy those equations. If then we differentiate 
(104) we have the equations 

a^d8 + a^di + ... + a^dz = 0, 

(107) ^,ds' + ^,dt' + ... +i8,rf/ = 0, 



\ds' + \dt' + . . . + \dz = 0, 
with which (106) must agree and coexist. 
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The number of the equations (107) is p, the number of 
the difierentials ds\ dt' ... dz' is q, and since from the nature 
of the case q is greater than p, we can from (107) find the 
values of p differentials in terms of the remaining q — p 
differentials. Let us suppose this elimination to be per- 
formed, an^ that the values of p differentials found in terms 
of the others are then substituted in (106) ; we shall thus 
have an equation in which the remaining q — p differentials 
will be independent, and the coeflScients of each of these 
q — p differentials will therefore be severally equal to zero. 
Without actually performing this process, in each particular 
ca«e we can arrive at the general result of such an elimina- 
tion and substitution as follows. Multiply the first equation 
of (107) by the indeterminate number K^, the second by 
K^y ... the last by iTp, and also the equation (106) by — 1, 
and form the sum of these products. Then if K^, K^ ... A7 
be determined so that p differentials shall disappear {Ait. 35), 
the final equation will contain only the remaining q — p 
differentials. But these being independent, their coeflScients 
will be severally equal to zero; and hence we have in all 
q conditions that the coefficients d8\ dt' .., dz' in that sum 
shall be each equal to zero, viz, 

a,K^^P,K, + ... +X,iirp-5r/ = 0, 
(108) a^, + ^^^,+ ...+\^p-5r/=6, 



a^K, + ^,K^ + ... + -K^p-g.z - 0. 



If now we multiply the first of these by -^ , the second 

9, 

by — , ... the last by -^, and add the products, we have by 
comparison with the first equation of (104), 

9 ' 9 
in which we have the usual notation for sums, or 
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9 9. 9t 9. 




^«0^«A.^. a 



17? 



etc. 



In like manner multiplying the first by — , the second 

8 

by —f etc, we form a second normal equation. Thus we 

have in all^ normal equations containing the jp new auxiliary 
unknown quantities, viz. 

9 ' 9 * 9 ' 

(109) 2'!^.ir,+S^.ir, + M.+s^ir, = iV" 



The solution of these equations will give the values of 
K^, K^...K^, which being substituted in the correlative 
equations (108) will furnish the values of the required cor- 
rections s\ ^' ... «' ; thus 

(110) «' = i(a^, + /S,2r, + ... + \ff-,), 

etc. etc, 

and these values will not only exactly satisfy the conditional 
equations (104) but will make the sum ^gv^ a minimum, and 
hence are the best attainable system of values. Adding 
these corrections to the observed values iS», l\,,Zy we have 
finally the adjusted values 8,t...Zy and these will exactly 
satisfy the conditional equations (103). 
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49. The Prohahle Errors of conditioned observations 
follow directly from the formulae of Art. 38, and the general 
solution of Art. 44. 

Since the n observation equations contain ^ unknown 
quantities, and by elimination from the p conditional equa- 
tions we reduce that number to q—p independent unknown 
quantities, w^ have only in the formula (87) to replace q by 
q—p. Therefore the probable error of an observation of the 
weight unity is 



(111) r = 0-6745. /-^ 

V w — 



^g^ 



and the probable errors of observations or values whose 
weights are g^, G^, etc. are 

(81) ,.-^. B,=^,et. 

If we have determined the probable values of 8,t...z 
by means of corrections 8',t\..z' to the observed values 
/S, 7 ... Z, the sum S^'ti' is simply 

^g^^g.8''+g,i!^+...+g,z'*, 

in which g„gf*g» are given by the observations themselves 
or by adjustment as in Arts. 28 and 39. The final weights 
(?,, Q^, etc. of the values «, t, etc. are found by exactly the 
same process as developed in Art. 39. An example illus- 
trating the operation in full is given in Part L 



The Discussion of Ohservations. 

50. The above methods constitute the whole science 
of Least Squares as applied to observations involving linear 
equations, to .which all others may be reduced (Art. 59) ; and 
the most common formulae and methods for investigating 
probable errors. The determination of probable errors of 
functions of dependent quantities has not been given, as 
such investigations are rarely needed in practical discussions, 
and would be out of place in an elementary text-book. A 
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few points in Chapter V* require perhaps further notice and 
explanation^ 

61. The Deduction of Empirical FormuUB by the discus- 
sion of observations, is one of the most extensive applications 
of the foregoing methods. Having given the numerical 
results of a series of physical observations, we have only to 
assume a general algebraic expression, which includes the 
law of the phenomena in terms of the observed quantities 
and undetermined constants. Then inserting the numerical 
results, we have a series of observation equations from which 
we deduce the normal equations, whose solution gives the 
most probable values of the constants. The choosing of the 
proper algebraic expression is usually the most difficult part 
of this process. For this, no general rules can be given. 
The best plan is to assume convenient horizontal and vertical 
imits and plot the results of the observations, thus obtaining 
a curve which represents them graphically to the eye. A 
comparison of this curve with similar curves whose equations 
are known, will then often enable us to determine the 
general form of a convenient algebraic expression. 

62. If the plotted curve resembles a parabola or hyper- 
bola, it may be represented by the equation 

(112) y=zA-\-Bx-¥Cx*-\- etc. 

in which the absolute term A may often be directly de- 
termined by choosing a proper origin for the values of y 
ando?. 

53. If however the plotted curve repeats itself like the 
curve of signs, the general equation 

(113) y^A^B^mi~x + B^ + Osin ^^2aj+C') 



•f etc. 



4 



will be applicable. Here also the constants A, B, G' may 
often be omitted by choosing a proper point as the origin 
of the co-ordinates x and y. The value of m is generally to 
be assumed from the inspection of the plotted curves, or 
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its probable value be found by successive approximations. 
If this formula be expanded, we have, considering only the 
terms involving A and B, 

^ . T> . 360'* jy , o 360* ' jy , . 
y^A + Bsm xcosB + B cos x sm B + etc., 

and if in this we place 

BcoaB^By JBsinjB'=-B,, 

it becomes 

f^^4^ ^ . D • 360° . o 360' , . 
(114) « = A + 5, sm x + B^ cos x + etc. 

which is a more convenient form for computation. Inserting 
in this the values of y and x from the observations, we form 
the normal equations, and deduce the probable values of 
Ay B^, and B^, Then by means of the above relations which 
furnish 

B, , . TV B^ 



B=-Jis,'^B^, cos^'^:^, and sinB'=^, 

the derived equation (114) can, if desired, be reduced to the 
form (113), which is often more convenient for subsequent 
discussion. 

54. Formula involving undetermined cpnstants like thp 
case of the pendulum given in Part I. occasionally arise in 
theoretical investigations ; and if observations enough exist, 
the constants may be deduced. 

The determination of the probable errors of such formulae 
is rarely necessary, as the comparison of the computed and 
observed results indicate their precision suflSciently well to 
enable us. to decide upon the degree of confidence to which 
they are entitled. The weights iand probable errors of the 
deduced constants can in all cases be found by the methods 
of Arts. 38 and 39. The probable errors of the results de- 
duced from such formula (for example, the probable errors 
of the values found for y in (113) after A, B, etc. have been 
determined) cannot however be found from the relation given 
in (102) because the separate terms are not independent For 
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methods applicable to such cases we must refer the reader 
to the larger and more complete treatises upon the subject, 
a list of which is given in Art. 64. 

55. The discussion of the probability of errors or of the 
accuracy of observations is of importance in delicate mea- 
surements. This we have hitherto done by means of the 
Probable Error, or the error such that is an even wager that 
the result is within that amount of the truth (Art. 16). It 
is perhaps unfortunate that this particular error has been 
chosen as the one for comparison, for the mind is better 
satisfied with considering an error such that the probability 

99 
of an error being less than it is =-7^7T or some higher fraction 

instead of ^ • Such comparisons are readily made by the 

Table given in Part I., and we have only to explain the 
manner in which it is calculated. In Art. 12 we have shewn 
that the expression 

(22) F^^Te-^'^d.hx, 

expresses the probability that an error will be included 
between the limits — x and -fa?. If in this we place 

, 0-4769a? ^ 

hx^ = t, 

r 

we may write it 

(115) F^^^jV^dt, 

and by the methods explained in the foot-note to equation 
(22), its value may be found for successive numerical values 
of t = hx. But if we wish instead of h to employ r, we can 

compute it for successive values of - = ^ .^^^ , that is, in 

the Table in Art. 13, Part I, we have only to divide the 
numbers in the column hx by 0*4769 in order^o reduce 
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them to the values - , Then by interpolation the Table is 
easily written as in Art. 65. 

If then we hav^ made n observations which give a mean 
z^ with a probable error r^ we can easily find the probability 
that Zq is within ±x oi the truth by taking from the table 

the fraction corresponding to — . And conversely, if we 

ask what is the error a?, such that it is a wager of 99 to 1 
that z^ is comprised within the limits — a? and «H-a?, we 

have only to take the number -- correspondiDg to P' = 0'99. 

56. The table may also be used to investigate the pro- 
bability of constant errors, and to discuss numerous questions 
arising in the study of statistics, into which, Lowever, the 
plan of our book forbids us to enter. 
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58. The elementary applications and the theory of Least 
Squares has now been given and exemplified. A few other 
applications and extensions valuable to the computer, and 
a brief notice of the history and literature of the subject, 
interesting to all who have studied the science, will next be 
presented. 



Observations involving non-linear Equations. . 

59. In all that precedes, we have supposed that the 
observations can be represented by equations of the first 
degree : if this is not the case, but higher equations are 
involved, they can readily be reduced to linear ones by the 
following method. 

Let the quantities to be determined be represented by 
5, t... z, and the measured quantities by M^, M^... M^ and 
the observation equations have the general forms, 

f^{s,t...z)^M^ 



f^{s,t...z)^M^ 
n being the number of observations. These may be jnritten 

(liey <l>,=^f,{s,t...z)^M,^0, 
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Now let approximate values oS s,t...z he found either 
by trial or by a solution of a sufficient number of these 
equations, and let them be denoted by 8, T.,.Z, and let 
s\ t' ... z' he the most probable system of corrections to these 
values, so that 

Developing then the expressions (116) by Taylor's 
theorem, we have, neglecting the products and higher powers 
of the corrections 8\ t' ...z\ 

*.=/.(s,r,..z)-if,+^,-+^<'+...+§,., 
*.=/.(« r...z)-jK;+ f.-H-f «•+...+*•/. 

(117) 

Designating the constant term f^{S,T,..Z)hjN'^ etc, 
these become 

(118) #.,^#.,^...^^^_,f__^_, . 

etc. etc. 

where -j-^ , -^ , etc. are simply the differential coefficients 

found by differentiating each of the equations (116) with 
reference to each of the variables and then substituting 
8,T...Z for 8,t...z, and are hence constants. Denoting 
them then by a^, h^, etc., we have 

a^s+hJ+.^ + l^z'^M^-N^, 

(119) a,s+h,H'+... + l,z'^M,--N^, 



^"^'■^^»*'-*--^^^' = ^-£k.GoogIe 
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in which all the letters except «', t\..z denote kno wii quanti* 
ties. These equations are exactly like those of (66) or (74), 
and from them we form the normal equations, whose solution 
gives us the most probable values of the corrections s\ t\,. z\ 
and hence the best system of values for the observed quantities 

If non-linear conditional equations are also given, wo have 
only to find approximate values for the unknown quantities, 
and assume a system of correctipns. Then the functional 
conditional equations may be developed as above by Tatt.oe'3 
theorem, and reduced to linear equations of the same form 
as (104), which may be treated by the method of correlatives 
and the most probable system of corrections determined, 
which applied to the approximate values will give the 
adjusted results. If these do not satisfy the original condi- 
tional equations with sufficient accuracy, a new system of 
corrections may be assumed and the process again reptsatcd. 

In Art. 46' is exhibited the reduction of a transcendental 
conditional equation to a linear one, by the use of the tabular 
logarithmic differences, which is more convenient than the 
treatment by Taylor's theorem. The latter, however, must 
be used for higher algebraic or exponential equations. 



Oauss^ Method for the Solution of Normal Equations. 

60. The formation and solution of normal equations is 
the most laborious part of the practical reduction of observa- 
tions. In dealing with large numbers of these equations, 
computers usually follow the method of Gauss, by which the 
work is reduced to a systematic routine. This method con- 
sists in solving the equations by substitution so as to preserve 
throughout the work the symmetry which exists in the 
coefficients of the normal equations. To illustrate it, it will 
be sufficient to consider a case involving but three unknown 
quantities arising from observations of equal weight. Let 
the n observation equations be 



(66) a,s + h,t + c,u^M^, 

a^8 + Iji + cjii = M^, etc. 
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The three normal equations formed from these will be 
Xa*.8 + SoJ . t + Sac . w = 2a3f, 

(73) 2a6.« + 26\«+S6c.ti = 26if, 

'Xac .8 + Xbc , < + 2d* . u = ^cM, 
in which 

(72) S6*=5,« + V+... + V, 

2Jc = hjCj^ + ftjC, + ... + J»c« etc. 

The coeflScients of the unknown quantities are symmetri- 
cal, the first horizontal and vertical rows being alike, the 
normal equation for s being distinguished by the presence of 
2a', that for t by 26* and that for u by 2c*. Now let us 
proceed to find the value of u from these equations. The 
value of 8 from the first equation is 

/ion\ ^^^ 2aJ . 2ac 

(120) ,«^-^.«-^,.«. 

Placing this value of ^ in the second and third equations, they 
become 

(121) 2,6* . t + tfic . u = 2,61f, 
2i6c.< + 2iC*. M = 2iCi^ 

provided that we place,. 

2,6* = 26*-|SSaJ, 
Xfic = 26c — .=r^ 26c, 

(122) 2,c* = 2c*-|^2a6, 

2,6ilf=26ilf-^26c, 
* 2a' 

2jCm SB 2c Jf - -y-T ^^- 
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The two equations (121) are then exactly Bimilar in form 
to the second and third original normal equations, except 
that the terms containing s have disappeared, and each 
coefficient is marked with the index 1. From the first of 
these we take 

and substitute it in the second, giving 

(124) \<f.u = X,cM, 

in which, as before, we preserve the symmetry of the coeffi- 
cients by writing 

(125) V = V-^2^, 
2,cJlf = 2,clf — ^-j- 2j5c, 

and hence we have for u, 

2.cJf 



(124) « = 



V 



Inserting this in (123) we have the value of t, and then 
from (120) the value of *. 

The expressions (125) and (122) which result from the 
abridged notation are called auxiliaries. To deal with 
equations by this method, we may as Art. 37 form the sums 
(72) and the normal equations (73), then having obtained the 
coefficients 2a*, etc. we insert them in (122) and have the 
first auxiliaries S^a', X^ahy etc., which being substituted in 
(125) give the second auxiliaries S^c*, etc., and lastly, from 
(124), (123) and (120) the values of w, t and 8, By the 
process as thus followed the work is reduced to a routine, 
which may be followed by one ignorant of the theory of the 
method. As a check upon the work, the normal equation 
for 8 may be taken as the final one and a new set of auxiliaries 
written, and the values deduced in the reverse order «, ^, u. 
Other checks upon the accuracy of the work are also afibrded 
by the properties of the normal eejuations, which will readily 



188 APPENDIXi 

occut to all cdmputers. In common calculations, however, 
the use of such routines is not to be recommended. 

61. In Art. S9 we have shown that the weight of any 
unknown quantity as t* is the reciprocal of the absolute term 
in the normal equation for u contained in the general value 
of i^. Now in solving the normal equations as above, it is 
evident that the absolute term in the normal equation for u 
is unafifected by substitutions of t and 8 from the other equa- . 
tions. Hence in the above example SjC* is the weight of the 
value of u. 

In this manner, the values of the unknown quantities 
and their weights may be found, and thus the results obtained 
by the process of Art. 39 be checked by an independent 
method. In fact, in all calculations where large numbers of 
observations enter, the solutions should be made by separate 
methods, and if possible by independent computers. 



Other Formulas for ProhahU Errors. 

62. The formation of the squares of the residuals involves 
in practice considerable labour, when the residuals are so 
large that tables of squares cannot be used. Formulae have 
hence been deduced for probable errors in which only the 
residuals themselves are employed. We give here, without 
demonstration, two such formulae for the common case of 
direct observations of equal weight (Arts. 23 — 26). Let ii 
denote the number of measurements, and 2v the sum of the 
residuals v^, v^ *.. r^, all being taken with the positive sign^ then 
if n is a large number, the probable error of a single observa- 
tion is, approximately, 

(126) r = 0-84!53 — , 
or more nearly, 

(127) r = 8453 ^^ 



^/r^(n-l)' 



Thus in the example of Art. 24 we have 24 observations, 
and the sum of the residuals aU taken positively is 37'''48. 



r 
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The formula (126) gives r = l"-32 and (127) gives r = 1'344, . r 
while from the stricter formula (42) we found r = 1'349, , 
With a larger value of n, a closer agreement might be - 
expected. The probable error of the arithmetical mean found 
by dividing r by the square root of w, will in the above 
example be practically the same by all three formulae. For 
values of n less than 24 it is best to hold fast to the more 
exact formula (42), and even that cannot for such cases be 
expected to give precise results, since the hypothesis of its 
development supposes that enough observations have been 
taken to exhibit the several errors in proportion to their 
respective probabilities. The formula (127) is due to Peters. . • 



The Mean Error, 

63. The choice of the probable error as a means of com- 
parison of diflferent series of observations is, as we have before 
mentioned (Art, 55), entirely arbitrary, although it seems to 
be the most natural one from its middle position in the series 
of errors (Art. 16). Another error very commonly employed 
for the same purpose is called the mecm error, whose defini- 
tion is, the error whose square is the mean of the squares qfall 
iJie errors. Hence, the mean error is the square root of the 

quantity — (59), and is consequently that part of the pro- 
n 

bable error included under the radical sign. If then e be the 

mean and r the probable error, 

(128) * = a6fe = l-*826n 

To transform then our expressions for probable error into 
those for mean error, we have only to omit the constant 
factor 0*6745. In Fig. 1. the abscissa OP denotes the pro- 
bable and OM the mean error. It is a probability of ^, or a 
wager of 1 to 1 that an error taken at random is less than 

682 
the probable error; it is a probability of :rjr7rx or a wager of 

214 to 1 that it is less than the mean error. 
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In this book we liave chosen to employ only the probable 
error as being the simplest in theory, and the one in most 
common use. For the sake of uniformity, it is certainly to 
be desired that the mean error should be discarded and the 
former only employed. 



List of Literature. 

64. The following list gives the titles of some of the 
most important memoirs and books upon the Method of 
Least Squares and the law of errors of observations. It 
is intended to include not only the best text-books, but also 
those works which are of the greatest historical value. The 
arrangement is chronological 

It would be easy to greatly extend the limits of this 
list. The titles have in fact been selected from a list of 
about four hundred which I hope sometime to publish, ac- 
companied by historical and critical notes. But as an aid 
to the general reader, a selection such as here given will 
prove of greater value than if the number were increased 
tenfold. 

1. Cotes. Estimatio errorum in mixta mathesi ...; 
ffarmonia Mensurarum (Cantabridgiae, 1722, 4to.), pp. 2 — 
22. 

2. BoscoviCH. De Uttera expeditione per Pontificiam 
ditionem ad dimetiendos duos meridianigradvs; Romae, 1755, 
4to., pp. xxii, 516. 

3. Simpson. An Attempt to show the Advantage of 
taking the Mean of a Number of Observations...; Misceir 
laneous Tracts (London, 1757, 4to.), pp. 64 — ^75. 

4. Laplace. Determiner le milieu que Ton doit prendre 
entre trois observations... . Mdva, Acad, Paris par divers 
savans [etrangers], 1774, VoL vi. pp. 634 — 644. 

5. Bernoulli (Daniel). Dijudicatio maxime probabilis 
observationum discrepantium atque verisimillima inductio 
inde formanda; Acta Acad, Petrop. for 1777, Pt. I. pp. 
3—23. 
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6. Laplaob. Chap^ V. Book ni. of Traits de meoomque 
celeste^ Paris, 1799, 4to. 

7. Leqendrk. N(mveV^ m^thodes p(mr la d^ermincUion 
des orbites dea comdtes; Paris, 1805, 4to., pp. viii, 80. 

8. Adriak. Research concerning the probabilities of the 
errors which happen in making observations ; Analy^y 1808, 
No. IV. pp. 93—109. See ATner. Jour. Bd. 1871, VoL i. 
p. 412. 

9. Gauss. Determinatio orbitsB observationibus quotcun- 
que quam proximae satisfacientis ; Theoria motus corporum 
codlestium (Hamburgi, 1809, 4to.), Lib. IL Sect in. pp. 
205—224. 

10. Laplace. Th^orie analytique dea Prohahilit^s ; 
Paris, 1812, 4to., pp. 464. Third ed. with supplements, 
1820, pp. c^ilii, 606, 34, 50, 36. 

11. Gauss. Bestimmung der Genauigkeit der Beobach- 
lungen; Zeitschr.f. Astr., 1816, Vol. L pp. 185 — 197. 

12. Paucker. JHe Anwendtmg der Methode der klein- 
sten Quadratmmmen auf physikalische Beobachtmigen: Mitau, 
1819, 4to., pp. 32. 

13. Gauss. Theoria combinationis observationum 
erroribus minimis obnoxiae ; Comment Soc. Ootting., 1819 — 
22, VoL V. pp. 33—90. 

14. Ivory. On the Method of the Least Squares: 
Phil. Mag., 1825, Vol. LXV. pp. 3—10, 81, 161—168; 1826, 
VoL Lxvin, pp. 161 — 165. 

15. Gauss. Supplementum theoriae combinationis ob- 
servationum... ; Comment Soc. Ootting., 1823 — 27, VoL vi. 
pp. 57—98. 
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On the History of the Method of Least Squares, 

65. In the following brief sketch constant reference is 
understood to be made to the preceding list of literature. 

The average or arithmetical mean has always from the 
earliest times been employed for the combination of direct 
observation* of equal precision made upon a single quantity. 
Out of this arises so naturally the idea of weights and of 
the general mean (Art. 27), that it is probable that both wer6 
in extensive use long before any attempt to deduce general 
rules based upon scientific principles. The first recorded 
discussion of indirect observations with a view of establishing 
a general method of adjustment seems to be that of Cotes, 
who about the year 1714 introduced the idea of observation 
equations, or, as they are often called, equations of condition. 
He took up the simple case of determining a quantity z 
from measurements upon the related quantity 3I=^az, and 
representing each measurement by an equation, obtained 
a^z — if J = 0, a^z — if, = 0, etc. "Without here entering into 
the details of his reasoning, we may say that his method 
consisted essentially in adding the several equations, and 
from their sum finding the value of z, which value he re- 
garded as the most advantageous or plausible. This process 
involves the principle that the algebraic sum of the errors 
shall be zero, and since this coincided with the practice 
of the arithmetical mean, when the constants a^, a,, etc. 
were equal, it seemed to be the proper adjustment of such 
cases. It afforded, however, no means for the combination 
of equations containing more than one unknown quantity. 
The next step appears to have been made by Pater 
BoscoviCH, who being deputed by the Pope to measure an 
arc of the meridian, published his results in 1755, and who 
introduced the plan of adjusting observations by introducing, 
in addition to the principle of Cotes, the condition that the 
sum of the errors all taken positively should be a minimum. 
These two principles rendered possible the adjustment of 
indirect observations involving several unknown quantities, 
and appears to have been considerably employed by subse- 
quent writers, particularly by Laplace, who in 1799 gave 
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an exhaustive discussion of observations upon the leagth 
of the second's pendulum and the resulting ellipticity of the 
earth. This method, although incorrect, served to prepare 
the way for a better one. 

The first published application of the method of Least 
Squares is due to Legendre, who in 1806 stated, withoafc 
demonstration, the principle that the sum of the square^^ 
of the errors should be a minimum, and gave an example 
illustrating its use in determining the orbits of comets. 
The honour of its discovery and introduction is however 
universally conceded to Gauss, who as early as 1795 had 
used the method in his computations, and had communicated 
it to his astronomical friends, and who in his Theona motus 
corporum, published in 1809, gave first its development and 
demonstration. This proof has been followed by the great 
majority of subsequent writers, and being the best adapted 
for an elementary presentation of the subject is also used 
in this book. To Gauss is also due the development of the 
algorithm of the method, the formulae for probable error, 
the determination of weights, the method of correlatives, 
and many other features of the subject, as well as numerous 
practical applications with which his writings abound. Very 
few branches, of science owe so large a proportion of subject 
matter to the labours of one man. In 1812 Laplace in his 
treatise on Probabilities took up the subject, and gave an 
entirely diflFerent demonstration of the fundamental principla 
that the sum of the squares of the errors should be a mini* 
mum, and in 1826 Gauss in his Supplementum, etc* again 
returned to the subject with a second demonstration. 

The method thus thoroughly established spread among 
astronomers with wonderful rapidity. The theory was sub- 
jected to rigid analysis and discussion by IvORY, Pc^tsson, 
Encke, and others, while the labours of Hansen, Bessel, 
and Gerling developed its practical applications to astrono- 
mical and geodetical observations. The works of Eksskl 
may be particularly mentioned as establishing the processes 
since universally employed in extensive trigonometrical 
surveys. Our partial list of literature indicates the thorough- 
ness with which the subject has been treated by German 
and French scientists. In this connection we^-™^ jbe 
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allowed to mention as a guide to the student, that Nos. 14, 
23, 27 and 46 in that list give interesting critical discussions 
of the theory of the subject, and that No. 25 is a popular 
and interesting presentation of the law of probability of 
error. The Method of Least Squares is now universally 
employed in all branches of physical science where ac- 
curate observations arise, and it is perhaps not too much 
to say that the precision of astronomical tables is due in 
great part to its use. 



Remarks on the Theory of Least Squares. 

66. The proof which we have given of the Method of 
Least Squares is the one presented by Gauss in 1809 and 
followed by the great majority of subsequent writers. It 
first establishes the law of the probability of error repre- 
sented by the equation (12) y = ce~'^*'*, from which the 
principle of the method immediately follows. The whole 
reasoning is thus dependent upon the theory of proba- 
bility and upon the particular law of error y^ce'^^. In 
the demonstration of this law of error in Art. 11 J, there are 
two defects, which I now proceed to point out. 

The first is the assumption that the average or arith- 
metical mean furnishes, for direct observations of equal 
weight, the most prohahle value of the quantity sought. In 
the strict mathematical sense of the words "most probable," 
as defined in Art. 9, this is not true, or, if true, has never 
yet been proved. In the common sense of the words the 
arithmetical mean may perhaps be regarded as most correct, 
most advantageous, most plausible, or even as most probable, 
but I am unaware that it has ever been shown that the 
a priori mathematical 'probability of the average is the 
greatest out of all the probabilities of all the assignable 
values. Hence in equation (8) it is not really known that 
the letter z represents the same quantity on both sides of 
the sign of equality. 

The second is the transition from equation (8) to equa- 
tion (9). Granting that in (8) the value of^r is the same 
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in both members, it follows that the corresponding terms 
are equal each to each. But the quantities {z -itfj (z — M^), 
etc. in that equation are not actual errors, since z is not here 
the trtie but the probable value of the measured quantity. 
Hence equation (9), in which a; represents the actual true 
error, does not strictly follow from (8) unless we grant that 
the same law of error which occurs in the particular case (4) 
obtains also in the general case (3), that is, unless the law 
which is true for the residual is true also for the error. 

These two difficulties have puzzled mathematical students 
since the year 1809, and they cannot be bridged over or 
avoided, but will always exist in this mathematical develop- 
ment of the law of probability of error. And I think for 
this reason ; it is distinctly conceivable and hence a priori 
possible that in different .classes of observations different 
laws of error might exist, therefore we ought not to expect 
that from merely theoretical considerations a single definite 
law of error should result. And so we introduce an element 
derived from eooperience, viz. the arithmetical mean, which 
we know in a large number of observations gives a very 
near approximation to the measured quantity, and which 
in an infinite number of observations would give us its true 
value. At the limit then when n is infinite, all the equa- 
tions of Art. 11 are in strict agreement, and one law of 
error obtains. 

The demonstration of Laplace given in his Thdorie 
analytique des ProhahilitSs employs very different reasoning. 
It shows that the sum of the squares of the errors must 
be made a minimum in order to furnish the most ad- 
vantageous (not the most probable) values of the quantities 
sought, without uny reference to the law of probability of 
error. The proof is stated by him and by other writers 
to be entirely general, whate\^er be the law of probability 
of error. This point I do not propose to discuss, but will 
mention as a curious circumstance that Airy (No. 36 of 
our list of literature), who follows that method, finds no 
difficulty in deducing an expression which agrees in every 
respect with our equation (19), y = At7r-*6-'^^. I may say 
also that Laplace's reasoning supposes the number of ob- 
servations to be infinite. This demonstration k followed 
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as far as I am aware by only three authors, viz. by PoissON, 
by De Morgan and by Airy; it is severely attacked by 
Ivory, who proves conclusively, as he says, that it involves 
the particular law of error deduced by Gauss' method, and 
it is defended by Ellis^ who however does not show that 
Ivory was wrong. 

Ivory himself gives no less than three demonstrations 
of the method of least squares, which are entirely indepen- 
dent of any idea of probability, and which although untenable, 
furnish extremely interesting reading and deserve to rank 
high as mathematical curiosities. 

A critical and exhaustive discussion of these and other 
proofs of the Method of Least Squares is given by Glaisher 
in a memoir published in 1872. (See No. 46 of the above 
list of literature.) 

In conclusion, I may say that another difficulty some- 
times found in Gauss' proof does not exist as here pre- 
sented. In the equation (19) y = AiV"* e^^^, the probability 
y is an infinitesimal if a; is a continuous variable, owing to 
the presence of the letter i, and is an abstract number as it 
ought to be (see Art. 12). The law as given by most writers 
is y = A7r~*e"^^, which is absurd, since u is then not only 
a concrete number but finite, and hence the probability of 
an error lying between any two assigned limits must be 
infinity, a conclusion which those authors take care to avoid 
by an ingenious but questionable artifice. 
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„ of a triangle, 73, 78, 88, 99 
„ of a quadrilateral, 76, 92 
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, , probable error of, 26, 147 
Axioms, 9, 124 

Base line, 33, 39, 63, 93 

Bessbl's observations, eta, 13, 195 

Borings at Grenelle, 102 

„ at Speremberg, 104 
BoscoviCH, historical notice of, 194 
Bowditch's formula for pendulum, 
109. 

Certainty, 4, 121 
Chemical analysis adjusted, 93 
Chronometer, rate of, 65 
Cleveland, Ohio, elevation of, 48, 62, 63 
Coast Survey Reports ^ 27, 85, 108 
Coins, throwing of, 3, 5, 6, 123 
Columbus, Ohio, elevation of, 48 
Compass bearings, 26, 31 
Compass, declination, 106 
Conditional observations, 22, 68, 172 
Constant errors, 3, 115, 120 
Correlatives, method of, 73, 88, 173 
Cotes, historical notice of, 194 
Criterion for doubtful observations, 

117, 196 
Curve of Probability, 11, 125 

Declination, magnetic, 106 

Direct observations, 22, 24, 41, 141 

Earth, temperature of, 102 



Empirical formulae, 100, 179 
Equations, conditional, 68, 172 
„ correlative, 73, 88, 173 
normal, 44, 74, 156, 161, 
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„ observation, 41, 153 
„ of prob. curve, 10, 131 
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„ accidental, 3, 120 
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„ see Probable Error 
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Formulae empirical, 100, 179 
Friction, coe£&cient of, 66 
Function, defined, 9 

Gabdner's report on levels, 47 

Gauss' law of errors, 10, 126 
„ method for weights, 59 
,» „ for correlatives, 73, 

88, 173 
„ historical notice of, 195 

General mean, 31, 148 

„ prob. error of, 36, 149 

Gravity at New York, 101 

Guessing, problem on, 116 

Harrisburg, Penn. elevation of, 48 
Hartford, Conn, magnetic declina- 
tion, 107 
History of Least Squares, 194 

Independent observations, 22, 46, 153 
Indeterminate multipUers, 157 
Indirect observations, 22, 46, 153, 
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Integrals determined, 133 n. , 135 n. 
Iron, Cast, coefScient of friction, 66 

„ Pig, analysis of, 93 
Iyobt, on pendulum, 110 

„ on theory of least squares, 198 

Lake Erie, elevation of, 48 

Lai lAci:, liiBtori(?al notice, 194, 197 

Lntitiido, dc't^^rmiDation of, 65 

Law of gravity, 101 

Lttwof probability of errors, 8, 124, 181 
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Leoflt B«iu!i.r^s, critical remarks, 196 
„ history of, 194 

ti method of, 14 

„ principle of, 17, 136 

1, literature of, 190 
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Literature^ Hat of, 190. 
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McEUi, mitLiiiiiji.iLiiil, 17, 25, 141 

„ general, 31, 148 
Mean error, 189 

Measure of Precision, 10, 19, 131, 138 
Meter, length of, 102 n., 110 
Mistakes, 2, 120 
"Most Probable," 1, 7, 123 

Normal Equations, 44, 46, 156, 161, 
185 

Observation defined, 1 

„ equations,41,42, 45, 153 

Observations, adjustment of, 14, 46 
„ comparison of, 17 

„ conditioned, 22, 24, 70, 

172 
direct, 22, 24, 41, 141 
„ discussion of, 100, 177 

„ doubtful, 116 

,, errors of, 2, 120 

„ independent, 22, 153 

indirect, 22, 153, 194 
„ kinds of, 22, 141 
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Pendulum, second's, 109 
Personal equation, 2 
Pittsburgh, Penn. elevation of, 48, 62 
Precision, measure of, 10, 18; 131, 

138 . 
Probable error, 17, 138 

„ of arith. mean, 29, 146, 171 
„ of general mean, 36, 149 
„ of indirect observations, 60, 
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„ of conditioned observations, 

18, 139 
„ of single observations, 26, 

142, 188 
„ of functions, etc. 64, 169 
Probability, defined, 4, 121 

„ curve of, 8, 10, 125 

Probabihty of errors, law of, 7, 124, 
181 

„ table of, 12,112 

Probability, problems on, 7 

Quadrilateral, adjusted, 80, 85, 92 

Besiduals, 27, 142 

SoHOTT, on magnetic declination, 106, 

108 
Sets of observations, 36, 148 

Temperature of Earth, 102 
Triangles, plane, 73, 78, 88, 99 

„ spherical, 92 
Triangulation, 84, 89 

United States Coast Survey Reports, 

27, 85, 108 
United States and English feet, 110 n. 

Water, volume and temperature of, 

105 
Weights, 30, 148 

Weights of averages of sets, 25, 148 
„ of indirect observations, 60, 
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